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Dear Dr. Talay, 

Attached are the final report materials for NASA Cooperative Agreement NCC 1-229 
entitled Improving Conceptual Design for Launch Vehicles” that was conducted by the Space 
Systems Design Laboratory (SSDL) at the Georgia Institute of Technology during the period 
September 23, 1996 to September 22, 1999. The most recent year of this cooperative agreement 
focused on the economic viability of a Bimese reusable launch vehicle when servicing future 
intercontinental fast package delivery missions. Several alternative configurations were considered 
and the results are summarized in the attached report. In addition, a summary of funded activities 
and a set of presentation-style charts from this project are also included. 

This report concludes the nominal period of performance of the three-year NCC 1-229 
activity. During this overall project, the SSDL students and researchers developed a new cost and 
business assessment model for estimating the economics of reusable launch systems. This tool is 
called CABAM, and it has since been made available to a number of advanced design 
organizations. In year 2 of this project, our efforts focused on computational frameworks for 
supporting multidisciplinary design environments. For this task, we created a web-based 
demonstration environment that integrated four disciplines for designing an RBCC SSTO 
spaceplane. Our most recent activities have focused on the assessment of a Bimese RLV for 
intercontinental fast-package delivery markets. 

I would like to thank you for your continued support of our research and educational 
activities in advanced space transportation system design. Partly with the support of this 
cooperative agreement, SSDL has grown from only 2 graduate students in 1996 to a more fully 
developed and capable research lab of 18 graduate students and 4 undergraduate students in 1999. 
We look forward to working with you and your organization on future projects. 
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Dr. John R. Olds ^ 

School of Aerospace Engineering 
Director, Space Systems Design Lab 
Georgia Institute of Technology 
Atlanta, GA 30332-0150 
404-894-6289 
john.olds@ae.gatech.edu 



A Umi of tht' University Suttcm ofGcorgiii 


An Lijuol LdUiiitiO}i toiil l ni}’lo\ttih’ii! i^'portunitv 




Final Summary of Cooperative Agreement Activities 
Improving Conceptual Design for Launch Vehicles 

Cooperative Agreement NCC 1-229, NASA - Langley Research Center 

Total Period of Performance: September 22, 1996 - September 23, 1999 
Current Period of Performance: September 23, 1998 - September 22, 1999 


Project Background: 

This report summarizes key activities conducted in the third and final year of the cooperative 
agreement NCC 1-229 entitled "Improving Conceptual Design for Launch Vehicles." This project 
has been funded by the Vehicle Analysis Branch at NASA's Langley Research Center in Hampton, 
VA. Work has been performed by the Space Systems Design Lab (SSDL) at the Georgia Institute 
of Technology, Atlanta, GA. 

Accomplishments during the first and second years of this project have been previously 
reported in annual progress reports. This report will focus on the third and final year of the three 
year activity. 


Summary of Project Accomplishments: 

From September 23, 1998 to September 22, 1999 this project focused on the assessment of 
the economics of a Bimese reusable launch vehicle when configured to support a intercontinental 
fast-package (fast cargo) delivery business. Several configurations were examined including a 
single Bimese RLV element, a single element with extra internal tanks in the payload bay, a single 
element with strap-on solid rocket boosters, and a single element with extra tanks and strap-on 
solid rocket boosters. For each configuration, the payload performance was calculated as a 
function of launch azimuth from Kennedy Space Center, Florida. A parametric business market 
was created to assess the economic viability of a Bimese vehicle operating as a fast-package 
delivery vehicle. The study results are summarized in the attached report and the accompanying 
presentation charts. 

A secondary (and less tangible) goal of the three-year project was to encourage and support 
the development of a research laboratory at Georgia Tech focused on the study of advanced space 
transportation. Toward that end, the project was also quite successful. With the support of 
government funding, the Space Systems Design Laboratory at Georgia Tech has grown from only 
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2 graduate students in 1996 to 18 graduate students and 4 undergraduate students in 1999. During 
this time, the lab has graduated 10 students with advanced degrees, developed several new tools 
and techniques applicable to the design of conceptual launch vehicles, and has published a number 
of technical papers documenting these accomplishments. 


Graduate Students Supported: 

For the current period of this project (year 3), two Georgia Tech graduate students were 
directly supported by funds derived from this cooperative agreement. These graduate students were 
awarded graduate research assistantships (GRAs) and received monthly stipends and tuition 
reimbursement. 

1) Jeffery Tooley (graduate student, GRA supported from 9/98 - 9/99) 

2) Brad St. Germain (graduate student, GRA, supported from 7/99 - 8/99) 

Jeff Tooley was supported during the entire third year of the project and is responsible for 
most of the work on the Bimese study. Brad St. Germain was supported for a summer period that 
he spent working as a summer intern in the Vehicle Analysis Branch at NASA Langley. While at 
VAB, Mr. St. Germain worked on KLIN-cycle propulsion concepts under the supervision of 
NASA's Roger Lepsch. 


Degrees Awarded: 

Two advanced degrees were awarded based partially on research work performed during this 
period of the project. 

1) Jeffrey Tooley, Master of Science in Aerospace Engineering, December 1999. 

2) Brad St. Germain, Master of Science, December 1999. 

Jeff Tooley completed his M.S. A. E. degree in the fall of 1999 and has entered the 
aerospace workforce. Brad St. Germain (who was supported by Georgia state funds during the 
academic year) earned an M. S. degree in the fall of 1999 and is continuing his education in the 
Ph.D. program at Georgia Tech. 
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Travel & Summer Activities: 

The following travel was taken in support of activities related to this project. 

1) Brad St. Germain spent 7 weeks in the summer of 1999 as a summer intern in the Vehicle 
Analysis Branch at NASA - Langley. Mr. St. Germain's travel expenses were partially 
offset by funds from the cooperative agreement. 

2) Dr. John Olds and Jeff Tooley visited NASA Langley on August 9 - 11 to discuss 
preliminary results of the Bimese project. Dr. Ted Talay provided some input and advise 
regarding the project. 

3) Dr. John Olds visited NASA Langley on November 3 and 4 in conjunction with the AIAA 
Spaceplanes and Hypersonic Technologies Conference in Norfolk, VA. Dr. Olds and 
several SSDL graduate students met with Dr. Ted Talay, Roger Lepsch, Doug Morris, 
Prasun Desai, and other NASA researchers. 


Plans for Continuing This Cooperative Agreement: 

A proposal for a continuation of this cooperative agreement has been submitted to NASA 
Langley and is currently in evaluation. If accepted, this continuation will investigate markets, 
design criteria, and vehicle concepts for economically successful Space Tourism ventures. This 
topic is of mutual interest to NASA and SSDL. 
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NOMENCLATURE 


AATe 

Architecture Assessment Tool enhanced 

CAB AM 

Cost And Business Analysis Module 

CSTS 

Commercial Space Transportation Study 

DDT&E 

design, development, testing and evaluation 

GEM 

graphite epoxy motor 

IRR 

internal rate of return 

KSC 

Kennedy Space Center 

LEO 

low Earth orbit 

LH2 

liquid hydrogen 

LOX 

liquid oxygen 

NASA 

National Aeronautics and Space Administration 

OMS 

orbital maneuvering system 

POST 

Program to Optimize Simulated Trajectories 

T/W 

thrust-to-weight 

TFU 

theoretical first unit cost 

SRB 

solid rocket booster 

SRM 

solid rocket motors 
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I. INTRODUCTION 

The Bimese architecture for low cost access to space is based on a philosophy of 
maximum commonality between reusable flight elements. This strategy is thought to lead 
to reduced development cost, fleet production cost, and operations cost. As envisioned by 
researchers at NASA Langley (Talay, et al), the nominal Bimese configuration consists of 
two identical autonomous LOX/LH2 rocket-powered wing-body RLVs mated together 
into a vertical liftoff TSTO configuration. In this nominal or “mated” configuration, the 
system is designed to deliver 60 klb of payload to low earth orbit from a launch site at 
NASA Kennedy Space Center. The “booster” RLV crossfeeds propellant into the 
“orbiter” RLV until the staging point at which the empty booster separates and performs 
an unpowered return to the launch site. The fully loaded “orbiter” continues on to orbit 
with the payload. In fact, the “booster” and “orbiter” in the mated configuration are 
identical rocket vehicles and can be interchanged from flight to flight. While this results 
in some non-optimal compromises in the design of the Bimese RLV (e.g. ordinarily, the 
booster of a TSTO system wouldn’t need a payload bay), the life cycle costs are expected 
to benefit from the high degree of commonality between the systems (e.g. same engines, 
same construction methods, same maintenance routines). In addition, the TSTO approach 
results in a smaller and less technologically risky development program relative to a 
SSTO configuration with similar technology levels. 

In addition to the nominal mated configuration, a number of complementary Bimese 
configurations have been considered to widen the payload delivery range of the 
architecture and increase payload delivery flexibility. That is, a large family of 
configurations (all using a single or multiple Bimese RLV flight elements in some way) 
could capture medium payload delivery missions (20 - 40 klb), heavy payload delivery 
missions (40 - 60 klb), very heavy payload delivery missions (>100 klb), and crew and 
cargo rotation missions to the International Space Station. Strap-on boosters, large 
expendable liquid core stages, and propellant drop tanks can be added to the basic 
Bimese flight vehicle to create these alternate configurations. 

The goal of the research outlined in this report is to consider the possibility of using 
elements of the Bimese architecture to capture part of the market in ultra-fast 
intercontinental package delivery . That is, can a single Bimese RLV (augmented or 
unaugmented) be used to fly a fast, suborbital trajectory to deliver high priority payloads 
from the United States (nominally KSC) to fictitious spaceports in European or Asian 
cities? If so, can a positive economic business case be made for this fast “point-to-point” 
market using a Bimese RLV? How much payload can be carried per flight? What type of 
annual flight rates are necessary to make an attractive economic case? 

This study was completed by researchers in the Space Systems Design Laboratory in the 
School of Aerospace Engineering at Georgia Institute of Technology between 1 998 and 
1999. 
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II. THE BIMESE LAUNCH VEHICLE 


The Bimese is a conceptual design for a fully reusable wing-body launch vehicle. It is 
the vehicle that will be used as the base element for all architectural development in the 
Bimese space transportation system, hence the name. In Figure II. 1 the Bimese is 
pictured in a three-view along with vehicle specifications. Although all of the given 
specifications are for a fully fueled vehicle with zero payload, not all configurations will 
use these parameters (e.g. in the case of ascent propellant off-load). Designed by NASA 
Langley Research Center, it was sized to deliver 60 klb of payload when launched from 
Kennedy Space Center (KSC) to a 100 nmi x 50 nmi at 28.5* inclination orbit in what is 
called the mated (Bimese) TSTO configuration (discussed in Section III.5). 





* 


93 ft 


* 



T 


25 ft 


A 


Gross weight (zero 
payload) =1,1 30,550 lb 
Dry weight = 163,100 lb 
Ascent prop. = 941,600 lb 
LOX wgtVLH2 wgt. = 6.9 


Figure II. 1 - Bimese Three- view 

With four conceptual liquid oxygen (LOX)/liquid hydrogen (LH2) engines the Bimese 
relies on the development of a new propulsion system with the parameters listed in Table 

ILL 
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Table II. 1 - Conceptual LOX/LH2 Engine Parameters 


Sea level thrust (lb) 

384,000 

Vacuum I sp (s) 

443 

Sea level T/W 

74.6 

Engine throttle 

30% 

Mixture ratio 

6.9 

Lifetime (flights) 

250 


Other design parameters that are important for the analysis of the vehicle are listed in 
Table II.2. 


Table II.2 - Design Constraints 


Acceleration limit (g) 3 

Maximum wing normal force (lb) 379,000 

Maximum dynamic pressure (lb/ft 2 ) 1,000 
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III. MISSION AND PAYLOAD OPTIONS 

The mission options of the Bimese transportation system will be analyzed in terms of 
trying to fill NASA and the commercial markets' demand for a wide variety of payload 
capabilities with minimum architectural development. In this study there are five 
varieties of the Bimese element being explored: single element, fuel-augmented, thrust- 
augmented, fuel/thrust-augmented, and the mated (Bimese) configurations. These five 
varieties are shown in Figure III. 1 - Other variants that exist, but are not being 
investigated, are the heavy lift concepts that are characterized by the addition of a large 
second or third expendable stage to any of the previous configurations. Quantifying the 
performance of each configuration will be done in terms of two parameters: due east 
point-to-point range capability for 1 klb from 28.5° latitude (KSC); and low Earth orbit 
(LEO) payload delivered to a 100 nmi x 50 nmi at 28.5 inclination orbit. For the LEO 
case there is enough orbital maneuvering system (OMS) propellant on the Bimese to 
circularize to 100 nmi x 100 nmi at 28.5 inclination orbit. 



Single element Fuel-augmented Thrust-augmented Thrust/fuel- Mated (Bimese) 

augmented 


Figure III. 1 - Bimese Space Transportation System Mission Options 


For these missions all of the trajectory analysis is done using three degree-of-freedom 
Program to Optimize Simulated Trajectories (POST). For the LEO missions POST is 
used to optimize the controls for maximum burnout weight. For point-to-point missions 
POST simulates a ballistic boost-glide trajectory, while optimizing alpha (limited to 40 
degrees) for maximum range. For fuel and thrust augmentation new component weights 
are analyzed using mass estimating relationships. More on the specifics of the vehicle 
trajectory and component weights will be introduced as each configuration is studied. 
The POST input and aerodynamic files are in Appendix A. 
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111.1. SINGLE ELEMENT BIMESE 

The single element configuration consists of the Bimese flown without any other 
components. 

111. 1.1. Single Element Bimese: LEO 

It is determined that the single element configuration cannot m ake it to LEO before 
burning out of ascent propellant. In the simulation the single Bimese must use about 
40,000 lb of non-propellant mass to make it to orbit. 

111. 1.2. Single Element Bimese: Point-to-Point 

For the single element point-to-point configuration 6,200 lb of propellant is off-loaded 
which corresponds to the OMS fuel needed to circularize and de-orbit. The single 
element point-to-point simulation shows that the Bimese can transport 1 klb of payload to 
28.5' latitude with a range of 3,900 nmi. A ground track of the trajectory is shown in 
Figure III.2. Also shown in this figure are the approximate landing locations for launches 
in all directions and the single element ranges loaded with 60 klb of payload. The figure 
depicts that the added payload weight reduces the range by about 30% and launching in a 
westerly direction (as compared to an easterly one) reduces the range by approximately 
30%. 



Figure III.2 - Ranging and Ground Track for Single Element Bimese 
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Five plots of the reference point-to-point trajectory are seen in Figure III.3a through 
III.3e. The plots show, in order, a time history of the vehicle’s altitude, range, angle of 
attack, acceleration, and velocity. Figure III.3a shows that the trajectory stops at an 
altitude of 50,000 ft; once the vehicle has reached this altitude and a Mach number of 
about 1 the trajectory is assumed to be an automated loiter and landing phase that is not 
simulated. Another feature of this plot is the skipping trajectory; by skipping across the 
upper atmosphere the vehicle can obtain maximum range. As can be seen from the range 
plot most of the ranging is done during this skipping phase. Unfortunately this skipping 
also corresponds to high aerodynamic heating. Also note that the vehicle travels 3,900 
nmi in 37 minutes, which gives an average speed of 8,000 miles per hour. This high 
average flight speed is important for applications of the point-to-point trajectories and 
will be discussed in more detail later. 

III. 1.3 Single Element Bimese Conclusions 

The single element Bimese has no LEO capability, but it does have a marginal trans- 
Atlantic point-to-point range with small payloads (under 1 klb). Trans-continental 
service between the coasts of the United States is also a possibility, setting aside the 
obvious problems associated with land overflight (noise, safety, and regulatory issues), 
but the performance of a single element Bimese falls short of enabling trans-Pacific 
service. To snnnort a more ambitious . long ranee f as t package delivery service, t he 
sinyle element Bimese would have to be sized up bv roughly 25% - 35% ; 

Given its small payload and range capability, the single element Bimese configuration 
might best be used to offer "charter style" trans- Atlantic service for extremely high value, 
but lightweight cargoes on an as-needed basis (see the economic analysis in Section IV). 
In addition, the suborbital trans-Atlantic trajectory flown by the simple Bimese would be 
excellent for testing a Bimese prototype, putting it through many of the extremes that an 
orbital vehicle would encounter. 
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III.2. FUEL-AUGMENTED BIMESE 

The fuel-augmented Bimese configuration consists of a single element with an extra LOX 
and LH2 tank in its payload bay. The payload bay is a cylinder, 50 ft long with a 15 ft 
diameter. If it is fully loaded with propellant, it could hold about 200,000 lb of 
LOX/LH2. 


III.2.1. Fuel- Augmented Bimese: LEO 

The fuel-augmented Bimese for the LEO mission for still cannot get any payl oad to orbit. 
Although the increased fuel added velocity increment capability, the fact that no 
additional thrust is added causes the overall thrust-to-weight at liftoff (T/W) to decrease 
as propellant is added. Even assuming that a liftoff T/W of 1 .05 provides adequate liftoff 
margin, the fuel-augmented Bimese must bum about 20,000 lb of non-propellant mass to 
get to orbit. 


III.2.2. Fuel-Augmented Bimese: Point-to-Point 


Simulation of point-to-point fuel-augmented trajectories show that even a payload bay 
stuffed full of fuel only increases the point-to-point range by about 400 nmi. Combine 
this with the fact that for minimum architectural developments, the same payload tanks 
will be used for the fuel/thrust-augmented vehicle (which is expected to have a capability 
of about 20,000 lb to LEO) leads to the choice of filling three-fifths of the payload bay 
with tanks. This leaves room for about 20,000 lb of payload and adds only 200 nmi in 
range to the fuel-augmented point-to-point trajectory. A line drawing of the fuel- 
augmented Bimese along with payload tank specifications is shown in Figure III.4. The 
fuel-augmented point-to-point range for 1 klb to 28.5' latitude is 4,100 nmi, no plots will 
be shown for this because they are very similar to the single element plots. 



Figure III.4 - Fuel-Augmented Bimese with Payload Tank Parameters 
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111.2.3. Fuel-Augmented Bimese Conclusions 

With the assumptions made on available fuel load, the fuel-augmented Bimese offers 
little advantage over the single element Bimese. It has no LEO capability and the point- 
to-point range is only increased by a few hundred nautical miles. In addition it has the 
operational disadvantage of having to install, fill, and purge during ascent the payload 
bay tanks. 

111.3. THRUST-AUGMENTED BIMESE 

The thrust-augmented Bimese configuration consists of a single element with Solid 
Rocket Motors (SRMs) strapped to the side of the Bimese. A new motor will need to be 
designed to fill this piece of the transportation system. Keeping the design realistic the 
new motor is modeled as a derivative of the Graphite Epoxy Motor (GEM), which is 
currently used on the Delta II 7925. The GEM is chosen because of its good performance 
and lightweight structure, which are shown in Table III. 1 . 

Table III. 1 - Reference GEM Performance and Design Parameters (Delta) 


Propellant mass (klb) 

25.8 

Gross mass (klb) 

28.6 

Sea level thrust (klb) 

99 

Sea level I sp (s) 

265 

Bum time (s) 

63.0 

Expansion ratio 

10.7 

Overall length (ft) 

42 

Core diameter (ft) 

3 


Scaling the GEM involves the use of simple scaling equations . To increase the bum time 
of the GEM more propellant is added, while the dry weight is scaled linearly with the fuel 
weight (dry weight is calculated to be -10% of the propellant weight). Linearly 
increasing die nozzle exit area and fuel flow rate while keeping a constant expansion ratio 
scales thrust. Specific impulse remains a constant for all of the scaling. 

A few changes are made to the single element trajectory for the simulations with SRMs. 
First a 10% drag rise is included while the motors are attached to capture some of the 
aerodynamic effects of the SRMs. Also because the SRMs provide added T/W the 
Bimese accelerates much faster resulting in violation of the dynamic pressure constraint 
listed in Table II.2. To alleviate this problem the main engines are throttled to a constant 
value while the SRMs are thrusting; this throttle value is optimized within the trajectory 
simulation. 

III.3.1. Thrust- Augmented Bimese: LEO 

With thrust augmentation the single Bimese can finally make it to orbit . In order to 
investigate the ability of the thrust-augmented Bimese to ferry payload to LEO, a design 


Georgia Tech SSDL 


Year 3 Final Report, Improving Conceptual Design for Launch Vehicles 





NCC 1-229 


The Bimese Concept: A Study of Mission and Economic Options 


14 


of experiments analysis is performed. Both SRM bum time and total SRM sea level 
thrust are varied and LEO payload is observed. Based on initial simulation experimental 
ranges of 75 to 125 sec. for bum time and 1,500 to 2,000 klb for total SRM sea level 
thrust are chosen. The lower limits are set by zero payload capability. The upper limits 
are set by throttle limit violation when trying to meet the dynamic pressure constraint. 
Table III.2 shows the results of the design of experiment (a negative payload indicates the 
vehicle cannot make it to orbit). A response surface with a mean square of 0.999 is 
generated and plotted in Figure II1.5. 


Table III. 2 - Experimental Design for Thrust- Augmented Bimese 


SRM bum 
time (s) 

Total SRM Sea 
Level Thrust (klb) 

Payload to 
LEO (lb) 

75 

1,500 

-7,896 

100 

1,500 

-1,727 

125 

1,500 

1,123 

75 

1,750 

-1,553 

100 

1,750 

5,742 

125 

1,750 

8,712 

75 

2,000 

3,651 

100 

2,000 

11,939 

125 

2,000 

17,077 



Payload to LEO 



Figure III. 5 - LEO Payload to Orbit for Thrust-Augmented Bimese 
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This figure illustrates that the thrust-augmented single Bimese can get anywhere between 
1 and 20 klb to LEO. It also shows that many SRM designs can fill a single payload 
requirement to LEO. That is, the solution for a given payload capability is a non-unique 
combination of total SRB sea-level thrust and bum time. To determine which SRM 
design is best (here, lowest SRB gross weight) for each payload capability, a plot of SRM 
gross weight for the 5 and 10 klb payload thrust-augmented Bimese is introduced in 
Figure III.6. 



Total SRM Sea Level Thrust klb) 

Figure III.6 - SRM Gross Weight Change for Constant Payload to LEO 

For a constant payload, the minimum SRB weight occurs with maximum allowable 
thrust. Therefore the maximum thrust value that is within the ranges of the experimental 
design will be used for thrust augmentation. That is, the total sea-level SRB thrust is set 
to 2,000,000 lb. in subsequent analysis and SRB bum time is varied to get the desired 
payload capability. 

The number of individual of SRB's for the 5 and 1 0 klb payload cases was chosen to be 4 
to yield a symmetric configuration and a reasonable SRB size (500 klb thrust at sea-level 
per SRB). Figure III.7 shows the GEM derivatives and the resultant bum times for each 
payload capacity (GEM- 10 for the 10 klb payload case and GEM-5 for the 5 klb payload 
case). The figure also compares the Bimese GEMs with the Shuttle SRB and the 
reference GEM from the Delta 7925. 
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A 


SRM name 

Shuttle SRB 

GEM-10 

GEM-5 

GEM 

Vehicle use 

Two on the Space 
Shuttle 

Four on 10 klb to 
LEO thrust- 
augmented Bimese 

Four on 5 klb to 
LEO thrust- 
augmented Bimese 

Nine on 
Delta 
7925 

Propellant mass (klb) 
Gross mass (klb) 

Sea level thrust (klb) 
Sea level I $p (s) 

Bum time (s) 
Expansion ratio 
Overall length (ft) 
Core diameter (ft) 

1,107 

1,300 

2,650 

267 

124 

7.5 

150 

12 

190 

210 

500 

274 

92.8 

10.7 

60 

6 

155 

175 

500 

274 

78.2 

10.7 

50 

6 

25.8 
28.6 

98.9 
274 
63.0 
10.7 
43 
3 


Figure III.7 - SRM Size and Performance Comparison Chart 

The Bimese with four GEM- 10s, seen in Figure III-8. will be used as the reference thrust- 
augmented Bimese. 



Figure III. 8 - Thrust-Augmented Bimese with Four GEM-1 Os 

A verification run was performed to obtain the actual payload (instead of the regressed 
payload) obtained when four GEM- 10s with the parameters in Figure HI-7. For the 
verification run, the thrust-augmented single Bimese was found to insert 9,740 lb of 
payload into LEO. 
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Altitude, throttle setting, dynamic pressure, and acceleration of the simulated ascent 
trajectory are seen in Figures III. 9a through d. In Figure III.9b notice the throttle setting 
for the main engines is set to a constant value at twenty seconds up until the solids are 
dropped. Also notice in Figure III.9d the high liftoff T/W of about 1.75; this is what 
causes the dynamic pressure violation and the need for throttling, 
a) 



Time (s) 


b) 



.£ 0.2 
<a 

5 0 


0 50 100 150 200 250 300 350 400 

Time (s) 




Figure III.9a-d - Time History Plots for Thrust-Augmented Bimese to LEO 
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III.3.2. Thrust-Augmented Bimese: Point-to-Point 

The thrust-augmented Bimese with four GEM- 10s cannot fly a ballistic point-to-point 
trajectory while meeting all of the constraints listed in Table III.2. Therefore the point- 
to-point trajectory is simulated with two GEM- 10s. For the point-to-point configuration, 
6,200 lb of propellant is off-loaded from the orbital configuration which corresponds to 
the OMS fuel needed to circularize and de-orbit. The single element point-to-point 
simulation shows that the thrust-augmented Bimese can transport 1 klb of payload to 
28.5° latitude with a range of 6,050 nmi. A ground track of the trajectory is shown in 
Figure III. 10. Also shown in this figure are the approximate landing locations for 
launches in all directions. The same graphic has range capabilities for the thrust- 
augmented Bimese with 60 klb of payload. The thrust-augmented Bimese is easily 
capable of delivering 1 klb or more to a variety of spots in Europe, but still falls short of 
trans-Pacific range. From the figure, carrying 60 klb of payload weight with the thrust- 
augmented single Bimese results in a range similar to that of the 1 klb single element 
Bimese point-to-point range (i.e. marginal trans-Atlantic range). 



Figure III. 10 - Ranging and Ground Track for Thrust-Augmented Bimese 

Four plots of the point-to-point trajectory are seen in Figure III.3a through III.3e. The 
plots show in order a time history of the vehicle's altitude, acceleration, angle of attack, 
and velocity profile. Another feature of the trajectory is the extreme altitude (80 km) to 
which the vehicle ascends. The trajectory takes about 45 minutes. There are abrupt 
peaks in the velocity profile (Figure II. 1 Id) which would indicate severe heating loads. 
Because the vehicle appears to be at the limit (or perhaps even beyond in the case of 
aeroheating) of many of its trajectory constraints, no higher thrust configurations will be 
examined for longer range point-to-point missions. 
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a) 




0 5 10 15 20 25 30 35 40 45 

Time (min) 




Figure III. 1 la-d - Time History Plots for Thrust-Augmented Bimese to PTP 
III.3.3. Thrust-Augmented Bimese Conclusions 

The thrust-augmented Bimese is a success both in terms of LEO payload and point-to- 
point range. With four GEM- 10 SRBs, the vehicle can deliver almost 10.000 lb to LEO . 
The point-to-point ability has si gnificant range and pavload with the ability to boost to 
such places as western Europe and the tip of South America from KSC in under an hour . 
Its drawback is the requirement of the use of large solids that will need to be purchased 
and integrated for every flight. 
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III.4. FUELA'HRUST-AUGMENTED BIMESE 

The fuel/thrust-augmented Bimese configuration consists of a single element Bimese 
with fuel tanks in the payload bay gnd SRMs strapped to the side. The same designs for 
the fuel-augmented Bimese payload tanks (99,900 lb. of propellant per Figure II1.4) and 
thrust-augmented Bimese SRMs are used on this vehicle for the LEO mission. As 
mentioned earlier, this configuration was not examined for point-to-point range since it 
would easily exceed the trajectory constraints on a programmed suborbital mission. Since 
this configuration can easily make low earth orbit, one might think of it (and the four 
SRB thrust-augmented Bimese in the previous section) as a global r^ngg fast-package 
delivery vehicle with a point-to-point payload equal to its orbital payload capability. 

III.4.1. Fuel/Thrust-Augmented Bimese: LEO 


Similar to that performed for the thrust-augmented Bimese, a design of experiments 
analysis is performed on the thrust/fuel-augmented Bimese to test its LEO capability. 
The same range of 75 to 1 25 sec. for bum time and 1,500 to 2,000 klb for total SRM sea 
level thrust are chosen. Table III.3 shows the results of the design of experiments matrix 
(a negative payload indicates the vehicle cannot make it to orbit). A response surface 
with a mean square of 0.985 is generated and plotted in Figure III. 12. 


Table III.3 - Experimental Design for Fuel/Thrust-Augmented Bimese 


SRM bum 
time (s) 

Total SRM Sea 
Level Thrust (klb) 

Payload to 
LEO (lb) 

75 

1,500 

-1,751 

100 

1,500 

4,507 

125 

1,500 

7,774 

75 

1,750 

3,562 

100 

1,750 

10,179 

125 

2,000 

15,566 

75 

2,000 

9,960 

100 

2,000 

19,672 

125 

2,000 

22,323 


Georgia Tech SSDL 


Year 3 Final Report, Improving Conceptual Design for Launch Vehicles 




NCCl-229 


The Bimese Concept : A Study of Mission and Economic Options 


21 



Payload to LEO 


•■•■O' 

Oklb 

— 

— 5 klb 

X • 

10 kb 

1 — 

- 15 kb 

A 

20 kb 


Total SRM Sea Level Thrust (klb) 

Figure III. 12 - LEO Payload to Orbit for Fuel/Thrust- Augmented Bimese 


This figure illustrates that the fuel/thrust-augmented Bimese can get anywhere between 1 
and 25 klb to LEO. For similar sized SRMs, it can deliver about five thousand more 
pounds of payload over the thrust-augmented Bimese without additional tanks in the 
payload bay. For commonality with the thrust-augmented Bimese, the fuel/thrust- 
augmented design will use four previously sized GEM- 10s (500 klb sea level thrust each 
with a 92.8 sec. bum time, from Figure III.7). Verifying the response surface prediction 
with an actual trajectory simulation with these SRMs gives a payload of 15,900 lb. A 
picture of this Bimese configuration is shown in Figure III. 13. 



Figure III. 13 - Fuel/Thrust-Augmented Bimese 
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Although there are tanks in the payload bay, there is still a 20 ft long 1 5 ft diameter space 
in the bay, which is enough room for the reference payload. Time history plots of the 
ascent trajectory are not shown because they are very similar to the thrust-augmented 
Bimese plots. 

III.4.2. Fuel/Thrust-Augmented Bimese Conclusions 

Fuel augmentation increases the payload of thrust augmentation by about 5,000 lb, which 
is attractive considering the only added costs are fuel tanks, operational complexity, and 
fitting the Bimese with proper payload feed lines. If a global range fast package delivery 
is desired, this configuration could delivery 15,900 lb. to orbit, and then deliver it to 
nearly any location on the globe. However, the four SRB's and the internal fuel tanks 
make this configuration more expensive to operate than any of the previous concepts 
considered. 


III.5. MATED BIMESE 

For completeness, the nominal mated Bimese configuration was examined. The mated 
Bimese configuration consists of two single elements attached to each other. Pictured in 
Figure III. 14 the mated Bimese has the following characteristics: propellant cross-feed, 
un-powered fly back, and commonality between booster and orbiter. 



Orbiter 

Gross weight (zero payload): 1,083,700 lb 
Ascent propellant weight: 909,700 lb 


Figure III. 14 — Mated Bimese 


III.5.1. Mated Bimese: LEO 

The Bimese was designed so that the mated configuration could lift 60 klb to LEO. Upon 
launch fuel is cross-fed from the booster Bimese to the orbiter Bimese, while all eight 
engines are ignited. The booster stages at Mach 3.2 (propellant must be off-loaded to do 
this) and then performs a hypersonic turn and glides back to the launch site. From here 
the orbiter ascends to orbit, with the benefit of a Mach 3.2 boost. Plots of the ascent 
trajectory are seen in Figure III.15a-c. 
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Figure III. 1 5a-c - Time History Plots for Fuel/Thrust- Augmented to LEO 


III.5.2. Mated Bimese Conclusions 

The large payload capability of the mated Bimese makes it a candidate for NASA 
missions to the space station or heavy lift missions such as the deployment of a large 
space telescope. The commercially market for this size payload is unproven, but it is 
possible it could be used for multiple manifesting or geo-stationary payload delivery. 
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III. 6. SUMMARY OF MISSION OPTIONS 

A fleet of vehicles has now been established that can deliver a wide range of payloads to 
orbit from Kennedy Space Center and (at least from a performance perspective) can 
deliver cargoes to terrestrial spaceports in support of a ultra-fast package delivery service. 
Except for the fuel-augmented Bimese each of the vehicle seems to have a niche where it 
could be used to enter the commercial or government launch market. The vehicles are 
shown in Figure III. 16 along with the capabilities. 



Mission 

Single 

element 

Fuel- 

Thrust- 

augmented 

Thrust/fuel- 
augmented 
(4 GEM- 10s) 

Mated 

augmented 

# GEM-1 Os 

(Bimese) 



2 

4 


LEO Capability 
(lb)" 

0 

0 



15,900 

60,000 

Point-to-point 
capability (nmi)* 

3,900 

4,100 


t 

t 

+ 


** LEO payload delivered to a lOOnmix 100 nmi at 28.5' inclination 
* Point-to-point range capability for 1 klb to 28.5° latitude 

f These vehicles have essentially global range with a point-to-point cargo equal to their LEO payload 


Figure III. 16 - Bimese Payload and Mission Options Summary 

For global range fast-package delivery, any of the configurations with orbital capability 
would meet the requirements. The thrust-augmented configuration with 2 GEMs 
provided the longest range and highest payload capacity in a suborbital configuration 
(although still short of trans-Pacific). The simplest and cheapest configuration, the single 
element Bimese, is too small to provide more than lightweight cargo delivery (< lklb) to 
extreme western Europe from KSC. 
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IV. ECONOMIC ANALYSIS 

One of the goals of this study is to investigate the economic prospects for a profitable 
ultra-fast package delivery business using components of the Bimese architecture. All of 
the following economic analyses will assume the creation of a startup fast-package 
delivery company named Bimese PTP, Inc. Four economic scenarios will be considered 
with assumptions ranging from conservative to very optimistic. The success of each 
scenario will be quantified in terms of an Internal Rate of Return (IRR) that is calculated 
from an internal cash flow analysis (revenues - expenses and taxes). In the aerospace 
sector, a business venture is ordinarily considered to be economically attractive if the IRR 
is greater than about 25%. Note that for most of the point-to-point scenarios considered, 
any operating synergies between Bimese PTP, Inc. and companies using Bimese 
elements in the LEO crew/cargo delivery market are ignored. That is, the Bimese PTP, 
Inc. company is considered a separate company and does not share personnel or flight 
vehicles with LEO missions. However, three of the four scenarios examined assume 
some reusable element production and design cost synergy with an unnamed LEO 
company. 

IV. 1. THE POINT-TO-POINT FAST PACKAGE MARKET 

For the emerging, intercontinental ultra-fast package delivery market, the advantage for 
using rocket over aircraft delivery lies in the former's speed of flight. The time for 
checkout, truck delivery, and port delay for the two modes are the same. Therefore to see 
any real advantage over aircraft, the Bimese must ship goods fast and far to places where 
there is a large difference between aircraft trip times and boost-glide trip times. Table 
IV.3 shows the advantage rocket delivery has over airplane delivery in terms of trip time 
for selected trans-Atlantic and KSC - South American city pairs. Representative prices 
from the current aircraft-based fast package market are also given (about $ 10/lb for a 2 
lb. package). 


Table IV. 1 - Time of Flight of Aircraft versus Bimese 


From 

To 

Approximate flight 
time (hr) 

FedEx next 
day air cost 

($/lb) 

Aircraft 

Bimese 

KSC 

Madrid 

8 

~0.7 

10 

KSC 

Los Angeles 

3 

-0.5 

7 

KSC 

London 

8 

-0.7 

10 

KSC 

Rio De Janeiro 

10 

-0.7 

12 

KSC 

Paris 

9 

-0.7 

13 


From Table IV. 1, a rocket-based system can reduce the flight time for package delivery 
by over 90% in some cases (longer ranges show better reductions). While this does not 
necessarily translate to significant door-to-door time reductions, it is clearly an advantage 
for rocket systems. Of course the Bimese has many disadvantag es compared to aircraft 
systems including: no existing spaceports for vertical launch, high turn around times. 
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high startup costs, land overflight issues at lower altitudes, uncertain reliability, and the 
use of an unproven vehicle concept. Economically these disadvantages will be hard to 
overcome for Bimese PTP, Inc. 

IV. 1 . 1 . The Size of the Ultra-Fast Package Delivery Market 

The basic assumption in this analysis is that shippers will pay a premium for speed. Just 
as companies such as Federal Express and Emory World Wide became economically 
successful by offering next-day service within the U.S. and Europe, many believe that 
there is a market for next-day (or even same-day) service between long range city pairs. 
The size of this market will depend on the price offered to shippers. That is, the ultra-fast 
package is assumed to be price elastic. Unfortunately, no universally accepted data exists 
for the ultra-fast package delivery market 

For the purposes of this study, the fast-package delivery data from NASA’s Commercial 
Space Transportation Study (CSTS) will be used for order of magnitude analysis. In 
general, the CSTS data can be summarized by Figure IV. 1. Note that this figure uses a 
log-log scale. Here one can see the upper and lower bounds of the CSTS predicted 
markets. 



If we assume that the average configuration used by Bimese PTP, Inc. will ship an 
average of 5,000 lb. per flight and set the lower limit on average package price at $ 10/lb. 
(current aircraft system prices) and the upper limit one order of magnitude higher 
($1 00/lb), then a region of the CSTS data can be bounded within the CSTS data. Note 
that the CSTS data for this market consists of an upper limit and a lower limit. The 
bounded region corresponds to a revenue of between $0.5 M and $5 M per flight and a 
total number of annual flights of the system (all routes) of between 200 and 8000 flights 
per year or about 0.5% to 50% of the current aircraft-based fast package delivery market. 
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Higher package prices (higher premiums over today's prices) result in very low traffic 
rates and seem unlikely to justify the required investment to create Bimese PTP, Inc. 
Offering lower package prices might be considered a less risky way to enter new market. 
That is, it is more reasonable to expect current shippers to pay a premium of 25% - 100% 
for faster speed (say, up to $20/lb.) than it is to expect shippers to pay prices a factor of 
10 higher than current prices. The former is more of an incremental market expansion, 
while the later approach would depend almost exclusively on enticing new, previously 
unserviced, shippers to enter the market (e.g. shipment of biological specimens, volatile 
pharmaceuticals, or other priority freight). However, it is unlikely that Bimese PTP, Inc. 
could offer prices close to or lower than that offered by today's aircraft given the extra 
expenses of operating a rocket relative to an aircraft. 

Based on these qualitative arguments, the ideal "sweet spot" of the CSTS fast-package 
market for a rocket system would seem to lie in the middle to lower right of the box 
indicated in Figure IV. 1 — roughly 2000 flights per year at a revenue of $1M per flight. 
Whether the system is actually profitable in that range will be determined in a subsequent 
section of this report. 

IV. 1.2. The "Ideal" Service Scenario 

Before directly considering the economics of Bimese PTP, Inc., an "ideal" service 
scenario was envisioned based on the point-to-point pavload performance of the Bimese 
configurations previously examined. 

In order to be successful, Bimese PTP, Inc. will likely need to offer regular scheduled 
service between selected city pairs or spaceports between economically developed parts 
of the world. Daily weekday service each way between two spaceports would yield 
about 500 flights per year. Four routes would be needed to serve the 2000 flight per 
market estimated for the CSTS sweet spot above. 

Based on performance, Bimese PTP, Inc. would prefer to offer a range of flight 
configurations. The longest range and heavy cargo route(s) would be serviced by the 
thrust-augmented single element Bimese with four GEM- 10 SRB's. This vehicle would 
in fact provide global range with up to almost 1 0,000 lb. of payload. The intermediate 
route(s) would be served by the thrust-augmented single element Bimese with only two 
GEM- 10 SRBs. This configuration provides a sufficient range to reach mid-Europe from 
KSC with payloads up to 40 klb. - 60 klb. (and could marginally reach the Pacific-rim 
from the west coast of the U.S. with small payloads up to 1 klb.). Relative to the four 
GEM- 10 option, the two GEM- 10 configuration is expected to be cheaper to operate and 
is therefore the preferred workhorse (for a performance point of view). The mated 
Bimese configuration requires two flight elements per mission and is thought to be too 
"overkill" for servicing the point-to-point market. 

Given its limited range, the single element Bimese by itself is undersized for use on a 
regularly scheduled service route. Therefore, it is more likely to be used in limited role. 
A standby "charter" mode for very high priority, but lightweight cargoes is envisioned. 
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From the east coast of the U.S., the single element Bimese might be used to fly cargoes 
up to 1 klb to extreme western Europe or to South America. No more than 50 or so of 
these charter flights are expected per year (once per week). 

While this "ideal" service scenario would provide the requisite level of service to enable 
the ultra-fast package delivery market, it remains to be seen whether it results in a 
economically attractive scenario for Bimese PTP, Inc. 

IV.2. BIMESE PTP, INC. COST AND OPERATIONS 

To estimate IRR for a business scenario, the cost and revenue cash flows must be 
predicted as a function of flight rate. Bimese PTP, Inc. will incur both non-recurring 
(start-up) costs and recurring costs (variable costs). Non-recurring costs typically consist 
of facilities construction costs, flight vehicle DDT&E costs (design, development, testing, 
and evaluation), and reusable flight vehicle procurement costs (related to TFU, 
theoretical first unit cost, and required fleet size). Recurring costs include procurement 
of expendable hardware, ground labor costs, maintenance hardware costs, propellant 
costs, and insurance costs. Financing costs (interest on debt) are also accounted for in the 
cash flow analysis. 

IV.2.1. Economic Assumptions 

Bimese PTP, Inc. is assumed to begin to develop the Bimese launch vehicle with 
government help in 2007; the entire Bimese program ends in 2037. The transportation 
system created is one full generation beyond the Shuttle Transportation System in terms 
of the technologies utilized (Gen2 RLV). Production and operations learning curves of 
85% have been assumed. All of the cost and business analysis is forecasted using Cost 
And Business Analysis Module (CABAM) developed at Georgia Tech s SSDL and all of 
the operations modeling is performed using Architecture Assessment Tool-enhanced 
(AATe) vl.O developed by NASA KSC. Monetary units will be given in terms of 
constant year 1998 United States dollars. 

The TFU and DDT&E for the Bimese launch vehicle were estimated using component- 
level weight-based NASA Air Force COst Model (NAFCOM) equations stored in 
CABAM. The results obtained for the Bimese airframe and engine are listed in Table 
IV.2. The values shown are prior to taking a deduction for government contribution and 
include a 20% cost margin for conservatism. The TFU listed for the Bimese liquid 
engine is for a single engine, not for a complete flight set of 4. The engines are assumed 
to have a lifetime of 250 flights before they need to be replaced and each airframe has a 
life of 1,000 flights before it needs to be replaced. For all business cases it will be 
assumed that, at the very least, the government pays for 20% of the airframe DDT&E and 
100% of the engine DDT&E. For the mission option with payload tanks, the tanks are 
assumed to be reusable, and the DDT&E and TFU for these tanks are also listed in Table 
IV.2. 
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Table IV.2 - Bimese Launch Vehicle Costs 


Component 

DDT&E (M$) 

TFU (M$) 

Airframe 
Main engine 
Payload tanks 

$ 6,950 
$450 
$95 

$ 1,431 
$ 109 
$20 


It is assumed that a separate company makes the GEM- 10s and charges Bimese PTP Inc. 
a fixed price for the SRB motors. The current GEM on the Delta has a fixed price of 
about $1 .2M. Taking into account an expected price benefit from volume production and 
also the fact that the GEM- 10 is scaled up by a factor of about five in terms of mass and 
thrust, the GEM-10's will be assumed to cost $2M each. As will be shown, this estimate 
has a significant effect on the economics of the thrust-augmented single element Bimese 
options. 

Using AATe (Architecture Assessment Tool from NASA KSC), the turnaround times 
between flights of the various Bimese configurations are estimated and listed in Table 
IV.3. AATe is grounded in Space Shuttle data and estimates the ground processing time 
of advanced space vehicle concepts based on system-by-system order-of-magnitude 
comparisons between the current concept and the Space Shuttle. The single element 
Bimese is the simplest of the configurations examined and is estimated to have an 
average turnaround time of just under two- weeks (leading to an annual flight of about 3 1 
flights per airframe per year). By comparison, a shuttle orbiter requires about 12-16 
weeks of processing between flights. For the Bimese configurations, there is a slight 
increase in turnaround time, as the integration becomes more complex. It is assumed for 
the mated Bimese that the preparations before launch and after landing can be done on 
the two elements simultaneously. Note that while a two week turnaround time is good 
compared to the current space shuttle, it is no where near the turnaround time of a typical 
long range aircraft (which can be measured in hours). This variable has a significant 
effect on required fleet size for Bimese PTP, Inc. To support the "ideal" scenario of 2000 
flights per year , the required fleet size would be about 70 airframes! 

Table IV.3. - Turnaround Times for Bimese Configurations 


Architectural Configuration 

Turn around 
time (days) 

Yearly 
flight rate 

Single element 

12 

31 

Thrust-augmented 

13 

28 

Fuel-augmented 

13 

28 

Thrust/fuel-augmented 

14 

26 

Mated 

16 

23 


Georgia Tech SSDL 


Year 3 Final Report, Improving Conceptual Design for Launch Vehicles 



NCC 1*229 


The Bimese Concept: A Study of Mission and Economic Options 


30 


Non-recurring construction of facilities cost are assumed to be paid for by a local or 
national government investing in a spaceport for possible economic benefits to a region. 
That is, the spaceport is assumed to be a comparable fixed infrastructure to a regional 
airport and Bimese PTP, Inc. is analogous to an airline operating from a government- 
owned spaceport. Operating from this spaceport, Bimese PTP Inc. is assumed to account 
for 15% of the total flights and pays a nominal "fee" of $50k - $100k per flight to use the 
facilities, launch pads, and runways of the spaceport. 

IV.2.2. Recurring Cost Results and Comparisons 

Expendable GEM- 10 cost, ground personnel labor costs, maintenance hardware cost, 
propellant cost, stage integration, and insurance costs sum up to total recurring costs per 
flight. Economies of scale are assumed to benefit recurring costs slightly. Figure IV.2 
shows the estimated magnitude of these variable costs for each of the five Bimese 
configurations examined. 



Yearly Flight Rate 

Figure IV.2 - Average Recurring Costs for Multiple Configurations 

It is quite clear the major effect on the recurring cost is the purchase cost of the GEM- 10 
SRBs at $2M each. Given that the total revenue goal from section IV. 1.1 is only on the 
order of $1M per flight, any configuration with the recurring cost over $1M is 
economically infeasible (before amortization of non-recurring costs are even considered). 
The fully reusable single element and mated Bimese configurations pass this first-order 
cost filter, but all configurations with expendable SRB's do not. This requires a major 
rework of the "ideal" service scenario of section IV. 1.2. 
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IV .2.3 The "Realistic" Service Scenario 

The discussions above create a dilemma for Bimese PTP, Inc. Due to range and payload 
requirements, a full global range fast-package service with regularly scheduled routes 
requires a thrust-augmented Bimese configuration. The expected ultra-fast package 
market only reaches appreciable size when package prices are around $20/lb. so revenue 
is limited to under $1 M/flight. However, because the thrust-augmented Bimese 
configurations use expendable SRB's, their cost quickly exceeds the expected revenues. 

Assuming that the single element Bimese design cannot be sized up to carry more 
payload or achieve a longer range and further assuming that the use of the mated 
configuration for fast package delivery will not be competitive since the fleet size will 
essentially be doubled at each site (doubling the non-recurring costs that must be 
amortized each flight), the following scenario is established. Bimese PTP, Inc. will offer 
only "charter-style" service on single element Bimese vehicles from KSC to points within 
its 3,900 nmi range. Payload will be limited to the 1 klb. capacity of this configuration. 
Note that the fuel-augmented version of the Bimese could be used where necessary to 
extend the range slightly. This option is not modeled explicitly, but the results would not 
be expected to change the "realistic" service economic analysis results significantly. 

Bimese PTP Inc.'s realistic service scenario is to offer chartered flights of the single 
element Bimese across the Atlantic Ocean . To begin the Bimese PTP, Inc. will operate 
from four launch sites and as the years go buy ramp up building sites inland on the two 
continents. Figure IV.3 shows a representative scenario for the Bimese fast package 
delivery system. 



Figure IV.3 - Possible Initial Routes for Bimese PTP, Inc. Fast Package Delivery 
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IV.3. POINT-TO-POINT BUSINESS ANALYSIS RESULTS 

When degrading the level of service from the regularly scheduled "ideal" case to the 
trans-Atlantic chartered flight scenario forced by performance and cost considerations, 
the CSTS price elasticity curves used in section IV. 1.1 can no longer be relied upon to 
estimate the number of flights for a given package price. That is, the number of chartered 
trans-Atlantic fast-package flights for the single element Bimese is expected to be 
significantly less than the global fast-package market in the CSTS model. However, the 
CSTS market demand can serve as an upper bound sanity check for further analysis. 

Given the uncertainty in the market size vs. package price, the following economic 
analyses treated the market elasticity parametrically. The package delivery price and the 
flight rate are varied to determine what they need to be to result in a positive economic 
result for the realistic Bimese PTP, Inc. operating scenario. The requirements are then 
compared to the traffic levels predicted for global service from the CSTS model. 

For estimating the cash flow and internal rate of return of Bimese PTP, Inc, four different 
scenario assumptions were considered ranging from conservative to very optimistic. 
Recall that the vehicle turnaround time has a major impact on fleet size. Non-recurring 
costs (DDT&E and reusable hardware production, related to TFU) must be recouped over 
all flights in the model and can be a significant contributor to total cost per flight. For 
example, a $25 B non-recurring investment for a system that flies 2000 flights per year 
for 30 years would contribute over $400,000 to die per flight cost of a fast-package 
delivery flight (not even counting the interest cost associated with the debt!). With these 
contributors in mind, the following four scenarios are examined: 

1 . standard non-recurring cost, standard turnaround time (baseline) 

2. reduced turnaround time case 

3. reduced non-recurring cost and reduced turnaround time 

4. zero non-recurring cost 

IV.3.1. Baseline Case Economic Results 

The economic results for the baseline case are shown parametrically in Figure IV.4. 
Here, the baseline assumptions that the government contributes 20% of airframe DDT&E 
and 100% of engine DDT&E are used. The AATe prediction of turnaround time of 12 
days for each single element Bimese is also used. Finally, the variable costs from Figure 
1V.2 are used for the recurring cost estimates. 

Recall that an acceptable IRR is around 25%. Figure IV.4 shows that for the baseline 
case, Bimese PTP Inc. requires about $7M per flight of revenue to achieve this goal when 
having a flight rate of about 5,000 per year (the highest flight rate considered in this 
analysis). At this high price ($7000/lb.), the CSTS market for the entire global fast 
package market predicts a demand of no more than 20 - 30 flights per year! Clearly, the 
baseline assumptions will not result in a positive business case. 
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One problem with the baseline assumptions is the slow turnaround time. To support 
5,000 flights per year over the entire scenario, extra Bimese RLV's must be purchased at 
each site to fill in for vehicles that are being processed for return to flight. For the 
baseline scenario, 171 single element Bimese RLVs would be required to be purchased. 
This leads to a fleet acquisition cost of $154 B! Even if the economic return on 
investment was attract, this startup cost is well out of the realm of feasibility for a private 
startup business venture. 

IV.3.2. Reduced Turnaround Time Scenario 

As a sensitivity, the second case analyzed assumes that the turnaround time of a single 
element Bimese RLV can be reduced significantly to only one day . This is comparable to 
assuming "aircraft-like" operations for a reusable launch vehicle. In this case, required 
fleet size is no longer governed by the length of time it takes to ready a flight vehicle for 
another flight (the processing queue), but rather by the life of the airframe. The life of 
the airframe was earlier assumed to be 1,000 flights before replacement (still low by 
aircraft standards, but an order of magnitude longer than the Space Shuttle). 

The IRR sensitivity for this scenario is shown in Figure IV.5. Again, the desire to obtain 
the lowest package price drives the flight rate up significantly. At 5,000 flights/year (the 
highest flight rate considered), the price associated with 25% IRR is $6M per flight. 
Again, 5,000 flights/year is well outside for the number of flights predicted by the global 
market for this price. 

At 5,000 flights/year, the assumption of a 1 day turnaround time reduces the fleet size to 
144 airframes for the nearly 30 year mission model (with a ramp up to full operating 
capacity between 201 1 and 2012). Compared to the baseline case, the fleet acquisition 
cost is lowered by $51 B, but the total cost is still well out the range of what can be 
expected from a private venture. 
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Figure IV.5 - IRR for Bimese Inc. with Reduced Turnaround Time 


In Figure IV.5, the benefits of fast turnaround time are more evident at low flight rates 
where turnaround time determines fleet size rather than airframe life. Lower flight rats 
will reduce overall program costs, and if they can be made economically attractive, are a 
preferred solution for Bimese PTP, Inc. However, at low flight rates, the DDT&E and 
other non-recurring costs become more of a contributor to the economic success. That is, 
DDT&E is a more significant component of life cycle cost at low flight rates. 

IV.3.3. Reduced Non-Recurring Cost and Turnaround Time Scenario 

As a third scenario of assumptions assumptions, the earlier predicted values of TFU and 
DDT&E are reduced. Bimese PTP Inc. now has 100% of airframe DDT&E paid by the 
government (or other entity) and TFU is assumed to have been reduced by 20%. These 
changes might be thought of as the result of a synergy between the single element Bimese 
used for fast-package delivery and a government configura tion used for an orbital 
mission (the mated Bimese!. For example, if the government pays for the development 
of the Bimese RLV to service the International Space Station, then Bimese PTP, Inc. 
might make use of the design work for free. Further, the government could pay a third- 
party company to create a production line for manufacturing Bimese RLVs for the 
government's mated configuration. Bimese PTP, Inc. could then purchase their single 
element Bimese vehicles from the third party at a discount. 

The results for this scenario are show in Figure IV.6. From this chart it is seen that the 
economics are becoming more reasonable, but the prices still do not match the demand. 
Here, charging a price per flight of $6M with only 250 fleet flights per year returns an 
IRR of 25%. The fleet size for this scenario is only 8 airframes and the fleet procurement 
cost falls to $15 M (DDT&E is zero to Bimese PTP, Inc.). Unfortunately, the predicted 
market size for the global market at a price of $6M ($6000/lb) is nearly an order of 
magnitude lower than 250 trans-Atlantic flights per year required to make this scenario 
successful. 
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Figure IV.6 - IRR for Bimese Inc. with Reduced DDT&E and TFU 


IV.3.4. Zero Non-Recurring Cost 

One final case is looked at with the extremely optimistic assumption of zero non- 
recurring cost. The government (or some other entity) was assumed to bear the entire 
DDT&E cost for the single element Bimese and pay the production cost to have them 
produced for Bimese PTP, Inc. The one day turnaround time assumption was also 
retained. In this scenario, Bimese PTP, Inc. pays only the direct recurring cost of 
operating the single element Bimese RLV on the trans- Atlantic charter flights. 

A plot of the IRR sensitivity for this assumption is seen in Figure IV. 7. The return on 
investment curve is relatively flat for this case (more closely related to the recurring cost 
curves) and reveals that a 25% return on investment can be obtained by charging about 
$800,000 per trans-Atlantic charter flight and flying 500 times per year or charging about 
$1M per flight and flying about 100 flights per year. The IRR = 25% curve actually 
crosses the CSTS curve for global demand for fast package services, but recall that the 
demand for charter-style trans-Atlantic flights is likely a tenth or less of the global CSTS 
market. In addition, the low payload capacity of the single element Bimese (1 klb) 
results in a package price of approximately $ 1000/lb. — a factor of 100 premium above 
current prices. Note that our previous "sweet spot" analysis of the CSTS market had 
assumed that the average Bimese configuration would carry 5,000 lb cargo globally, 
therefore $1M per flight would have been only $200/lb. The need to avoid the expensive 
SRB configurations forced the selection of the less capable single element Bimese for the 
point-to-point missions. Unfortunately, the recurring cost alone are still to expensive to 
make even the trans-Atlantic fast package market economically attractive for Bimese 
PTP. Inc. 
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IV. 3. 5. Economic Analysis Summary 

A more detailed breakdown of the economics of the four economic scenarios considered 
is given below in Table IV.4. Even with government contributions of $7.4 B and zero 
fleet acquisition costs, the price and annual flight rates that Bimese PTP, Inc. requires to 
make a 25% IRR clearly exceeds the number of flights estimated to be available at that 
price. 


Table IV.4 - Economic Indicators for Four Cases 



Baseline 

Reduced 
turn around 
time 

Reduced DDT&E, 
TFU, and turn 
around time 

Recurring 
Cost Only 

Required Flights/year 

5,000 

5,000 

250 

100 

Price/flight (M$) 

$7 

$6 

$6 

$1.0 

Approximate trans- 
Atlantic flights avail. 

<10 

<15 

<15 

<25 

IRR 

- 25% 

- 25% 

~ 25% 

~ 25% 

Fleet size 

171 

144 

8 

4 

Fleet acquisition (B$) 

$154 

$103 

$ 15 

$0 

Turn around time (days) 

12 

1 

1 

1 

DDT&E Government 
Contribution (M$) 

$1,850 

$1,850 

$7,400 

$7,400 

TFU (M$) 

$1,900 

$1,900 

$1,500 

$0 
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In the most optimistic case, the recurring costs alone prove to be too great for the 
predicted trans- Atlantic market. Here, a package price of $ 100/lb (a factor of 10 
premium over today's prices) would generate a revenue of only $100,000 per trip. For 
the single-element Bimese, direct recurring costs alone contribute over $400,000 of cost 
per trip. Paying the launch site fee of $100,000 and adding enough profit to pay taxes, 
interest, and still return a 25% IRR raises the required price to over $900,000 per flight. 
It does not appear that the Bimese will lead to a successful economic business case given 
the assumptions used in the study. 
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V. CONCLUSIONS 


This study performed an assessment of the Bimese launch vehicle architecture for use as 
a vehicle for intercontinental, ultra-fast package delivery. Specifically, the point-to-point 
payload performance of several configurations of the rocket-powered Bimese RLV 
element were assessed. Then an economic assessment was made to determine the 
revenues, costs, and internal rate of return (IRR) for a fictitious Bimese PTP, Inc. 
company created to provide the ultra-fast package delivery service. Unfortunately, the 
study was unable to identify any attractive business scenario for a Bimese-derived fast- 
package service. 

The performance analysis examined the point-to-point range capability of a single 
element Bimese RLV, a fuel-augmented Bimese RLV, a thrust-augmented Bimese 
RLV(augmented with 2 or 4 expendable GEM SRBs), a thrust and fuel-augmented 
Bimese RLV, and the nominal "mated" Bimese configuration. Of these configurations, 
the single element Bimese was found to have about a 3,900 nmi. range with a 1 klb 
payload — marginally sufficient to reach western Europe from KSC, but insufficient to 
reach the Pacific rim from the continental U.S. The configuration augmented with two 
GEM- 10 SRB's can carry significant payload on trans- Atlantic routes and might even be 
capable of trans-Pacific service from the west cost of the U.S. This two SRB 
configuration was preferred from a suborbital performance point-of-view. The larger 
configurations were all capable of making orbit, and therefore could be used to support 
long/global range or heavy cargo routes. 

A subsequent economic analysis of the expected ultra-fast package delivery market 
identified a "sweet spot" at a revenue per flight of around $1M per flight for a vehicle 
capable of carrying an average of 5,000 lb. per flight in a 2,000 flight per year global 
market. This price point represents a premium of about 100% above today's (aircraft 
system) fast package prices and would result in a service that could cut 90% of the flight 
time required for long routes. Unfortunately, the preferred mix of Bimese configurations 
relied heavily on GEM SRB thrust-augmentation to achieve the required payloads and 
ranges. At an estimated price of $2M per booster, the cost of any SRB augmented 
configurations was well above the target price of $ 1 M per flight. A case was made for 
using the single element Bimese RLV in a limited "charter-style" trans-Atlantic service, 
but even with the assumption that Bimese PTP, Inc. paid only the recurring costs of 
operating the Bimese on this route, no economically attractive solution was found. Given 
the more limited range and payload of the single element Bimese, available revenue 
would be reduced to only about $100,000 per flight and even the fully reusable single 
element Bimese proves more expensive to operate than the available revenues. Therefore 
it was concluded that no attractive economic scenario exists for the assumptions of this 
study. 

Several key parameters in the current study are responsible for this negative conclusion. 
The following observations are offered as potential areas of improvement: 
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1 . Reduce Direct Recurring Costs of Rocket Vehicles . The direct recurring cost (labor, 
propellant, hardware spares, insurance, and site fee) of a flight vehicle form the 
ultimate "price bottom" for any fast package system. Rocket systems are expected to 
be more expensive than aircraft to operate on these routes due to higher loads, higher 
complexity, smaller margins, etc. While the direct cost of a single element Bimese 
RLV at approximately $400,000/flight is already very optimistic by rocket standards, 
if rockets are to be competitive in future ultra-fast package delivery missions, 
recurring costs must be reduced bv an order of magnitude over the already aggressive 
assumptions used in this analysis . 

2. Increase Size and Range of Single Element Bimese RLV . It was shown that 
expendable SRBs could not be used for thrust-augmentation and still meet the 
stringent ultra-fast package economic goals. A fully reusable element was the only 
attractive solution . Unfortunately, the current single element Bimese (sized for the 
mated TSTO mission) is too small in both range and pavload capacity for servicing a 
global ultra-fast package market . A larger, almost single-stage capable, vehicle with 
5,000 lb. of global range payload capacity would be more optimum for this mission. 

3. Reduce or Share Initial Non-Recurring Costs . The economic scenarios examined 
highlight the critical importance of low non-recurring cost to the economic success of 
an ultra-fast package delivery market. A company such as Bimese PTP, Inc. has little 
chance of economic success if its startup costs (including vehicle DDT&E, facilities 
costs, and fleet procurement) are excessive. For a new rocket-based system, this can 
easily be 10's or even 100's of billions of dollars for ambitious global flight rates. If 
the government can provide a mechanism for reducing this non-recurring cost 
(perhaps by a synergistic development of a version of the point-to-point flight vehicle 
for civil or military orbital missions ! then the chance of economic success is greatly 
increased. Toward this end, the highly flexible Bimese architecture is an excellent 
concept for a synergistic development program. State or federal development of 
common spaceport infrastructures (analogous to regional airports) is also important. 

4. Streamline Rocket Vehicle Ground Operations . The economic scenarios examined 
highlighted the need to have very low processing times between flights of a point-to- 
point ultra-fast package vehicle. Even using optimistic assumptions, current models 
predict turnaround times of around two weeks for rocket-style point-to-point vehicles. 
However, to keep fleet sizes low and make maximum use of reusable assets to 
support a regular flight schedule, "aircraft-like" turnaround times on the order of a 
single dav are required . 
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APPENDIX A. POST DECKS 
A.l. BIMESE AERODYNAMICS 
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-10.0 1532, 
0,00*57, 

10.0 1147, 

20.0 26*9, 

30.0 5*73. 

40.1 0*43. 

50.1 7335. 

60.2 3629. 

70.2 9518, 
6. 

-20.02935, 
-10,0 1337. 
0,0 076*, 

10.0 097. 

20.0 229*. 

30.0 5245, 

40.0 9993, 


50.1 6241. 

60.2 2*95, 

70,2*552, 

10 . 

-20,0 273, 

-10,0 1233, 

0,0 0724, 

10.0 0*76. 

20.0 2109, 

30,04973, 

40.0 9643, 

50.1 5798, 

60.224*, 

70.2 *078, 

12 , 

-20.0.2717. 

-10,0 1222, 

0,0.072, 

10,00*67, 

20.0 2099, 

30,04964. 

40,09625, 

50.1 5754, 

60.2 2477, 

70.2 7915, 

15, 

-20,0 272. 

-10,0 1208, 

0,0 0704, 

10.0 0*43, 

20 , 02101 , 

30,049*6. 

40,09639. 

50.1 5742. 

60.2 2443. 

70.2 7652. 

20 . 

-20,0 269*. 

- 10,0 11 ?*. 

0,0 071, 

10.0 0827. 

20.0 2075, 

30.0 4976, 

40.0 96*8, 

50.1 5782, 

60.2 2446. 

70.2 7169. 

25, 

-20,0 266, 

-10,0 12. 

0.0 0739, 

10,00*51. 

20.0 205. 

30,04978, 

40.0 9646. 

50.1 5772. 

60.2 2428, 

70.2 6603, 

S 

ISub t*ble-6hcmat ,2,6ha)ph* ,6hm*ch ,10,15,8*1, 

0, 

-10.0 0327, 

-5.0 0152, 

0,-0 0026, 

5.-0 0207. 

10,-0 0392, 

15,-0 05*. 

20,-0 0772, 

25.-0 0967, 

30.-0 1165, 

35,-0 1367, 

03, 

-10,00327, 

-5,0 0152, 

0,-0 0026. 

5, -00207, 

1 0.-0 0392, 

15,-0058, 

20,-0 0772. 

25.-0 0967, 

30,-0 1165, 

35,-0 1367, 

06, 

-10,00379, 

-5,0 0175, 

0,-0 0032, 

5.-0 0242, 

10,-00455, 

15,-0 0672, 

20,-0 0*92. 

25.-01117. 

30.-0 1345. 

35.-0 1578. 
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095. 

-10,00592. 
-5,00248. 
0,-0 0098, 


20.00108, 
30,0 0056. 
40, -0 0095, 
50,-0 0303, 


5,-0 0447, 60, *0 0566, 

10,-008, 70,-0 0923, 


15,-0 1157, 


15 , 


20,-0152, 


-20,-0 0494, 


25.-0189, 

30,-02266, 

35,-0265, 

11. 

-10,0 1043, 
-5,00532, 
0,00021. 

5, -0 049, 
10,-01001, 
15,-01512, 
20,-0 2023, 

25.4) 2534, 
30,-0 3045, 

35.4) 3555, 
2. 

-20,-0 0274, 
-10, -0 021, 
0,-00211, 
10,-00258, 

20.4) 0296, 
30, -0.0296, 

40.4) 0318, 

50.4) 0584, 

60.4) 0775, 

70.4) 1124, 

3, 

-20,-00364, 

-10,-00282, 

0,-00197, 

10.4) 0134, 
20,-0 0103, 
30,-00111, 
40,-00187, 
50,-0 0342, 
60,-00591, 
70,-0 093, 

4, 

-20,4)0419, 

-10,4)0322, 

0,4)0192, 

10.4) 0076, 

20.4) 0014, 

30.4) 0026, 

40.4) 0122, 

50.4) 029, 

60.4) 0556, 
70,-0 0878, 

6, 

-20,4)0468, 

-10,4)036, 

0,4)019, 

10.4) 0024, 

20.0 0063, 

30.0 0039, 

40.4) 0085, 

50.4) 0277, 
60,-0 057, 
70,-0 0867, 
8. 

-20,-0 0486, 
-10,-0 0377. 
0.4)019, 

10.4) 0001, 

20.0 0092, 

30.0 0056, 

40.4) 0083, 
50,-0 0284, 

60.4) 0572, 
70.-0 0873, 
10 , 

-20,-0 0493, 
-10,4)0385, 
0.-0019, 

10.0 0011 . 

20.0 0104. 

30.0 0059, 

40.4) 0088. 

50.4) 0295, 

60.4) 0577, 

70.4) 0902, 

12 . 

-20.4)0495. 

-10,4)0389, 

0,4)019, 

10,00018, 


-10,-0 039, 

0,-0 0189. 

10.0 0023, 

20.0 011, 

30.0 0051, 

40,-00104, 

50,-0 0312, 

60,-0 0573, 

70,-0 0965, 

20 , 

-20,4)0489. 

-10,4)0397. 

0,4)0193, 

10.0 0029. 

20.0 0111 , 

30.0 0043, 

40,-0 01 19, 

50,-0 0324, 

60.4) 0583, 

70.4) 1045, 

25, 

-20,-0 0494, 

-10,-0 0398, 

0,4)0193. 

10,0003, 

20.0 0113. 

30,00039, 

40.4) 0125. 

50.-0 0334, 

60.4) 0599, 

70,-0 1138, 

S 

PStafc table~6hcdbdli ,2,61ulpha ,6hmach .10,15,8*1, 
0, 

-10,4)021666, 

-5,-0 021766, 

0,-0 021666, 

5.4) 021666, 

10.4) 0228, 

15,-0 0228, 

20.4) 0228, 

25.4) 0227, 

30.4) 0228, 

35.4) 0228, 

03, 

-10,4)021666, 

-5,4)021766. 

0,4)021666, 

5, -0 021666, 

10.4) 0228, 

15,-0 0228, 

20.4) 0228, 

25,-0 0227, 

30.4) 0228, 

35,-0 0228, 

06. 

-10,4)020004, 

-5,-0 020104, 

0,4)020104, 

5.4) 020004, 

10,-0 034, 

15,-0 034, 

20.-0 0341, 

25,-0034, 

30.4) 034, 

35.4) 0341, 

095, 

-10,4)051882, 

-5.-0.051882, 

0,4)051882, 

5.4) 051882, 

10.4) 0619, 

15.4) 0619, 

20.4) 0619. 

25.4) 0619, 

30.4) 0619, 

35.4) 0619, 

I 1. 

-10,4)056865, 

-5,4)056865, 

0,4)056865, 

5.4) 056865, 

10,-0 0676, 

15.4) 0676. 

20,-0 0676. 
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25.-0 0676 . 
30,-0 0676 . 
35 .- 00676 . 

2 , 

- 20,-0 0273 * 4 , 
- 10.-0 02 * 7 * 4 , 
0,-0 0292 * 4 . 

10 . 4 ) 02 * 7 * 4 , 
20 , *0 0302 . 
30,-0 027 *. 
40 ,- 0 . 0246 . 
50 ,- 0 . 0207 , 

60, -0 0161 . 
70,-0 on . 

3 , 

- 20 .- 0019262 . 
- 10 ,- 0019962 . 
0,-0 020162 . 
10 ,- 0 . 019962 , 
20 ,- 00147 , 
30 ,- 00135 . 

40 , - 001 19 , 
50,-0 01 , 
60 ,- 0007 *. 
70,-0 0054 . 

4 , 

- 20 ,- 0011751 . 
- 10,-0 Ot 2051 , 
0.-0 012251 . 
10,-0 012051 , 
20,-0 00 * 5 , 
30 .- 00077 , 
40.-0 0069 . 
50,-0 005 *. 
60 ,- 00045 . 
70,-0 0031 , 

6 . 

- 20.-0 006393 , 
- 10 ,- 0006593 , 
0,-0 006593 , 
10 .- 0006593 , 

20 . 4 ) 003 *, 
30.-0 0034 , 
40 ,- 0003 . 
50,-0 0026 , 
60,-0 002 , 
70,-0 0014 . 

5 , 

-20.-0 002, 
-10,-0 0022. 
0.-0 0022, 
10,-0 0022, 
20.-0 0021. 
30.-0 0019 , 
40,-0 0017 , 
50,-0 0014 . 
60 ,- 00011 , 
70,-0 0008 , 

10 , 

- 20, -0 0013 , 
- 10,-0 0014 . 
0,-0 0015 . 
10,-0 0014 , 
20.-0 0012 . 
30 .- 00011 . 
40,-0 0009 , 
50.-0 0008 . 
60,-0 0007 , 
70,-0 0005 , 
12 . 

- 20,-0 0007 , 
- 10 . 4 ) 001 , 
0 ,- 0001 , 
10 ,- 0 . 001 , 
20 ,- 0 . 000 *. 
30,-0 0006 , 
40,-0 0005 . 
50,-0 0005 . 
60.-0 0004 . 
70,-0 0003 , 
15 , 

- 20,-0 000 *. 
- 10 , 4 ) 0006 , 
0 . 4 ) 0006 , 

10 . 4 ) 0006 , 

20 . 4 ) 000 *. 

30 . 4 ) 0001 , 

40 . 4 ) 0001 . 

50 . 4 ) 0002 . 

60 . 4 ) 0002 , 

70 . 4 ) 0002 . 

20 . 

-20.-0 0002, 


- 10 , 00001 , 
0,0 0002, 
10,0.0002, 
20.-0 0002, 
30.-0 0003 . 

40 . 4 ) 0003 , 
50,-0 0002 . 
60 , 0 , 

70 . 4 ) 0001 , 
25 , 

- 20 . 0 , 

- 10 , 0 . 

0 , 0 , 

10 , 0 . 

20 , 0 , 

30 . 0 , 

40 . 0 , 

50 . 0 , 

60 . 0 , 

70 , 0 . 
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A.2. SINGLE ELEMENT BIMESE: LEO 


PSSEARCH 

* 

C LOX/LH2 SINGLE BIMESE W/RS2I00 Engine (reviled 1/1/99) 

C WB-003-C Bimese 0 Klb to 1 00 nmi 
C 2* $ DEGREES INCLINATION 
C LATEST WB-OO I AERO ADDED 5/6/97 

* 

c 

SRCHM - 4, / PROJECTED GRADIENT METHOD 
c IPRO •-!,/ PRINT FINAL TRAJECTORY ONLY 
IOFLAG - 0, / ENGLISH IN * ENGLISH OUT 
CONEPS(l)“ 89 9, / CONVERGENCE TOLERANCE 
CONEPS<2) - 00001, / MIN ALLOWED PERCENTAGE CHANGE 
CONEPSO)- 00001,/ 

CONEPS(4) - 00001, / 

CONEPS(5) - 00001,/ 

PCTCC- 001, / 

IDEB-0. / No Detail 
C 

•••••••••••• ••• 

C OPTIMIZATION VARIABLE 

* 

C 

OPTVAR - 6HW EIGHT. / MAXIMIZE FINAL BURNOUT WEIGHT 
OPT - I, / MAXIMIZATION FLAG 

OPTPH - 500, / PHASE AT WHICH WEIGHT IS MAXIMIZED 
WOPT - 000003, / OPTIMIZATION VARIABLE WEIGHTING 
MAXITR -40, / MAXIMUM ITERATIONS ALLOWED 
C 

* 

C 

C CONSTRAINT VARIABLES 

C 

* * •••••••••• 

c 

NDEPV - 6, / NUMBER OF DEPENDENT VARIABLES USED AS 
CONTROLS 

DEPVR - ShgdaJt, 6HGAMMA1, 3H1NC, 5HXMAX1, SHXMAX2, 5HXMIN2, 
/ INSERTION ORBIT CONDITIONS + CONSTRAINTS 
DEPVAL - 303805 0, 0 0, 28 5, 1000 , 379165 , -379165 , / INSERTION 
ORBIT VALUES + CONSTRAINTS 
DEPTL-500, I, I, 1 100, 100, / DESIRED ACCURACIES 
DEPPH - 500, 500, 500, 500, 500, 500. / EVENTS 
IDEPVR - 0, 0. 0, I, 1,-1, 

C 

* * 

c 

c 

C CONTROL VARIABLES 

C 

* * 

c 

nindv - 16, 

indvr - I5*6hpitpc2.3heil, 

indph - 15,18,20.30.40.60.65,70,75,80,85,90.95.100,105.10, 
u- -1.1121 10516660E-KW, 9 59606 1841 384E-02. 2 4057203 1201 2E-02, - 
8 6478607 1 5809E-0 1 , -8 524468987230E-0 1 . 

-7.365I35456373E-01, -6 185493375 13 IE-02, *1 78091 56563 36E -01. - 

2 827569253803E-01, -2 1 171 5323 1269E-01 . 

-I 9777225994 79E-03, -6 98 1303256862E-02, - I 94965587895 IE-02. - 

3 901 7533 14859E-01 . - I 333034849750E-01 , 

8 62I963895029E+01 , 

u- -9.63 1 1 4661 2932E-0 1 , 1 206853378424E-01, 3 371 722786330E-02, 

4 335495035323E-OI, -4 798 1 4 900764 2 E -01, 

-4 000001681097E-01, -4 34802005 3465 E-02. - 1 41 54 1 95 1 5 1 40E-0 1 , 

-4 0362450 17758E-01, -9 91 599864371 2E-02, 

-1 85 12988967 1 7E-03, -6 98 1 303256862E-02. - 1 94965587895 1 E-02. 

*3 9017533 14859E-01, -1 33 3034849 750E-01, 

8 .62 1 963895029E+0 1 , 

u- -2 04303487 1238E+O0, 1.22536924 18 72E-01, 2 I 29087327276E-02, - 
2 688018798002E-01, -8 078862072536E-01, 

■6.01 18868 94 844 E-01, -7 6066683 13 577E-02, -2 2365061 1 1735E-01, - 

4 402526852204E-0 1,-1 4060 1 565487 1 E-0 1 , 

-I 878724563661 £-03, *1 0732455001 3 5E-0 1, -2 122 784772 1 97E-02, - 
3.4841 24575335E-01 , - 1 235092248 1 24E-01 , 

• 62 7 362883464 E+0 1 , 

u- -2 048482830768E+00. I 270951672930E-01, I 952917063025E-02, - 

2 808866024 1 97E-0 1, -8 9029623 26908E-01, 

-7 600675483939E-01, -1 057582836078E-0I, -2 6743 7255295 1 E-0 1 , - 

3 I88541758559E-OI, -I 783906646 563 E-Ol, 

• I .895282842 109E-03, -I .3960423006S2E-OI, -2 200578 107558E-02, - 
3.73401 7639366E-0I, -I 230747709286E-0I, 

8.627366205750E+01, 

u- -1 925681 15271 1E+00, I 27591 1 55494 1 E-OI, 1 952723724261 E-02, - 
2 79124S629003E-01. 4 70 10055428 1 4E-0 1, 

-7.1 9798785 9824E-01, -1. 0368024448 56E -01, -2 6 568 3 1 800408 E-0 1, - 
3.293637O48943E-01. -1.871700580567E-OI, 

-1 896853573847E-03.-1 4680662804 1 8E-0 1 . -2 22093023241 7E-02. - 

4 36 1 27076523 IE-01, -1 3 28583934429E-0 1 . 


8 62526626 1 937E+01, 

u- -I 897293827039E+00, 1 2777755 1 8953E-0 1, 1 95281 1 7608 79E-02, - 

2 775129774345E-OI, 4 517441085340E-01, 

-7 04 50 1 443 1 499E-0 1, -I 0371 84763853E-0 1 , -2 674 1 1 5 1 24 1 93E-01 , - 

3 3582245 14068E-01, -I 900500823 104E-01 , 

-I 897 1 88763099E-03, -I 480803 1 8023 5E-01, -2 223470725970E-02. - 

4 43732404473 IE-01, -I.336372477044E-01, 

• 6244O0839367E+OI, 

it- - 187707398723 1 E+OO, I 2793601 7561 5E-01, I 95 3 28344 2428 E-02, - 

2 757516305779E-01, 4.349255442527E-01, 

4.872640468 1 40E-01, -1.034497473644E-OI, -2 679953020693 E-01, - 

3 39931 8594345E-01, -1 921090903 175E-OI, 

-1 89742999203 5E-03, -1.490770425276E-OI, -2 225582901 147E-02, - 

4 510385562812E-01, -1.344942922373 E-01. 

8.623542756604E+0 1 , 

PERT • 16*1 0E-4, /PERTURB 

5 

••• * ••••*••• 

C 

C TRAJECTORY SIMULATION INPUTS 

C 

* 

PSGENDAT 

TITLE - OH* LOX/LH2 BIMESE 28 5 deg delivering 60 Klb *. 

EVENT - 1 0,/ FIRST EVENT 

WGTSG - 1 390000 , / 1NITAL LIFTOFF WEIGHT witfi ADDED FLUIDS 

FESN - 500, / FINAL EVENT NUMBER 

MAXTIM - 1000 , /MAXIMUM TIME 

ALTMAX - 10000000 , / MAXIMUM ALTITUDE ALLOWED 

ALTMIN - -5000 , / MINIMUM ALTITUDE ALLOWED 

PRNC - 0, /letup plot file 

MONX - 4HDYNP, 4HFAZB, / MONITOR Q & NORMAL FORCE 
PRNT(97) - 5HXMAXI, 5HXMAX2, 5HXMIN2, 6HALTTTO, 4HMACH, 
4HDYNP, 

4HAMYB. 5HTTMYB, 5HPSTOP, / 

IGUID( 1 ) - I, / INERTIAL EULER ANGLE OPTION 
IGUID(2) - 0, / SAME STEERING OPTION FOR ALL CHANNELS 
IGUID<4) - 1, / CUBIC POLY WITH CONSTANT TERM SET BY INPUT 
PITPC1 1 ) - 0 , / INITIAL PITCH RATE 
YAWPC( 1 ) - 0 . / INITIAL YAW RATE 
ROLPC( 1 ) - 0 , / INITIAL ROLL RATE 

NPC( 1 ) - 2, / PRINT CONIC BLOCK AT EACH PHASE CHANGE 
PINC - 20, / PRINT INTERVAL EVERY 20 SECONDS 
NPC(2) - 1, / FOURTH ORDER RUNGE KUTTA 
DT - 1. / INTEGRATION STEP SIZE 

NPC(3) - 4, / PLANET RELATIVE INPUT ON VELOCITY VECTOR 
VELR - I , / INITIAL VELOCITY AT LAUNCH PAD 
GAMMAR - 90, / INITIAL PLANET RELATIVE FLIGHT PATH ANGLE 
AZVELR - 0 , / INITIAL AZIMUTH ANGLE OF VELR 
azl -90, 

NPC(4) - 2, / PLANET RELATIVE INPUT ON ALTITUDE VECTOR 
ALTITO - 0 , / INITIAL ALTITUDE 
GDALT - 0 . / INITIAL ALTITUDE 
GDLAT - 28 5, / LATITUDE OF LAUNCH SITE 
LONG - 280 0, / LONGITUDE OF LAUNCH SITE 
LONG I - 280 0. / LONGITUDE EAST OF PRIME MERIDIAN 
NPC<5) - 5, / 1976 US STANDARD ATMOSPHERE 
NPC(6) - 0, / NO WINDS 
NPC(7) - I. / ACCELERATION LIMIT SET 

ASMAX -30./ UMIT ACCELERATION THROUGHOUT TRAJECTORY 
LREF - 162 6, /BODY LENGTH FROM NOSE TO BASE 
SREF - 3500 7, /REFERENCE AREA FOR ORB I TER ONLY 
NPC(8) - 2, / CL, CD A CM TABLE INPUTS 

NPC(9) - I, / ROCKET ENGINE WITH THRUST TABLE AND ISPVAC 
NPC( 10) - 0, / STATIC TRIM IN PITCH 

GXP( 1 ) - 160 2. /X-LOCATION OF ENGINE GIMBAL MEASURED FROM 
NOSE 

ITRIM - I, /TRIM WITH ENGINE DEFLECTIONS ONLY 
IENGT - 1, / CALCULATE ENGINE INCIDENCE ANGLES FROM STATIC 
TRIM EQUATIONS 

NPC( II ) - 0. / NO FUNCTIONAL INEQUALITY CONSTRAINTS 
NPC< 1 2) - 3, / CALCULATE CROSSRANGE AND DOWNRANGE BASED 
ON ORBIT 

ALTREF - 1000. / 100 NMI CIRCULAR MILE REFERENCE ORBIT 
AZREF - 0 , / AZIMUTH REFERENCE 
NPC(1 3) - 0, / DO NOT JETTISON PROPELLANT 
NPC< 14) - 0, / NO HOLDDOWN BEFORE LIFTOFF 
NPC( 1 5) - 0, / DO NOT CALCULATE AEROHEATTNG 
NPC< 16) - 0, / USE OBLATE EARTH GRAVITY MODEL 
c NPC( 1 7) - 2, / JETTISON BOOSTER 
NPC< 18) - 0. / DO NOT TERMINATE TRAJECTORY 
C NPC( 1 9) - 0. / DO NOT PRINT INPUT CONDITION SUMMARIES 
NPC(20) - 0, / DO NOT USE ANY SPECIAL DT CALCULATION 
NPQ21 ) - 0, / DO NOT CALCULATE FUEL AND OXIDIZER WEIGHTS 
AND VOLUMES 

NPC(22) - 0. / DO NOT CALCULATE THROTTLING 
NPC(23) - 0, / DO NOT COMPUTE VELOCITY MARGINS 
NPC(24) - 0, / DO NOT COMPUTE ANY PARAMETER INTEGRALS 
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NPC<25) - 3. > COMPUTE VELOCITY LOSSES AND PRINT AT EACH 
TIME 

NPC(26) - 0, t NO SPECIAL AEROHEATING CALCULATIONS 
NPC(27) - 0, / DO NOT INTEGRATE ENGINE FLOWRATES 
IWPF( l ) - I, /INTEGRATE BOOSTER FLOW RATES 
NPC(2S) " 0, / TRACKING STATION OPTION NOT USED 
NPC<29) - 0. / DO NOT COMPUTE VACUUM IMPACT POINTS 
NPC(30) - 0. / USE N-STAGE VEHICLE WEIGHT MODEL 
IENGA(I) - I , / THROTTLE ORB ITER 
IENGMF( I) - 1. / ORBITER ROCKETS ARE ON 
|WDF( I ) * 1, / CALCULATE ORBITER FLOW RATES AS A TABLE 
LOOKUP 

NENG- I./ THE NUMBER OF THRUSTING ENGINES 
NPC(3|) - 0, / NO VERNAL EQUINOX, SUN SHADOW, SUN ANGLE 
CALCULATIONS 

NPC(32) - 0, / NO PARACHUTE DRAG 

NPQ33) - 0, / NO ARC LENGTH CALCULATIONS 

NPQ34) - 0, / NO KEPLER1AN STATE CALCULATIONS 

NPC(3S) - 0. / NO SENSED VELOCITY INCREMENT CALCULATIONS 

NPC<36) - 0, / DO NOT ACTIVATE SUNLIGHT OPTIONS 

NPCQ7) - 0. / DO NOT ACTIVATE DATE OPTION 

NPQ3I) - 0. / NO ATMOSPHERIC TURBULENCE 

NPC(39) - 0. / NO ATMOSPHERIC GUSTS 

NPQ40) » I,/ 

$ 

PSTBLMLT 

TVCIM * 3.3351. / Thiuet multiplier for Orbiter 
WD1M - 3 3351, /Flow me multiplier for Orbiter 
AEIM - 3 335 1. / Exit im multiplier for Orbiter 
$ 

PSTAB 

TABLE - 5HTVC1T, 0,53 1 500 , / VACUUM THRUST OF Orbiter engine 
S 

PSTAB 

TABLE - 4HWDIT.0.1 199 77, / FLOW RATE OF Orbiter engine 

S 

PSTAB 

TABLE - 4HAE1T. 0,33 5 5S, / EXIT AREA OF Orbiter engine 
S 

PSTAB TABLE - 4HXCGT, l .6HWE1GHT.2. 1 ,- 1 . 1 , 

2006450,123 5 1,1 34*70,11 5 50, 

S 

PSTAB TABLE-5HXREFT.0.IO9 90. 

S 

PSTAB TABLE- 5HZREFT.0.0, 

S 

•include Thome/aadl l/Jtooley/poet/BlMESEybtroeae aero 
ENDPHS- I. 

$ 


event-1 5,critr-4htime,value-07,endph*-l , 

iguid(4>“0, 

S 

ISgendat 

event- 1 t.cri tr-4hti me, value- 1 7, end phi- 1 , 
S 


event-20,critr-4htime,value-30,cndphs- 1 , 

S 


event- 30, cntr-4hbm«,value-40,endphi- 1 , 

S 


event— 40,cri tr-4hbme, vahie-60.endphs- 1 , 

S 


event^O, critr-4hrime,value-S0.endphs- 1 , 

S 


event-«5,cntJ-4htime,value-l 10.endphi-l, 


event- 70,cno-4hti me, value- 1 50,endphs- 1 , 

S 


ISgcndet 

event-7 5 ,critr-4htime,vahie- 1 70,endp*»- 1 , 

S 

ISgendat 

event-SO.cn tr-4btt me, v*k*e-200.»><lphr- 1 , 

S 

ISgendat 

•vent-8 5 ,cntr-4ht) me, value- 230,endph»- 1 , 

S 

ISgendat 

event-90,critr—4hti me, value- 260,endphj- 1 , 

S 

ISgendat 

event-95 ,critr-4h6me,value-2S0,eodph*- 1 , 

S 

ISgendat 

event- 1 00,cri tr-4hbme. value- 300,endptw- 1 , 

S 

ISgendat 


event- 1 0 5 ,cri tr*4hti me, value- 3 20,endphi- 1 , 

S 

ISgendat 

c this event marts injection into orbit 

event-500,cntr-4hvdi,value-25*55 000, 

ENDPHS - 1. 

ENDJOB- 1. 

ENDPRB - I, 

S 
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A.3. THRUST-AUGMENTED BIMESE: LEO 


pSaearch 


c LOX/LH2 Single Bimese w/solids 
c wb-003-c bimete 
c latest wb-001 aero added 5/6/97 

..•••*••••* 

c 

c Ascent of: 

c Single Bimeae w/sdids to a 100 nmi x SO nmi x 21.5 degrees inclination 
c Jeffrey Tooley. Sept 1999 
c 

c Optimization Variables 


archrn - 4, / projected gradient method 

ioflag* 0, /engjish in - english out 

aptvar - 6h weight, / maximize final burnout weight 

opt - I. / maximization flag 

optph * 500, / phase at which weight is maximized 

wopt * .000003, / optimization variable weighting 

maxitr * 40, / maximum iteration* allowed 

coneps(l)-S9 9, /convergence tolerance 

coneps(2) - 00001. / min allowed percentage change 

conepa(3)- 00001,/ 

ooncpe(4) * 00001, / 

cooeps(5)- 00001,/ 

pace-. 001,/ 

idcb-0, / no detail 


c Dependent Variables (constraints) 


ndepv -II,/ number of dependent variables 

depvr - Shgdalt, 6hgammai, 3hinc, Shxmaxl, Shxmax2, 5hxmin2, Shvdi , 
5hxmax4, ShxmaxS, 5hxmin5,5hxmin6, 

depvaJ - 303805, 0 0,2* 5. 1000,379165,-379165, 25*55, 3,1 0001, 

0 19999, 0 19999, 

deptl - 500, 1, 1, 1 , 100, 100, 100, 1, 001, 001. 001,/ desired accuracies 
depph - 500, 500, 500, 500, 500, 500, 500, 500, -14, 500, 500/ everts 
idepvr • 0, 0, 0, 1. 1, -1,0, I. I. -1.-1. 


c Independent Variables 


nindv - 15, 

indvr - 1 2 *6hpi tpc2,3hazl ,6hwstpd2,6heupc I , 
indph - 1 5,1 6,20,30,40,60,6 5, 70,7 5,80, 85,90, 1 0, 1 0, 1 4, 
pert - lSM 0e-4, / perturb 


c initial guesses for independent variables (last set only) 


/za gem motor 

/ 

Ais for T - 1250000 and tb - 50 

u- -1.01661 14830S3E +00, -8 73492I04O992E-OI.-7 298162085384E-OI,- 
1 837842834798E-0I.-I 200160273233E-01, 

-2 . 7093 52 765022E-0 1,-1 056 161*1 3457E-0 1 , -3 *99*92 1 28789E-01 , - 
1 663 1 5262939 1 E-0 1 , - 1 1 560002062 2 6E-01 , 

-9 IIS640989624E-05, -* I00390628086E-03, I 630*739277 19E +01, 

11 74*42 158622E+06, 7 75283S032258E-0I, 

/ 

AisforT- 1250000 and tb - 75 

u- -1 016487 161 331 E+OO. -7 661 696981 75 3E-0 1, -9 30S445422486E-01. - 
1.812502263332E-0I, -l 161 50403 564 2E-01, 

-2 . 5547389623 28E-0 1 , -I 0127*4101 59IE-01, -3 42068988 5 57 IE-01, - 
1 5*0161 341920E-OI, -I 114993 191605E-01, 

-9.1 1 *4530422 50E-05, -8 0354601 1W65E-03, 8 63 20024676 50E +01 , 

I IM216263756E+06, 9 0900369*002 5 E-0 1 , 

/ 

/us for T - 1250000 and tb- 100 

u— *1 013127II4I02E+00, -7 62 13 87946 283 E-0 1, -9 336093040501 E-0 1, - 
1 822677980694E-0 1 ,-1 1 706772 79692E-0 1 , 

-2.61 50085607 14E-01, -1 025580027732E-01, -3 569933474 1 43E-01, - 
l 5903I6364633E-01.-1 11 9607963 05 5 E-0 1 , 

-9 1 184S2452200E-05, -* 04363 1635838E-03, 1 632301 1 12353E+01, 

I 1M248485805E+06, 9.4I5103497269E-01. 

/ 

/us for T- 1250000 and tb - 125 

u- -1 2305339041 97E+00, -5 740527778979E-01, -7 399*4206 1065E-01, - 
5 3122747012I9E-OI.-1 49551 0*95 9 19E-01, 

-1 720920461 344E-01, -1 99601 1203752E-01, -2.384473650279E-01 . - 
1 59900537 1693E-01, -9 *67*6895 823 0E-02, 

-9.1 171468907 1 0E-05, -7 617017970988E-03 , 8 63 209504694 9E +01, 

I 1911197 1231 3E+06, 9 9*43 7*0323 58E-0 1, 

/ 

Ass for T - 1 500000 and tb - 50 

u- -I 021 634553664E+00, -8 9522492852 19E-0 1,-7 0**01 6202*25E-01, - 
1 SI301838I815E-01, -1 1*8169964731 E-0 1. 


-2 668948 103490E-01. -I 052673956727E-01, -3 *8458 56078 56E-01, - 
I 660365 592 1 92E-0 1,-1 1 54965926578E-0 1 . 

-9 118623 132244E-05. -8 091639257352E-03, 8 6308 563028 53 E+Ol, 

1 1795 1361 1878E+06. 6 723I80739283E-01, 

/ 

/us for T - 1 500000 and tb - 75 

u- -1 024092069 128 E+OO, -7 44421 1046506E-01, -9 75 1303063 896E-01, 

1 75928577782 IE-01, -I I38439419397E-01, 

-2 451 072795 132E-01, -9 971374051328E-02, -3.285 1 67753 522E-01. - 

1 5696571 703B5E-0 1,-1 109599670444E-0 1 , 

-9. 1 1*4 1 6294708 E-05 , -* 02467863 S726E-03, 8 63200033 7253E+01, 
1II9412235I05E+06, *454451 7*3*91 E-0 1. 

/ 

Ais for T- 1500000 and tb - 100 

u- -l 077524926498E+O0. -64944I1690905E-O1, -1 03000747 1972E+O0, 

2 5390629*6 1 69E-01 , -I 54 1 083409W7E-0 1 , 

-1 4*14294801 79E-0 1, -1.871 145333546E-01, -2 4*7*0449393 7E-01, - 
1 98565*66 1540E-01, -1 001 712559494E-01, 

-9 11 7326948552E-05, -7 673542691 382E-03, * 63 290 560977* E +01. 

1 19S670608738E+O6, 9 808437135977E-01, 

/ 

/us for T- 1 500000 and tb - 125 

u- -l 083483 70481 7E+00, -6.299342591 12OE-01, -7 5793 266*3 569E -01, 
4 3480360061 34E-01. -1 557763536568E-OI, 

-I 792066*062 23E-01, -2 257504916082E-01. -2 637245432764E-0I, - 
1 98 1 2 1 3733257E-01 , -9 97271 7688620E-02, 

-9 1 1 72 19021 4 1 2E-05, -7 639 1 70094409E-03, * 63 1 032487455E+O I . 

11 9893 740555 1E+06, * 690272270267E-OI, 

/ 

/us for T - 1750000 and tb - 50 

u- -1 030890331 1 52E+00, -8 819237775802E-0I, -6 04 1 52764727 1E-0 1 . 

3 058207695465E-O1, -1 0*8 1 4249*3 72E-01, 

-1.71 2496699280E-01 , -2 44 1 458892383E-0 1,-3 1 02 5 7440009 7E -01 . - 
1 83760191 3347E-0I.-1 376444393257E-OI, 

-9 11 84543 31 076E-05, -8 10033*6086 11 E-03, 8 630230586672E+O1. 

1 1 83879 192962E+06, 4 7532231 I648IE-01. 

/ 

/us for T - 1750000 and tb - 75 

u« -1 012561089280E+00, -* 109375803400E-01, -6 883514888500E-0I, 
3 082536426446E-0 1,-1 0747 1 3990997E-01 , 

-1 6958779 1 6238E-01 , -2 4 1 461 2908 1 23E-01 , -3 066562 1 63292E-0 1 , - 
1 770676 798990E-01.-1 38594204S407E-01, 

-9 1 18512509890E-05, -8 133636325475E-03. 8 629974 593443 E +01, 

1 194725675 189E+06, 5 686873842838E-01, 

/ 

AisforT - 1750000 and tb - 100 

u* .11 1 9644404 139E+00, -7.2992567420 1 4E-01, -6.290610811 593E-01, 

2 9982622 79759E-01, -1 089327620667E-0I, 

-1 71 5299868060E-0 1 . -2 672650478202E-01 , -3 62004 1 7 1 44 77E-0 1 , - 
1 939629 147078E-01, - I 49360250829 1 E-0 1 , 

-9 118639508832E-05. 4.2216*226273 IE-03, I 62* 2 5492054 2E +01, 

1 201342294671 E+06, 5 33*83 1737613E-01. 

f 

AisforT- 1750000 andtb- 125 

u- -1 03449062 1 985E+00. -6 296869070388E-01. -7 476293 309793 E-0 1, 

3 B66599777706E-01, -I 520678057466E-01, 

-1 7995673 20494E-01. -2 4313961 18968E-01, -2 8907501 53 747E-0 1. - 

2 5823 1 3830897E-0 1,-1 0450332 1 1 370E-01 , 

-9 1 1 7250966781 E-05, -7 642675498768E-03, 8 6295 774 78 56 5 E +01 . 

1 20672946463 7E+06, 7224I3691480SE-01, 

/ 

Ais for T - 2000000 and tb * 50 

u- -1 041 700574636E+00, -9 1 1 3060589251 E-0 1 . -6 0* 1 74842 795 5E-0 1 , 

3 0O4653S60676E-OI, - 1 07445*3 53958E-0 1 . 

-1 663580530485E-01, -2 331437685358E-01, -2 920105955044E-01, - 
1 S12598477304E-0I, -I 36278232I997E-01 , 

-9 1 1 *40 1 60720 1 E-05, -* 0*3 5*76703 84E-03, * 63065752054 1 E+0 1 , 

1 1 88357998283E+06, 4 579755706971 E-0 1 , 

/ 

Am for T - 2000000 and tb - 75 

u- -9 93S943556588E-OI, -7 I96945397596E-01. -* 3*3 3 2296445 7E-01, 

3 5 1 606978 5664E-0 1,-1 3930930*6263 E-0 1 . 

-1 526 15*53 1641 E-0 1. -I *9 II 47042*51 E-0 1,-2 250*1 4302084E-01, - 

2 42303 1 02407BE-0 1,-1 01408*8 1 6368E-0 1 . 

-91171 07480 505E-05, -7 61112 I03M06E-03, * 630657 1 20469E+0 1 , 

1 201I22540204E+06, 5 877265285799E-01, 

/ 

/ua for T - 2000000 and tb * 100 

u- -1 0253 59*06951 E+OO, -6 6364925 1 9973E-01,-* 18*55250 1 225E-01, 

3 676142214755E-OI, -1 440I072I9982E-01, 

-1 603*95421 557E-OI, -2 1 17I95648724E-01, -2 492677779574E-01, - 

2 4721583 13 723E-01, -1 022I2*4*46*4E-0I, 

-9 11716744 5845E-05, -7 626472284749E-03, *6303470991 1 1 E+O I, 

I 210*351 11653E+06, 6 429308137273E-0I, 

/ 

An for T - 2000000 and tb - 125 

u- -l 0244647 13350E+00, -6 31 1 56342447 1 E-0 1, -7 6986*7*83663E-01, 

3 *366673 59B59E-0 1, -I 503 17362975 1 E-0 1, 
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.| 772M9669686E-0I, -2.362454461379E-01. -2 *2707*427946E-01. - 
2 S46**4I12394E-0I, *1 0324O51122WE-OI. 

J9 U72316*6I69E-05.-7 63499I3374*9E-03. 1 629733 1 *3355E+0 1 , 

I 2I4971202342E+O6. 7 023*25295 1 25E-01. 

/ 

Au for T - 2000000 tad lb - 100 

u- -1 0*34*3 704* 1 7E+O0, -6 299342391 120E-01, -7 5793 26683 569E-01, - 
4 34*0360061 34E-01, -1.55776353656SE*01, 

•1 792066*06223 E-01 , -2.257S04916082E-01, -2.637245432764E-01, - 
» 981213733257E-01, -9 97271 7688620E-02, 

-9 1 172I9021412E-05, -7,639 170094409E-03, 1 631032487455E+01. 

J 19*93740555 1 E+06, * 6902722 70267E-01, 

S 


e BEGIN EVENT U STING * Vtpctary amuls&on inputt 


tide - Oh* LOX/LH2 Bimese 2* 5 deg delivering 60 klb *. 
event - 10, / fin* event 


^cOO -0,/ no functional inequality conraiints 
feen - 500. / fine! event number 
npc(1) - 2. / print conic Mock tt each phase change 
pine - 2, / print interval every 20 seconds 
npc(2) • 1, /fourth order runge kutta 
dt - I, / integration step ex* 
npc(20) - 0, / do not use any *>eaal dt calculation 
npc(25) - 3, / compute velocity I oases and print at each time 
maxtm - 1000 . / maximum time 
litmax - 10000000 , / maximum altitude allowed 
altmin - -5000 , / minimum altitude allowed 
pmc - 0, /setup plot file 

esonx - 4hdynp, 4hfazb, Shalpha. 4hasmg. 4heta . 4hetaJ./ monitor q k 
normal force 

pmt(97) - Shxmax 1 , ShxmaxZ Shxmin2, 6haltno, 4hmach, 4hdynp, 3hwdl, 
3hwd2, 

pn»t(105^ Jhwprpl, 5hwprp2, 4hamyb, Shttmyb, 5hxmax5. 5hxmin5, 
5hxmin6y 


c Guidance setup 


igutd(l ) - 1,/ inertial euler angle option 
iguid<2) - 0, / seme steering option for all channels 
iguid(4) - 1 , / cubic poly with constant term set by input 
pitpc(l) - 0 . / initial pitch rate 
yawpc< 1 ) * 0 / imtial yaw rate 
rolpc<l) - 0 . / initial roll rate 


c Earth relative Components 


npc(3) - 4, / planet relative input on velocity vector 
vdr - .1 , / initial velocity at launch pad 
gammar - 90., / initial planet relative flight path ingle 
arvelr - 0 ,/ initial azimuth angle of vdr 
ad -90, 

npc(12>- 3, / calculate croaartnge and downrange baaed on orbit 
altref- 100 0, / 100 nmi circular mile reference orbit 
azref - 0 , / azimuth reference 


c 

c 

c Earth Launch Site 

— — 

npc(4) - 2. / planet relative input on altitude vector 
alrno - 0 , / initial altitude 
gdaJt - 0 , / initial altitude 
pdlat - 2* 5, /latitude of launch site 
long * 2*0 0, / longitude of launch die 
longi -2*0 0,/ longitude east of prime meridian 


c Earth etmoaphere inputs 

c — — 

npc<5) - $./ 1976 us standard atmoaphert 
npc(6) - 0. / no winds 

c 

c 

c Earth Gravity 

c 

npc(!6) - 0. / use oMalc earth gravity model 


c Limit acceleration 


npc(7) - 1 , / do not limit eocderatkm 
asmax -30,/ limit acceleration throughout trajectory 


c 

c 


Aerodynamic and Trim Data (Tables down below) 


npc(8) - 2, / d, cd k cm table inputs 
lref - 1626. /body length from nose to base 
irtf - 3500 7, /reference tree for orbiter only 
npc(IO) -0, /static trim in pitch 

gxp( 1 ) - 160.2, /x-locabon of engine gjmbai measured from noee 

itrim - 1, /trim with engine deflections only 

npc()5) -0, /do not calculate acrohceting 

npc/26) -0, /no special aeroheating calculations 

npc(38) - 0, / no atmospheric turbulence 

npc(39) - 0, / no atmospheric gusts 

npc(40) -1,/ 


c Propulsive Data 


npc(9) - 1, / rocket engine with thrust cable and i*>vac 
iengt - 1, / calculate engine incidence angles from static trim equations 
npc(l3)-0,/ do not jettison propellant 
npc<22) - 2, / aat ett to value inputted by etape 
etapc(l) - 1.0, /throttle value 
npc(27) - 0,/do not integrate engine flowrMw^ 

/the throttling is kind of complicated 

/the engine can be throttled to meet dynamic preeure and wing loading 
/constraints during ascent, but the solid motori cannot be throttled 
/once the vehicle gets to 3g's the post is instructed to limit the acc to 3 
ienga(l)- -1,0. /throttle orbiter using input eta, but not strap on* 
iengmflt)- 1,1, /orbiter and strip rockets are on 


I I A— ■ 



wprp(l) - 
men stp - 


1, 

2. /highest number of engines 

1 ,2, /map each engine to a specific tank 
7 529Se5, /weight of sit propellent 

\X 


c 


c Weight model 


npc( 1 7) - 2, / jettison booster 
npc(30) - 3. / use n-suge vehicle weight model 
wstpd/l) - 7529 5e5. /dry weight of srb 

wstpd(2) - 121*2*3 ,/groasweight of orbiter 


c Random Stuff 


npc(14) - 0, / no holddown before liftoff 
npc( 18) - 0, /do not terminate trajectory 

npc(2 1 ) - 0, / do not calculate fuel and oxidizer weights and volumes 

npc(23) - 0, / do not compute velocity margins 

npc(24) - 0, / do not compute any parameter integrals 

npc(2S) - 0, / tracking station option not used 

npc(29) - 0, / do not compute vacuum impact points 

npc(3l) - 0, / no vernal equinox, sun shadow, sun angle calculations 

npc(32) - 0, / no parachute drag 

npc<33) - 0, / no arc length calculations 

npc(34) - 0, / no kepicrian state calculations 

npc(35) - 0, / no sensed velocity increment calculation* 

npc(36) - 0, / do not activate sunlight options 

npc(37) - 0, / do not activate date option 


c Tables 


* 

pStblmlt 
edm — II. 

tvc2m - 3 3351, / thrust multiplier for orbiter 
wd2m - 3 3351 ,/ flow rate multiplier for orbiter 
ae2m - 3 335 1, / exit area multiplier for arbiter 
$ 

pStab 

table - Shtvc2t, 0,531 500., / vacuum thrust of orbiter engine 
$ 

p$tab 

table - 5htvc1t, 0.2 2406e6. / vacuum thrust of strap-oru engine 
S 

pStab 

table * 4hwd2t,0,1 199 77, / flow rate of orbiter engine 
S 

pStab 

table - 4hwdlt,0,BI40. / flow rate of strap-ons engine 
I 

pStab 

table - 4hae2t, 0.33 551, /exit area of orbiter engine 
S 

pStab 

table - 4hael t, 0,1 13 73. / cut area of stnjMmi engine 
S 

pStab taMe - 4hxcgt, 1 ,6hweight,2, 1 .* 1 . 1 . 

2006450,123 31.134*70,1 15 50. 

S 
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pStab tab»e-5hxreft,0,109 90, 

S 

pStab table-5hzreft,0,0. 

$ 

•indude 'Awmc/aadl l/)tooley/poat/BIMESE/bimeae aero' 
endphs - 1, 


$ 

ISgendat 

event-M,l,critT- Shtime , 

value *20, /change throttle to acceleration limiting 

xmax(S) - 0, /react throttle ao can make aure etape/ I ) isn't aet to above one 

«*apc(l)- 5, 

endphs - 1, 

$ 

ISgendat 

event- 1 5,critr-4htii»e,vahie-07 .endphs- 1 , 
igui<MH>, 

S 

ISgendat 

event- 1 6,cri tr-4htime, vahie-20, endphs- 1 , 

S 

ISgendat 

event- 1 1, 1 ,cri tr-4hasmg, 

value - 2.999, /change throttle to acceleration limiting 

i*ui<M)-2, 

npo(7) - I, / limit acceleration 
asmax - 3.0, / limit acceleration to 3 g't 
npc(22) - 2, 

ienga(l)- 0,1, /throttle arbiter tiling acc limit, but not atrip on i 
etapc(l)“ I, 
endphs - I, 

S 

ISgendat 

event-1 9, 1 ,chtr-5hwprp I ,value-0, 

value -0, /turnon new engine when 1st stage lud is empty 
natpl -2, 

icngmfll) - 0, /strap ons off 
npc(7) - 1, / limit acceleration 
asmax -30,/ limit acceleration to 3 g's 
npc(22) - 2, 

icnga( I ) - 0, l , / throttle arbiter using acc limit 
etapc<1)- I. 

S 

pStblmlt 
edm - 1.0, 


pStab 

table - 5htvc2t, 0,53 1 500 , / vacuum thrust of orbiter engine 
endphs - 1 , 

ISgendat 

event-20,critr-4htiine, value-5 1 .endphs- 1 , 

ISgendat 

event-30,critr-4htime, value-80, endphs- 1 , 

ISgendat 

event-40,chtr-4htinic,value«1 1 0.endphs- 1 , 

ISgendat 

event-60,critr-4htime, value- 1 40, endphs- 1 , 

ISgendat 

event-65 ,critr-4ht»me,value-l 70, endphs- 1 , 

ISgendat 

even t-70,cntr-4hti me, vaJue-200. endphs- 1 , 

ISgendat 

event-75 > chtr-4htime,vaJue-230 l endphs- 1 , 

ISgendat 

event-30,critr»4htimc, value-240, endphs- 1 , 

ISgendat 

cfvent-S 5,critr-4htime,vaiue-250,endphs- 1 , 

ISgendat 

event-90, cntrMhti me, value- 260, endphs- 1 , 

ISgendat 

c this event marls injection into orbit 
event- S00.cn tr-5hwprop,value-~99973 3, 
endphs - 1 , 
endjob - 1 , 
endprb - I, 

S 
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A.4. MATED BIMESE: LEO 


PSSEARCH 

€•••*••• ••••••••*• 

C LOX/LH2 BIMESE W/RS2100 Engine (reviaed 1/1/99) 

C WB-003-C Bimeae 60 Klb to 100 nmi 
C 21 5 DEGREES INCLINATION 
C LATEST WB-001 AERO ADDED 5/6/97 

* 

c 

SRCHM - 4, / PROJECTED GRADIENT METHOD 
C IPRO - -1. / PRINT FINAL TRAJECTORY ONLY 
IOFLAG - 0, / ENGLISH IN - ENGLISH OUT 
CONEPSO) - *9 9, / CONVERGENCE TOLERANCE 
CONEPS<2) - 0000 1 . / MIN ALLOWED PERCENTAGE CHANGE 
CONEPSO) - 00001,/ 

CONEPS(4) - 00001, / 

CONEPSO) - 00001,/ 

PCTCC-001, / 

IDEB-0. / No Detail 

C 



c OPTIMIZATION VARIABLE 

** 

c 

OPTVAR - 6HWEIGHT, / MAXIMIZE FINAL BURNOUT WEIGHT 
OPT - 1, / MAXIMIZATION FLAG 

OPTPH - 500, / PHASE AT WHICH WEIGHT IS MAXIMIZED 
WOPT - 000003, / OPTIMIZATION VARIABLE WEIGHTING 
MAXITR - *1, / MAXIMUM ITERATIONS ALLOWED 
C 



c 

C CONSTRAINT VARIABLES 

C 

***** 

C 

NDEPV - 6, / NUMBER OF DEPENDENT VARIABLES USED AS 
CONTROLS 

DEPVR - 5HGDALT. 6HGAMMA1, 3HINC, 5HXMAXI, 5HXMAX2. 
5HXMIN2, / INSERTION ORBIT CONDITIONS «■ CONSTRAINTS 
DEPVAL - 303805 0, 0 0. 2* 5, 1000 , 379165 . -379165 , / INSERTION 
ORBIT VALUES + CONSTRAINTS 
DEPTL - 300 0, I, 1,1., 100. 100, / DESIRED ACCURACIES 
DEPPH - 500. 500. 500. 500, 500. 500 . 1 EVENTS 
©EPVR - 0, 0. 0. I. l.-l. 

C 

* * 

C 

C CONTROL VARIABLES 

C 

C MM 

c 

mndv - 16, 

indvr - 1 5*6hpitpc2,3hazJ, 

indph - 15, 18, 20.30.40.60, 65, 70. 75,80.8 5. 90,9 5, 100. 105, 10. 

u- *9.631 1466 1 293 2E-01, I 2068533 784 24E-0 1, 3 371722786330E-02. 

-8 3 3 5495035323E-OI, -4 798I49007642E-01, 

-4 000001681097E-01. -4 3480200 5346 5E-02. -I 4 1 541951 5I40E-0I. 

-4 0362450 1 7 758E-0 1.-991 5998643 7 1 2E-02, 

- 1 85 1 2988967 1 7E-03 , -3 262486724374E-02, -2 64803 11553 82E-02. 

•4 657733088827E-01, -1 3 1291730I909E-01, 

00, 

u- *1 1121 105 16660EXX), 9 59606184 1384 E -02, 2 4057203 1201 2E-02, * 

8 6478607 1 S809E-0 1,-8 5 2 446898 7230E-0 1 . 

-7 365 1 35456373E-01, -6 1 85493375 131 E-02, -1 78091 5656336E-01 , - 
2 827 569253803 E-0 1 . *2 61715323 1 269E-01. 

-1 977722 599479E-03, -7.98I303256862E-02, -3 94965587895 1E-02, - 
5 9017533 14859E-01, -1 333034849750E-0I. 

8.621 963895029E-KM. 

PERT - 16*1 0E-4, /PERTURB 

8 

* * 

C 

C TRAJECTORY SIMULATION INPUTS 

C 

• ********" 

PJGENDAT 

TITLE - OH* LOX/LH2 BIMESE 2* 5 deg delivering 60 Klb *, 

EVENT - 10, / FIRST EVENT 

WGTSG - 2274254 , / INITAL LIFTOFF WEIGHT with ADDED FLUIDS 

FESN - 500, / FINAL EVENT NUMBER 

MAXTIM - 1000 , / MAXIMUM TIME 

ALTMAX - 10000000 . / MAXIMUM ALTITUDE ALLOWED 

ALTMIN - -5000 . / MINIMUM ALTITUDE ALLOWED 

PRNC - 0, /setup plot file 

MONX - 4HDYNP, 4HFAZB. / MONITOR Q ft NORMAL FORCE 
PRNT(97) - 5HXMAX1, 5HXMAX2, 5HXMIN2. 6HALTITO, 4HMACH, 
4HDYNP.6HNET1SP, 

4HAMYB, 5HTTMYB, 5 HP STOP, / 

1GU1D( 1 ) - 1, / INERTIAL EULER ANGLE OPTION 


)GUID(2) - 0. / SAME STEERING OPTION FOR ALL CHANNELS 
IGUID(4) - 1, / CUBIC POLY WITH CONSTANT TERM SET BY INPUT 
PITPC( 1 ) - 0 . / INITIAL PITCH RATE 
YAWPC{ 1 ) - 0 , / INITIAL YAW RATE 
ROLPC<1 ) - 0 , / INITIAL ROLL RATE 

NPC( l ) - 2, / PRINT CONIC BLOCK AT EACH PHASE CHANGE 
PINC - 20, / PRINT INTERVAL EVERY 20 SECONDS 
NPC(2) - I , / FOURTH ORDER RUNGE KUTTA 
DT - I. / INTEGRATION STEP SIZE 

NPQ3) - 4, / PLANET RELATIVE INPUT ON VELOCITY VECTOR 
VELR - 1. / INITIAL VELOCITY AT LAUNCH PAD 
GAMMAR - 90 , / INITIAL PLANET RELATIVE FLIGHT PATH ANGLE 
AZVELR - 90 , / INITIAL AZIMUTH ANGLE OF VELR 
NPQ4) - 2. / PLANET RELATIVE INPUT ON ALTITUDE VECTOR 
ALTTTO - 0 . / INITIAL ALTITUDE 
GDALT - 0 , / INITIAL ALTITUDE 
GDLAT - 28 5, / LATTTUDE OF LAUNCH SITE 
LONG - 280 0, / LONGITUDE OF LAUNCH SITE 
LONGI - 280 0, / LONGITUDE EAST OF PRIME MERIDIAN 
NPC(5) - 5, / 1976 US STANDARD ATMOSPHERE 
NPC(6) - 0, / NO WINDS 
NPC(7) - 1, / ACCELERATION LIMIT SET 
ASMAX -30./ LIMIT ACCELERATION THROUGHOUT TRAJECTORY 
LREF - 162 6, /BODY LENGTH FROM NOSE TO BASE 
SREF - 7001 4, /REFERENCE AREA 
NPC<8) - 2, / CL, CD ft CM TABLE INPUTS 

NPC(9) - I, / ROCKET ENGINE WITH THRUST TABLE AND ISPVAC 
NPC(IO) * 0, / STATIC TRIM IN PITCH 

GXP( I) - 160 2. /X-LOCATION OF ENGINE GtMBAL MEASURED FROM 
NOSE 

ITRIM - I, / TRIM WITH ENGINE DEFLECTIONS ONLY 
IENGT - I, / CALCULATE ENGINE INCIDENCE ANGLES FROM STATIC 
TRIM EQUATIONS 

NPC( 1 1 ) - 0, / NO FUNCTIONAL INEQUALITY CONSTRAINTS 
NPC( 12) - 3, / CALCULATE CROSSRANGE AND DOWNRANGE BASED 
ON ORBIT 

ALTREF - 100 0, / 100 NMI CIRCULAR MILE REFERENCE ORBIT 
AZREF - 0 , / AZIMUTH REFERENCE 
NPC(13) - 0, / DO NOT JETTISON PROPELLANT 
NPC( 14) - 0, / NO HOLDDOWN BEFORE LIFTOFF 
NPC( 1 5) - 0, / DO NOT CALCULATE AEROHEATTNG 
NPC( 16) - 0, / USE OBLATE EARTH GRAVITY MODEL 
NPC( 1 7) - 2, / JETTISON BOOSTER 
NPC( 1 8) - 0, / DO NOT TERMINATE TRAJECTORY 
C NPC(I9) - 0, / DO NOT PRINT INPUT CONDITION SUMMARIES 
NPC(20) - 0, / DO NOT USE ANY SPECIAL DT CALCULATION 
NPC<2 1 ) - 0. / DO NOT CALCULATE FUEL AND OXIDIZER WEIGHTS 
AND VOLUMES 

NPC(22) - 0, / DO NOT CALCULATE THROTTLING 
NPCI23) - 0, / DO NOT COMPUTE VELOCITY MARGINS 
NPC(24) - 0, / DO NOT COMPUTE ANY PARAMETER INTEGRALS 
NPC(25) - 3, / COMPUTE VELOCITY LOSSES AND PRINT AT EACH 
TIME 

NPC(26) - 0, / NO SPECIAL AEROHEATING CALCULATIONS 
NPCX27) - 0, / DO NOT INTEGRATE ENGINE FLOWRATES 
IWPF( I ) - 0, / DO NOT INTEGRATE MAIN ENGINES 
IWPF{2) - 1, / INTEGRATE BOOSTER FLOW RATES 
NPC(28) - 0, / TRACKING STATION OPTION NOT USED 
NPC(29) - 0, / DO NOT COMPUTE VACUUM IMPACT POINTS 
NPC(30) - 0, / USE N-STAGE VEHICLE WEIGHT MODEL 
IENGA( 1 ) - 1 , / CALCULATE BOOSTER THROTTLE NECESSARY TO 
LIMIT ACCELERATION 
IENGA(2) - 1. / THROTTLE ORB I TER 
IENGMF( l ) - I , / BOOSTER ENGINES ARE ON 
IENGMF(2) - 1 , / ORB ITER ROCKETS ARE ON 
IWDF(I) - I. / CALCULATE BOOSTER FLOW RATES AS A TABLE 
LOOKUP 

IWDF<2) - 1, / CALCULATE ORB ITER FLOW RATES AS A TABLE 
LOOKUP 

NENG - 2. / THE NUMBER OF THRUSTING ENGINES 
NPC<3 1) - 0, / NO VERNAL EQUINOX SUN SHADOW. SUN ANGLE 
CALCULATIONS 

NPQ32) - 0, / NO PARACHUTE DRAG 

NPC(33) - 0, / NO ARC LENGTH CALCULATIONS 

NPC(34) - 0, / NO KEPLERIAN STATE CALCULATIONS 

NPC(35) - 0, / NO SENSED VELOCITY INCREMENT CALCULATIONS 

NPQ36) - 0, / DO NOT ACTIVATE SUNLIGHT OPTIONS 

NPQ37) - 0. / DO NOT ACTIVATE DATE OPTION 

NPQ38) - 0. / NO ATMOSPHERIC TURBULENCE 

NPC(39) - 0, / NO ATMOSPHERIC GUSTS 

NPC<40)« I./ 

S 

PSTBLMLT 

TVCIM - 3.3351, / THRUST MULTIPLIER FOR BOOSTER 
WD1M - 3 3351, /FLOW RATE MULTIPLIER FOR BOOSTER 
AEIM - 3 3351, /EXIT AREA MULTIPLIER FOR BOOSTER 
TVC2M - 3 3351,/ Thrust multiplier for Ortwter 
WD2M - 3 3351, / Flow rate multiplier for OrtHtcr 
AE2M - 3 3351, / Exit area mulbplier for Orb iter 
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s 

PSTAB 

TABLE - 5HTTVCIT,0,531500 , t VACUUM THRUST 

S 

PSTAB 

TABLE - 5HTVC2T, 0,531 500 , / VACUUM THRUST OF Orbiter engine 

S 

PSTAB 

TABLE - 4HWDIT.0.1 199 77. / FLOW RATE 

S 

PSTAB 

TABLE - 4HWD2T.0J 199 77, / FLOW RATE OF Orb iter engine 
S 

PSTAB 

TABLE - 4HAEIT,0,33 558, / EXIT AREA 

S 

PSTAB 

TABLE - 4HAE2T, 0,3 3 5 54, / EXIT AREA OF Orbiter engine 
S 

PSTAB TABLE - 4HXCGT, 1 .6HWE1GHT.2, 1,-1, l. 

2006450, 1 23 .5 1 , 1 34870, 11 5 50, 

S 

PSTAB TABLE-5HXREFT.0.109 90. 

S 

PSTAB TABLE-SHZREFT.0,0, 

S 

•include Thome/asdl l/jtooiey/poet/BIMESEAnmeae aero’ 

ENDPHS - 1, 

$ 

ISgendat 

event* 1 5,critr-4htime,vnlueK)7.endphs- 1 , 
iguid<4)-0, 

$ 

ISgendat 

event- 1 8,critr-4htime, value- 1 7,endphi- 1 . 

$ 

ISgendat 

evcnt-20,critr-4htime,value-30,endphi- 1 . 

S 

ISgendat 

event-30,critr-4htime,value-40,endplit- 1 . 

S 

ISgendat 

event-40,critr-4litime,vahie^0,endpli»- 1 . 

S 

ISgendat 

event^0,critr-4htinie,value-80,endph»- 1 , 

S 

ISgendat 

event-bS.critr-dlitime, value- lOO.endphs- 1 , 

$ 

ISGENDAT 

EVENT - 35, CRJTR - 6HWE1CON, VALUE -909679 . 

WJETT -174031,/ JETTISON DRY Booeter 
1ENGMF( 1 ) - 0. / TURN Booeter OFF 
JENGMF(2) - 1 . / TURN Orbiter ON 
IWDF( I ) - 0, / FLOW RATES OFF 
IWDF<2) - 1 , / FLOW RATES ON 
IENGA( 1 ) « 0, / 

IENGA(2) - 1,/ 

IWDF( 1 ) - 0, / 

IWDF(2) - I./ 

LREF - 162 6, /Orbiter 
SREF - 3500 7, /Orbiter alone 
ENDPHS - 1, 

S 

ISgendat 

event- 70,cn tr-4ho me, value- 1 45,endph>- 1 , 

S 

ISgendat 

event-75.critr-4hbme, value- 1 85,endpb*- 1 , 

$ 

ISgendat 

event-80,critr-4htime, value- 220,endplii- 1 . 

S 

ISgendat 

event-85.cntr-4htime,value-255,endphi- 1 , 

S 

ISgendat 

event-90,critr«4htime, value-290, endpht- 1 , 

S 

ISgendat 

event-95, cntr-4htime,vaJue-330,endph»- I, 

S 

ISgendat 

event- 1 00,critr-4hti me, value- 360,endphs- 1 , 

S 

ISgendat 

event- 1 05,cntr-4htime,vahie-390,endphi- 1 , 

S 

ISgendat 

c this event marks injection into orbit 
event-500,critr-4hveli,value-25855 000, 


ENDPHS - 1, 
ENDJOB - I, 
ENDPRB - I, 

S 


Georgia Tech SSDL 


Year 3 Final Report, Improving Conceptual Design for Launch Vehicles 



NCC l«229 


The Bunese Concept: A Study of Mission and Economic Options 


52 


A.5. SINGLE ELEMENT: POINT-TO-POINT 


PSSEARCH 

■ 

C LOX/LH2 BIMESE W/RS2 100 Engine (revised 1/1/99) 

C WB-003-C Bimeee point to point 
C *agle Orbuer cannot make it to orbit 
C Goal it to maximise crossrenge of arbiter 
C Simulates the ascent (but not to orbit) and reentry 

C The Gliding entry phase is simulated down to 10,000 ft and a rd vel of 2 500 
life 

C At tbit point energy management phase would take over, but it n not 
■mutated 

C Energy management phase usually remit* in a 100 nmi croatrange 
C LATEST WB-001 AERO ADDED 5/6/97 


c 

SRCHM - 4, / PROJECTED GRADIENT METHOD 
c IPRO - *1 . / PRINT FINAL TRAJECTORY ONLY 
IOFLAG - 0. / ENGLISH IN - ENGLISH OUT 
CONEPSd ) - 19 9, / CONVERGENCE TOLERANCE 
CONEPS/2) “ 0000 1 . / MIN ALLOWED PERCENTAGE CHANGE 
CONEPS<3)- 00001./ 

CONEPS<4)« 0000 1 . / 

CON EPS/ 5) - 00001./ 

PCTCC-001,/ 

IDEB-0. / No Detail 
C 

* 

c OPTIMIZATION VARIABLE 

c* MM * * 

c 

OPTVAR -6HDPRNGI. /MAXIMIZE DISTANCE TRAVELED 
OPT - 1, / MAXIMIZATION FLAG 

OPTPH - 500. / PHASE AT WHICH WEIGHT IS MAXIMIZED 
WOPT - 000003, / OPTIMIZATION VARIABLE WEIGHTING 
MAXTTR - -1, / MAXIMUM ITERATIONS ALLOWED 
C 

c* ••• 

C 

C CONSTRAINT VARIABLES 

C 



c 

NDEPV - 5, / NUMBER OF DEPENDENT VARIABLES USED AS 
CONTROLS 

DEPVR “ 5HXMAX3.5HXM AX 1 . 5HXMAX2. 5HXMIN2. 5HXMAX4. 
DEPVAL - 4. 1000 , 379165 , *379165 .40, t INSERTION ORBIT VALUES + 
CONSTRAINTS 

DEPTL- I. 1. 100, 100. I, /DESIRED ACCURACIES 

DEPPH - 500,500,500.500,500/ EVENTS 

IDEPVR - I, 1, 1,-1. I. 

* * 

c 

c 

c CONTROL VARIABLES 

C 

* 

C 

tabl(l) “4bpin,4bpitt,4hpitt.4hpitl,4HpttL4hp»tl.4hpttt,4hpin, 
tabl(9) 

-6haIphat,6halphat,6haJphai,6h*lphat,6haJph*t,6hilphat,6halphat,6halphat, 

tably( 1 )r 1 ,2,3. 4, 5, 6,7, 1, 
tably<9)-3.4,5.6, 7.8.9.10. 

mndv - 16, /In this deck banking angles are not needed to meet constraint* 
indvr - 5htabH,5htabl2,5htabl3,5htabl4.5htabl5,5btabi6,5Ktabl7 1 5htabl8, 
indvr/9) 

-5htabi9,6htab! 1 0,6htaW 1 1 ,6htaW 1 2,6htabl 1 3 ,6htabl 1 4,6htabl 1 5.6htabl 1 6, 
indph - 15,15,15.15.15.15,15,15. 
indph(9) - 20.20,20,20,20,20,20,20, 

PERT - 16* l OE-4,1, / PERTURB 

4 2741 1 1864380E-03.-2 1 33260052939E+0 1 ,*2 39766697477 5E+0 1 , 

-5 6265 5 58263 70E+O l ,-7 36 1 997959382E+0 1 .4 6222274 1 886 1 E+0 1 , 

-9 33681695769SE+01,-9 8 7396151B200E+01, 1 627310476026E+01, 

u- 4 274245 165 569E-03.-2 1391 10649241 E+0 1.-2 42082430495 1 E+0 1. 

-5 6I8058820848E+01, -7 3465160141 14E+01, 4 1 7178251 3255E+OI, 

-9 367542853691 E+0 1 .-9.860263846201 E+0 1 , 1 627300997425E+01 , 

/u'» that result in opt with only alpha control violates normal force 

*1^4 2742 1 74 1 3 973 E-03 .-2 1 2 1 879899720E+0 1 ,-2 4 1 3003298457E+0 1 , 

-5 462767091 147E+01.4 475214444848E+OI.4 I52759020379E+01, 

-9 2792 1 8168809E+0 1 ,-9 86026384620 1 E+0 1 , 8 65 1 244888 1 45E+0 I , 

40.5 951 168900506E+01, 4 140549 13 794 1 E+0 1, I 06546U4I837E+OI. 

5 929888736206E+00, 5.058832037600E+00. I 32763 1 52423 1 E-02,- 
60.60.60.40,40,20. 

u- 4 2742 1 74 1 3973E-03.-2 1 2 1 879899720E+OI .-2 4 1 3003 2984 57E+0 1 , 

-5 462767091 147E+01.4 4752 14444848E+0 1,4 1 527590203 79E+0 1, 
-9.2792 ) 8 1 68809E+0 1 ,-9 860263846201 E+0 1, 8 65 12444881 4 5E+01, 

40.5 951 168900506E+01. 4 1 40549 13 7941 E+0 1. 1.065461 14I837E+01, 


5 9298B8736206E+00, 5 058832037600E+00, 1 .327631 52423 IE-02,- 10,- 10,- 
10.-10.-10,-10, 

/converged solution without constraining final velocity 

u» 4 2742 1 74 1 3973E-03.-2. 1 2 1 879899720E+0 1 ,-2 4I3003298457E+OI. 

-5 46276709 1 1 4 7E+0 1 .4 4752 1 4444848E+0 1 ,4 1 527590203 79E+0 1 , 

-9 2792 1 8 1 68809E+01 , -9 86026384620 1 E+0 1 , 7 1 7638 5 530439E+0 1 . 

4 00000850281 5E+OI, 6. 10I418902565E+01. 4 .258546163 1 30E+O1, 

1 0636S5205035E+01, 5.913 129650736E+O0, 5.0452321 93 202E+00, 

1 32763 1425463E-02, 0 9997745204 59E+00, 0 98477375 11 25E+O0. 
0976526577I65E+00, 09951 85 582387E+00, 0 0000058 761 46E+01, 
OOOOOOOOOOOOOE+Ol. 

/converged with velocity constraini 

u- 4 .2742 1 866413 8E-03.-2. 1 2234854973 7E+0 1,-24 1 5495329997E+01 , 

-5 532985 1 69107E+01 ,4 484468 1 8 1437E+0 1 ,4 1 490728794 58E+0 1 . 
-9.22508 7779792E +01 .-9 849M9546882E+0 1 , 

4 072 1 971 86392E+01 , 4.04649924860 1 E+00, 4 0982706*09 1 E+00, 

4.027990 176507E+0 1, 

3 027990176507E+OI.2327990I765O7E+OU 327990I76507E+OU 327990176 
507E+01, 

/use u's for 1,000 lb payload 

u- 4 2742790397 59E-03.-2 I420H663S03E+01.-2 36I182370361E+OI, 

-5 6I8851516036E+0I.-9 697944809427E+01 .-6 83925B291462E+01, 

-9.3660371 761 18E+01.-9 99042473 8464E+0 1 . 4 00001 33 1 1904E+01 . 

4 000000031 147E+01, 3 091 1 9389643 7E+01, 3 O47280667995E+0I. 

2 0*74 1 8722764E+0 1 , 2 3208440 1 1004E+0 1,1315381 826792E+0 1 , 

0 317517262700E+01, 

/use u’s for 1,000 lb payload 

u- 4.27428 5458 1 1 6E-03.-2 144825875449E+0 1 , -2 3590790 1 1 73 5E+0 1 , 

-5 5 1 96497 1 8889E+01 .-9.249077673 248E+0 1,-7 11 69562 74 567E+0 1 , 

-9 258437443429E+0 1 ,-9 99(M24738464E+Ol , 3 9 1 0694 1 25902E+0 1 . 

3 9275608664 27E +01 , 3 0481 726042 14E+0 1, 3 01 581 000080 7E +01, 

2 035214171299E+01, 2 30969461 7097E+01, 1 3I7472828524E+01. 

3 1 76752890923E+00, 

u- 4 274285001 159E-03, -2 1 337409 74 347E +0 1 , -2 37238262048 1 E+0 1, - 
5.502854357052E+01, -9 204608000024E+01, 

-7 1 274493 78298E+0 1. -9 .26017192 836 7E+C1 . -9 9904 24 73 8464 E+0 1, 

3 910694125902E+01, 3 927560866427E+01, 

3 048 1726042 14E+01, 3 OI58I0000807E+01, 2 035214171299E+01. 

2 0096946 1 7097E+0 1 , I 3 17472828524E+01 , 

3 17675 289092 3 E+00, 

u- 4 274284650071 E-03, -2 42752 1889897E+01, -2 5744043603 1 6E +01, - 
5 500979 194040E+01, -9 2024261 47050E+01, 

•7 1 28267284684E+0 1 , -9 2602634 74492E+0 1, -9 990424738464E+01, 

3 910694125902E+01, 3 927560866427E+01 , 

3 048I726O4214E+OI, 3 OI5810000807E+01, 2 0352141 71299E+01 , 

2 0096946 17097E+01. 1 317472828524E+01, 

3. 1 76752890923E+00, 

u- 4 2742774053 5 2E-03, -2 2034655 5 2405E +01, -2 646959I48998E+0I. - 
5 486751 3591 44E+01, *9 155510109285E+01. 

-7 1 307356743 52E+01, -9 21 1244742081 E+0 1, -9 990424738464E+0I, 

3 9046982976 19E +01, 3 914151458512E+01, 

3 039242625757E+01, 3 009 199393551 E+0 1, 2 032308900060E+0!, 

2 007854985308E+01, I 3 1 7695 197695E+0 1 , 

3 1769307894 1 6E+00, 

u- 4 . 274292444634 E-03 , -2 2I3526273097E+01, -2 655677633 584E +01, - 
5 520896242270E+0 1 . 4 704952345146E+01, 

-7 487853389 1 45E+01, 4 69452 1 8373 72E+01, -9 990424 7 3 8464 E+0 1. 

3 794469689704E+0 1 , 3 618706367882E+01, 

2.926431 1 3601 5E+0 1. 2 942776111 155E+0 1, 1 992022902 4 36E+01 . 

1 99 1 2682 1 0402E+0 ! , 1.321 554806599E+0 1 . 

3. 1 80024775060E+00, 

u- 4 2743026O9566E-O3, -2 19810973 1 806E +01, -2 6874523 57529E+01, - 
5 54532 1 538596E+0 1 , 4 475943 1 1 2996E+0 1 , 

-7 36522691 3249E+OI. 4 5861901 3842 1 E+01 , -9 990424738464E+01 . 

3 683948762582E+01, 3 333347120303E+01, 

2 847403 324806E+0 1 , 2 890967570875E+0I, I 973721 71 2864E+01. 

I 98877925 1295E+01, I 32539944671 1 E+01, 

3 18306249909 1 E+OO, 

u- 4 2743005309 1 6E-0 3 , -2 193239408865E+01. -2 6669284 3 383 9E +0 1, - 
5 97240943776SE+01. 4 5705 347505 54E+01, 

4.2151 47334882E+0 1 , -7 8891 39786 174E+01, -9 26746530 183 2E +01, 

3 4901 642343 76E+0 1, 3 2229 10009984 E+0 1 . 

2 82 1 906569241 E+0). 2 •7I116523180E+0I, 1 967036S25028E+O1, 

1 9878341 37037E+O l, 1 3301 955345 12E+01. 

3 1864894 13846E+00, 

u- 4.274286253070E-03, -2 1 86 1 6067983 2E +01, -2 669959867244E+0 1 . - 
5 722351 248990E+O1. 4 73082557923 1E+01. 

-7 5643524O0382E+O1. 4 335736962521 E+Ol, -9 693 56 54 30695 E+01. 

3 49101 0638 185E+01, 3 223600706208E+01. 

2 82205294 5986E+01. 2 87 123771 43 34E+01, 1 967077 5 8043 7E+0 1. 

1 987840475437E+01, I 330172986096E+01, 

3 . 1 8647308 5384E+00, 

u- 4 2743005309 16E-03, -2 l9323940t865E+01, -2 6669284 33 839E +01, - 
5 97240943 7765E+0 1 , 4 5705347505 54E+01. 

4 2151 47334882E+01 . -7 8891 39786 1 74E+01 . -9 26746530 1 832E+0 1 , 

3 490 1 64234 3 76E +01, 3 222910009984E+OI, 

2 82 1 906569241 E+Ol. 2 871116523 180E+01, I 9670368 2 5028 E+01, 

1 987834I37037E+01, 1 330I955345I2E+01. 
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3 18648941 3846E+00, 

u- 4 27430080763SE-03, -2 193378482791 E+01, -2 66*8844 2078 7E+01, - 
5 980620048027E+0 1 , 4 58464648 3 286E +01, 

-8.23918874351 IE+01, -7 8906181 37000E+01, -9 25861 771I327E+01. 

3 490I64234376E+O1, 3 2229 10009984 E +01, 

2 82 1 906569241 E+01, 2 871 1 1 6523 1 80E+0 1 , 1 96703682SO28E+O1, 
1.9878341 3703 7E+01, 1.33019553451 2E+0 1, 

3 1864894 13846E+O0, 

u- -8.27430142791 4E-03, -2 19367767 1608E +01, -2 665854637207E+01, - 

5 9546959948 1 8E+01 , -8 5 1 8 1 1 578 1 1 98E+0 1 , 

4 22893975 1839E+01, -7 89871938 13 59E+01, -9 261 47626 1244E+01, 

3 .486356771 167E+01, 3 2I9810154606E+01, 

2 I21252839576E+01, 2 8705 74049 722E+0 1, 1 966858499644E+01. 

1 9878058 5 2482E+01 , 1 330297O94262E+O1. 

3 1 86562801 2 10E+00, 

U- 4 274329468 128E-03, -2.207572164824E+01, -2 648082538327E+01, - 

6 1 240298 54938E+01, -7 949261935 1 83 E+01, 

4 65408689 1 64 7E+01, -7 85740145351 5E+01. 4 81830731 9575E+01, 
3.469786 1 33216E+01, 3 2064 10790 120E +01, 

2.818501554105E+01, 2 868246733 162E+01. 1 96608404072 tE+Ol, 

1 9876183 74684E+01 , 1 330730628828E+OI, 

3 18687623 7782E+00, 

w- 4 27433308008 IE-03. -2 .2088933251 19E+01, -2 643193837210E+01, - 
6.097657 1 1 9546E+0 1 , -7 .86807 1 8 7206 1 E+0 1 . 

4 660725274903E+01, -7 992194878060E+01, 4 83801 1470419E+01, 

3 .427335034 189E+01, 3 1724408 3 7729E+01, 

2 8I1822072577E+OI, 2 8623 80238 5 84E+01, 1 964171632818E+01. 

1 987402867 1 45E+01, 1 331822394620E+01, 

3 18766405 1 246E+00, 

u- 4 274329468 1 28E-03, -2.207572 1 64824 E+0 1, -2 648082538327E+0I. - 
6 124029854938E+01, -7 94926 1935 183 E+0 1, 

4 65408689 1 647E+01 , -7 85740 1 453 5 1 5E+0 1 , 4 8 1 83073 1 9575E+0 1 , 
3.4697861 332 16E+01, 3 2064 1 0790 1 20E+0 1 . 

2.8 18501 554 105E+01, 2 868246733 162E+01, 1 966084040721 E+01, 

1 987688374684E+01, I 3 30 7 3 06288 28E +01, 

3 18687623 7782E+00, 

/63s chump 

u- 4 27433308008 IE-03, -2.2088933251 19E+0 1, -2 643193837210E+01. - 
6.097657 1 1 9546E+01, -7 86807 1 87206 1 E+0 1 . 

4 6607252 74903 E+01, -7 992194878060E+O1, 4 83801 1470419E+01, 
3.42733 5034 189E+01, 3 172440837729E+01, 

2 811822072577E+01, 2 8623 80238 584E+01, 1 9641 7 163281 8E+01, 

I 987402867 145E+01, 1 331822394620E+OI, 

3 18766405 1246E+00, 

u- 4.27433361 3270E-03, -2 20258549712 1 E+O I. -2 645259999024E+O1.- 
6.0633 76525407E+01, *7 8048444608 16E +01, 

4.62359873007 1 E+0 1 , 4.0071 983 285 77E +01 , 4 83801 1470419E+01, 
3.41958 1771260E +01, 3 1642 12737 1 07E +01, 

2 810181328856E+OI, 2 860985929624E+01, l 963703285971 E+O l . 
t 987333933870E+01, 1 3 3 20499292 77E +01, 

3 1 8782204493 7E+00, 

u- 4.274333080081 E-03, -2 2088933251 1 9E+0 1, -2 6431938372 1 0E+0 1, - 
6 097657 119546E+0 1 , -7 86807 1 87206 1 E+0 1 , 

4 660725274903E+01, -7 992 1 94878060E+0 1 , 4 83801 1470419E+0I. 

3 427335034 189E+01, 3 172440837729E+01, 

2.81 1822072577E+OI, 2 862380238584E+01 . I 964171632818E+01, 

I 987402867 145E+0 1, 1 33I822394620E+01, 

3 18766405 1246E+00, 

$ 

C* * 

c 

C TRAJECTORY SIMULATION INPUTS 

C 

* * •••• 

PSGENDAT 

TITLE - OH* LOX/LH2 BIMESE 28 5 deg delivering 60 Klb *, 

EVENT - 10, / FIRST EVENT 
WGTSG - 1122050, 

FESN - 500. / FINAL EVENT NUMBER 
MAXTIM - 3000, / MAXIMUM TIME 
ALTMAX - 10000000 , / MAXIMUM ALTITUDE ALLOWED 
ALTMIN - -5000 , / MINIMUM ALTITUDE ALLOWED 
PRNC - 0, /setup plot file 
c PRNCA - 0, /setup delimited plot file 
MONX - 4HDYNP, 4HFAZB, 4HASMG, 5HALPHA, 6HGAMMAI, / 
MONITOR Q * NORMAL FORCE 

PRNT(97) - 5HXMAX1, 5HXMAX2, 5HXMIN2, 6HALTITO, 4HMACH, 
4HDYNP.5HXMAX5, 

4HAMYB, 5HTTMYB, 4hf«ib, / 

IGUID( I ) - 1,/ INERTIAL EULER ANGLE OPTION 
1GUID(2) - 0. / SAME STEERING OPTION FOR ALL CHANNELS 
IGU1D<4) - 1. / CUBIC POLY WITH CONSTANT TERM SET BY INPUT 
PITPC( 1 ) - 0 , / INITIAL PITCH RATE 
YAWPC(l) - 0 , / INITIAL YAW RATE 
ROLPQ1) - 0 , / INITIAL ROLL RATE 

NPC( 1 ) - 2, / PRINT CONIC BLOCK AT EACH PHASE CHANGE 
PINC - 20, / PRINT INTERVAL EVERY 20 SECONDS 
NPC(2) - 1, / FOURTH ORDER RUNGE KUTTA 
DT - l. / INTEGRATION STEP SIZE 

NPC(3) - 4, / PLANET RELATIVE INPUT ON VELOCITY VECTOR 
VELR - 1, / INITIAL VELOCITY AT LAUNCH PAD 
GAMMAR - 90, / INITIAL PLANET RELATIVE FLIGHT PATH ANGLE 
AZVELR - 0 , / INITIAL AZIMUTH ANGLE OF VELR 


■zl -63, 

NPC(4) - 2, / PLANET RELATIVE INPUT ON ALTITUDE VECTOR 
ALTTTO - 0 , / INITIAL ALTITUDE 
GDALT - 0 , / INITIAL ALTITUDE 
GDLAT - 28 5, / LATITUDE OF LAUNCH SITE 
LONG - 280 0. / LONGITUDE OF LAUNCH SITE 
LONG1 -280 0,/ LONGITUDE EAST OF PRIME MERIDIAN 
NPC(5) - 5, / 1976 US STANDARD ATMOSPHERE 
NPC(6) - 0, / NO WINDS 
NPC(7) - I , / ACCELERATION LIMIT SET 
ASMAX -30,/ LIMIT ACCELERATION THROUGHOUT TRAJECTORY 
LREF - 162 6. /BODY LENGTH FROM NOSE TO BASE 
SREF - 3500 7, /REFERENCE AREA FOR ORB ITER ONLY 
NPC<8) -2, /CL, CD A CM TABLE INPUTS 

NPC<9) - 1 , / ROCKET ENGINE WITH THRUST TABLE AND ISP VAC 
NPC( 10) - 0, / STATIC TRIM IN PITCH 

GXP(1) - 160 2, /X- LOCATION OF ENGINE GIMBAL MEASURED FROM 
NOSE 

(TRIM - 1 . / TRIM WITH ENGINE DEFLECTIONS ONLY 
IENGT - 1 , / CALCULATE ENGINE INCIDENCE ANGLES FROM ST A TIC 
TRIM EQUATIONS 

NPC( 1 1 ) - 0, / NO FUNCTIONAL INEQUALITY CONSTRAINTS 
NPC( 1 2) - 3, / CALCULATE CROSSRANGE AND DOWNRANGE BASED 
ON ORBIT 

ALTREF - 100 0, / 100NMI CIRCULAR MILE REFERENCE ORBIT 
AZREF - 0 , / AZIMUTH REFERENCE 
NPC( 1 3 ) - 0, / DO NOT JETTISON PROPELLANT 
NPC( 14) - 0, / NO HOLDDOWN BEFORE LIFTOFF 
NPC( 15)- 0,/DO NOT CALCULATE AEROHEATING 
NPC( 1 6) - 0, / USE OBLATE EARTH GRAVITY MODEL 
c NPC( 1 7) - 2, / JETTISON BOOSTER 
NPC( 1 8) - 0. / DO NOT TERMINATE TRAJECTORY 
C NPC( 19) - 0, / DO NOT PRINT INPUT CONDITION SUMMARIES 
NPC(20) - 0. / DO NOT USE ANY SPECIAL DT CALCULATION 
NPC<2 1 ) - 0, / DO NOT CALCULATE FUEL AND OXIDIZER WEIGHTS 
AND VOLUMES 

NPC(22) - 0, / DO NOT CALCULATE THROTTLING 
NPC(23) - 0, / DO NOT COMPUTE VELOCITY MARGINS 
NPC( 24) - 0, / DO NOT COMPUTE ANY PARAMETER INTEGRALS 
NPC(25) - 3, / COMPUTE VELOCITY LOSSES AND PRINT AT EACH 
TIME 

NPC(26) - 0, / NO SPECIAL AEROHEATING CALCULATIONS 
NPC(27) - 0, / DO NOT INTEGRATE ENGINE FLOWRATES 
IWPF( 1) - », / INTEGRATE BOOSTER FLOW RATES 
NPC(28) - 0, / TRACKING STATION OPTION NOT USED 
NPC<29) - 0, / DO NOT COMPUTE VACUUM IMPACT POINTS 
NPC<30) - 0, / USE N-STAGE VEHICLE WEIGHT MODEL 
IENGA(1)- 1, /THROTTLE ORBITER 
IENGMF(1)- 1, / ORBITER ROCKETS ARE ON 
IWDF( 1 ) - 1 , / CALCULATE ORBITER FLOW RATES ASA TABLE 
LOOKUP 

NENG - I . / THE NUMBER OF THRUSTING ENGINES 
NPC0 1 ) - 0, / NO VERNAL EQUINOX, SUN SHADOW, SUN ANGLE 
CALCULATIONS 

NPC(32) - 0. / NO PARACHUTE DRAG 

NPC(33) - 0, / NO ARC LENGTH CALCULATIONS 

NPC(34) - 0, / NO KEPLERIAN STATE CALCULATIONS 

NPC(35) - 0, / NO SENSED VELOCITY INCREMENT CALCULATIONS 

NPC(36) - 0, / DO NOT ACTIVATE SUNLIGHT OPTIONS 

NPC(37) - 0, i DO NOT ACTIVATE DATE OPTION 

NPC(38) - 0, / NO ATMOSPHERIC TURBULENCE 

NPC<39) - 0, / NO ATMOSPHERIC GUSTS 

NPC(40) - I,/ 

s 

PSTBLMLT 

TVCIM - 3.3351, / Thrnet multiplier for Orbiter 
WDIM - 3 335 1, /Flow rate multiplier for Orbiter 
AEIM - 3 3351, /Exit area multiplier for Orbiter 
S 

PSTAB 

TABLE - 5HTVC1T, 0,531500 , / VACUUM THRUST OF Orbiter engine 
S 

PSTAB 

TABLE - 4HWDIT.0,! 199 77, / FLOW RATE OF Orbiter engine 
S 

PSTAB 

TABLE - 4HAE1T. 0,33 558, / EXIT AREA OF Orbiter engine 
S 

PSTAB TABLE - 4HXCGT, 1 ,6HWEIGHT,2. 1 ,- 1 . 1 , 

2006450,123 51,134870,1 15 50, 

S 

PSTAB T ABLE- 5HXREFT.0, 109 90, 

S 

PSTAB TABLE-5HZREFT.0.0, 

'include '/hofne/asdl l^tooley/poet/B!MESE/PT2PT«bimeee aero' 

ENDPHS- 1, 

S 

ISgendat 

event- 1 5 ,cri tr-4hti me, value-07 , 
iguid(4)-2, 

S 

IStblmh 

S 
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lStabtabie-4hpitU.6htime ,1.1,1. 1. 

10, 90, /all inertial pitch angles arc ind vara 
30, 90. 

50, 90, 

100. 90, 

ISO, 90, 

200. 90. 

2*0,90, 

320.90, 

«ndphr-l, 

S 

ISgendat 

title - Oh* Ramp up alpha*, 

event- II, critr-5hwprop, valua-*94 1 600, /Maximum allowable propellant 

iengmfl; I )-0, Zaira all engines off 

igmd-O.O.J, /use aero angle* with dalpha aet in next phase 

dalpha - 40, /independent variable 

endpht - 1, 

S 

Uganda! 

tide - Oh* alpha reaches specified value*. 

•vent- 1 9,critr-6hgammai,value-0, 
npc( 1 7>-2, / weight jettison option 

wjen-4925, /ventreeiduals 

endphs - 1, 

S 

ISgcndai 

title - Oh* Begin Decent*, 

•vent- 20, critr-dhgammai .valued), 
dt -4, /increase step siae during long descent 
iguid-0,0.2. /use aero angle tables 
npc(l7)-2, / weight jettison option 

wjett-4295, f vent residuals 

xmax(5)-0, / met max gamma monitor 

S 

IStMmlt 

S 

I Stab tabk-6halphaU ,4hvelr, 10, 1 .- 1 . 1 , 

25000, 40, 

1*000. 40, /Hypersonic trim at 40 degrees 
16000, 0, /control imparled to optimizer 

12000 , 0 , 

10000. 0. 

*000. 0. 

6000. 0, 

4000, 0, 

2000, 0, 

0 , 0 , 

endphs- 1, 

S 

ISgcndai 

title - Oh* Deed nab on scheived*. 

•vent- 500, critr-Shgdalt, value- 50000, 

ENDPHS - 1. 

ENDJOB- 1. 

ENDPRB * 1. 

$ 
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A.6. FUEL-AUGMENTED BIMESE: POINT-TO-POINT 


PSSEARCH 

c*** ••••••• ••••••••• 

C LOX/LH2 BIMESE with fud augmentation W/RS2 100 Engine (revi Bed 
1/1/99) 

C WB-003-C Bimese point to point 
C angle Orbiter cannot make it to orbit 
C Goal is to maximize crotsrenge of orbiter 
C Simulates the ascent (but not to orbit) and reentry 

C The Gliding entry phaae is simulated down to 10,000 ft and a ret vd of 2500 
ft/s 

C At this point energy management phaae would take over, but it is not 
simulated 

C Energy management phase usually results in a 100 nmi crossrange 
C LATEST WB-001 AERO ADDED 5/6/97 



c 

SRCHM - 4, / PROJECTED GRADIENT METHOD 
c IPRO - -I, / PRINT FINAL TRAJECTORY ONLY 
lOfLAG - 0, / ENGLISH IN - ENGLISH OUT 
CONEPS<!) - 89 9, / CONVERGENCE TOLERANCE 
CONEPS<2) - 00001 , / MIN ALLOWED PERCENTAGE CHANGE 
CONEPS<3)-.0000l,/ 

CONEPS<4) -.00001./ 

CONEPS<5) - .00001 , / 

PCTCC-.001. / 

IDEB-0, /No Detail 
C 



C OPTIMIZATION VARIABLE 

* •*•••••• 

C 

OPTVAR - 6HDPRNGI, / MAXIMIZE DISTANCE TRAVELED 
OPT - 1, / MAXIMIZATION FLAG 

OPTPH - 500, / PHASE AT WHICH WEIGHT IS MAXIMIZED 
WOPT - 000003, / OPTIMIZATION VARIABLE WEIGHTING 
MAXITR -40 , ! MAXIMUM ITERATIONS ALLOWED 
C 



C 

C CONSTRAINT VARIABLES 

C 



C 

NDEPV - 5, / NUMBER OF DEPENDENT VARIABLES USED AS 
CONTROLS 

DEPVR - 5HXMAX3,5HXM AX 1 , 5HXMAX2, 5HXMIN2, 5HXMAX4, 
DEPVAL - 4, 1000., 379165 , -379 165 ,40, /INSERTION ORBIT VALUES + 
CONSTRAINTS 

DEPTL - I, 1. 100, 100, I. /DESIRED ACCURACIES 

DEPPH - 500,500,500,500,500,/ EVENTS 
IDEPVR - l, 1, I. -I, I, 



C 

C 

C CONTROL VARIABLES 

C 



C 

tabl( 1 ) -4hpitt,4hpitt,4hpitt,4hpitt,4hpitt,4hpitt,4hpitt,4hpitt, 
tabl(9) 

-6halphat,6haJ phat,6halphat,6halphai,6haiphat,6haJ phat.6hal phat, 6ha) phtt, 
tably( I )- 1 ,2,3,4,5,6,7.8, 
tably(9^3,4, 5,6,7,8.9,10, 

nindv - 16, /In this deck banking angles are not needed to meet constraints 
indvr - 5htabll l 5htabl2,5htafcl3,5htabl4,5htabl5,5htabl6.5htabi7,5htabl8, 
indvr(9) 

-5htabl9,6htabl I0,6hubt 1 1 ,6htafcl 1 2,6htabl I3,6htabl I4,6htabl 1 5,6hub1 1 6, 
indph - 15.15,15,15,15,15,15,15, 
indph(9) - 20,20,20,20,20,20,20,20, 

PERT - 16*1.0E-4,1, / PERTURB 

u- -8.274 1 1 1 864380E-03.-2 1 33260052939E+0 1 .-2 39766697477SE+0 1 . 

-5 626S5S826370E+OI,-7 361997959382E+01.-8 6222274 1 886 1 E+O I, 

-9 3368 1 695769SE+0 1,-9 87396 1 5 1 8 200E+0 1 , 8 6273 1 04 76026E+0 1 . 
u- -8.274245 1 65 569E-03, -2 1391 1064924 1 E+O I, -2 42082430495 IE+01, 

-5 61 80S8820848E+01 ,-7 3465 1601411 4E+0 I .-8 1 7 1 7125 1 3255E+0 1 , 

-9 36754285369 1 E+O I, -9 86026384620 1 E+O l, 8 62730O997425E+O1, 

/u's that result in opt guidance with only alpha control violates normal force 
const 

u- -8 2742 1 74 1 3973E-03.-2 1 2 1 879899720E+0 1 ,-2 4 1 3003298457E+OI . 
-5.462767091 I47E+01.-I.475214444848E+01.-I I 527 590203 79E+01 , 

-9 2792 1 8 1 68809E+0 1 ,-9 86026384620 1 E+0 1 , 8 65 1 244888 1 45E+0 1 , 

40.5 951 I68900506E+01.4. 1 405491 37941 E+0 1, I 065461 141837E-K)!, 

5 929888 736206E+00, 5 058B32037600E+00. 1 .327631 52423 1E-02.- 
60, 60.60 ,-40, -40, 20. 

u- -8 2742 1741 3973E-03.-2 1 2 1 879899720E+OJ.-2 41 3003298457E+O1, 
-5462767091 147E>01, -8 475214444848E+01, -8 1 527590203 79E+01, 

-9 2 792 1 8 1 68809E+0 1 ,-9 860263 84620 1 E+0 1 . 8 65 1 244888 145E+OI, 
40,5.951 168900506E+0 1, 4 140549 13794 1 E+0 1, I 065461 I41837E+01. 


5 929888736206E+O0, 5 05883203 7600E+00, 1.3 2763 1 52423 IE-02,- 10,-1 0,- 

10,-10,-10,-10, 

/converged solution without constraining final velocity 

u- -8 2742 1 74 1 3973E-03.-2 1 2I879899720E+01.-2 4 1 30032984 5 7E+0 1 . 

-5 462767091 147E+01.-8 4752 1 4444 848E+01, -8 1527590203 79E+0 1, 

-9 2792 1 8 1 68809E+01 , -9 860263846201 E+O 1, 7. 176385530439E+01, 

4 0000085028 1 5E+01, 6 I0I418902565E+O1, 4 258546163 130E+01, 
I.063683205035E+OI, 5 913I29650736E+00, 5 045232 193 202E+00, 

1 32763 1425463E-02, 0 9997745204 59E+O0. 0 98477375 1 125E+00, 

0 976526577I65E+00, 0995I855S2387E+00, 0 000005876 1 46E +01, 
OOOOOOOOOOOOOE+OI, 

/converged with velocity constraint 

u- -8 2742I8664838E-03.-2 1 2234854973 7E+0 1,-2 41 54953 29997E+01, 

-5 532985 169I07E+OI, -8.484468 181437E+01, -8 149072879458E+01, 

-9 22 508 7779792E+0 1 ,-9 849049546882E+0 1 , 

4 072197I86392E+0J, 4 04649924860IE+00, 4 0982 706409 1 E+00. 

4 0279901 76507E+01, 

3 02 79901 76507E+01.2 3279901 76507E+OU 3279901 76507E +01, 2 327990176 
507E+O1, 

/use u's for 1 ,000 ib payload 

u» -8 274279039759E-O3.-2. 1 4208 1 663503E+0 1 , -2 361 1 8237036 1 E+O I . 

-5 6 1 88 5 1 5 1 6036E+0 1 ,-9 697944809427E+0 1 ,-6 8392 5 829 1 462E+0 1 , 

-9 3660371761 18E+0L-9 990424738464E+01, 4 OOOOI33 1 1904E+01, 

4 00000003 1 147E+01, 3 091 193B96437E+01, 3 04 728066799 5 E+O 1, 

2 0474 1 8722764E+01, 2.32084401 1004E+OI, 1.315381 B26792E+0 1, 

0 3175 17262700E+01. 

/use u’s for 1,000 Ib payload 

u- -8 2 742854581 16E-03,-2.l44825875449E+Ol,-2 35907901 1735E+01, 

-5.5 196497 1 8889E+01 , -9 249077673 248E+0 1.-7 1 1 6956274567E+01, 
-9.258437443429E+O1.-9.990424738464E+O1, 3 91 0694 12 5902 E+O 1. 

3 927560866427E+OI, 3 048I726O4214E+OI. 3 015810000807E+0!, 

2 035214171 299E+0 1 , 2.3096946 1 7097E+0 1 .1.31 7472828524E+0 1 . 

3 1 76752890923E+O0, 

u- -8 274285001 159E-03, -2 133740974347E+OI, -2 37238262048 1 E+O I, - 
5 502854357052E+01, -9 204608000024 E+O I, 

-7 12744937829SE+01, *9.2601 71928367E +01, -9 99042473B464E+01 , 

3 910694125902E+O1, 3 92 756086642 7E+01. 

3 0481 726042 14E+01. 3 015810000807E+01, 2 0352I4171299E+01, 

2 0096946 17097E+01, 1.3I7472828524E+01. 

3 1 7 67 52890923 E+00, 

u- -8 27428465007 IE-03. -2 427521889897E+01. -2 5744043603 16E +01. - 
5 500979 194040E+01, -9 202426I47050E+01, 

-7 I 28267284684E+OI, -9 260263474492E+01, -9 9904 2 473 8464 E+O 1, 

3 9I0694125902E+OI, 3 927560866427E+0 1 , 

3 048 1726042 1 4E+01, 3 0I58I0000807E+0!, 2 035214171299E+01, 

2 0096946 1709 7E+0 1, 1 3 1 7472828524E+0 1 , 

3 1 7675289O923E+O0, 

u- -8 2 742 774053 52E-03 . -2 2034655524O5E+O1, -2 646959148998E+OI. - 
5 486751359I44E+01, -9 I555101092S5E+O1, 

-7 1 30735674352E+O1 , -9.2 1 124474208 IE+01 . -9 990424 73 8464 E+O 1 , 

3 9046982976 19E+01, 3 914151458512E+01, 

3 039242625757E+01 , 3 009I99393551E+01. 2 032308900060E+01, 

2 007854985308E+0 1 , 131 7695 197695E+0 1 , 

3 1 769307894 16E+00, 

u« -8 274292444634E-03, -2.213526273097E+OI. -2 655677633584E+01, - 
5 5 208962422 70E+0 1,-8 704952345 1 46E+0 1 , 

-7 487853389 145E+01, -8 694521837372E+OI, -9 9904 24 7 3 8464 E+O 1. 

3 79446968 9704E+01, 3 618706367882E+01, 

2 926431 136015E+OI, 2 9427761 1 1155E+01. I 992022902436E+0 1 . 

1 99 12682 10402E+01. 1.32I554806S99E+OI, 

3 1 80024 77 5060E+00, 

u- -8 274302609566E-03, -2 I9BI09731806E+OI. -2 6874523 57529E +01 , - 
5 54532 1538596E+OI. -84759431 1 2996E+0 1, 

-7 36522691 3249E+0 1 . -8 586 19013842 IE+01, -9 9904 2473 8464E+0 1 , 

3 683948762582E+OI, 3 333347 120303E+01, 

2 847403324I08E+01. 2 89096757087SE+OI, 1.97372 171 2864E+0 1, 

1 98877925 1295E+01, 1 32539944671 1E+OI, 

3 18306249909! E+00. 

u- -8. 2743005 30916E-03. -2 193239408865E+01, -2 66692843 383 9E+01, - 
5 97240943 7765E+01, -8 570534750554E+01, 

-8 215 147334882E+OI , -7.889] 39786 1 74E+01, -9.267465301 832E+0 1 , 

3 490164234376E+01, 3 2229I0009984E+01, 

2 82 190656924 IE+01, 2 871 M6523I80E+01. 1 96703682 502SE +01, 

1 987834 1 3703 7E+01, I 3301955345I2E+OI, 

3 1 864894 1 3846E+O0, 

u- -8 274286253070E-03, -2 186 1606798 3 2E+01, -2 669959*67244 E+O I, - 
5 72235I248990E+01, -8 730825579231E+0I, 

-7 564352400382E+01, -8 33573696252 IE+OI, -9 693565430695E+01 , 

3 491010638 185E+01. 3. 223600706 208E+O l . 

2 822052945986E+O1, 2 87I2377I4334E+01, 1 96 70775804 3 7E +01, 

1 987840475437E+01. 1 330172986096E+O1, 

3 1 1647308 5 3 84E +00. 

u- -8 2743005309 1 6E-03, -2 193239408865E+01, -2 66692843 3 839E +01 , - 
5 97240943 7765E+01, -8 570534750554E+01. 

-8 215147334I82E+01, -7 IB9I39786I74E+0I. -9 26746530 1832E+01, 

3 490I64234376E+01, 3 222910009984E+OI, 

2 82 190656924 IE+OI. 2 B71II6523180E+01, I 96703682 502 8E +01, 
l 987834 137037E+0 1, I 330I95534512E+0I, 

3 1 864894 13B46E+O0, 
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U- 4 27430080763SE-03, -2 1 9337*4*2791 E+01. -2 666884420 78 7E+0 1, - 
5 9SO62OO48O27E+0I. 4 584646483 286E +01. 

4 23918874351 1 E+0 1, -7 49061 *1370006+01, -9 2SB6I7718327E+0I, 

3 490164234376E+-01, 3 2229I0009984E+0I. 

2 82 190656924 IE+01, 2 87 1116523 180E+01. 1 96703682 S028E+01. 

1 9I7834137037E+OI. 1 3301955345 1 2E+01, 

3 1 864894 13846E+00. 

u- 4 274301 42791 4E-03, -2 19367767 1 608E+01. -2.66585463 7207E+O1. - 

5 95469S99481IE+01, -8.5181 15781 J98E+OI, 

4 22193975 1839E+0 1 . -7 89871938I359E+01, -9 26 14 7626 1 244E +01. 

3.4863 56778 1 67E+01, 3.2 198101 54606E+01. 

2 821252839576E+01, 2.B70574049722E+01. 1 966858499644E+01. 

1 987805852482E+01 , 1330297094262E+01, 

3 1 86562801 210E+00, 

u- 4.274329468 128E-03. -2 .207572 164824E+01, -2 6480825 3832 7E+01 , - 

6 I24029854938E+01, -7 94926 1935 1 83E +01. 

4.65408689 1647E+01, -7 85740145351 5E+01. 4 8 1830731 95 75E+01. 

3 4697861 332 16E+01, 3 2064 1 0790 1 20E+0 1 . 

2 8 18501 554 105E+01, 2 86824673 31 62E+01, 1 96608404072 1E+01, 

1 987688374684E+01. 1.330730628828E+01, 

3 18687623 7782E+00, 

u- 4 27433308008 1 E-03. -2.2088933251 19E+01. -2 643 1938372 10E+01, - 
6 0976571 1 9546E+0 1 , -7 86807 1 87206 1E+01. 

4 660725 274903 E+0 1 . -7.992I94878060E+01. 4 83801 1470419E+01. 
3.427335034 1 89E+01, 3 172440837729E+01. 

2.811822072577E+01, 2 862380238584E+01. I 964 171 6328 1 8E+0 1 . 

1 987402867 145E+01, l 331822394620E+01. 

3 18766405 1246E+O0. 

u- 4.274329468 128E-03, -2.207572 164824E+01. -2 6480825383 27E+0 1 . - 
6 1 24029854938E+0 1 . -7 94926 1 93 5 1 83E+0 1 . 

4 65408689 1647E+01. -7 857401453 51 5E+01, 4 81 83073 19575E+01. 

3 469786 13321 6E +01. 3 2064I0790120E+01. 

2 8IB501554105E+OI. 2 868246733 1 62E+01, I 966084040721 E+0 1. 

1 9876883 746I4E+0 1 . 1 330730628828E+01, 

3 18687623 7782E+00, 

/63i chump 

u- 4 274333080081 E-03. -2 2088933251 19E+01. -2,643 1 93 8372 10E +01, - 
6 0976571 19546E+01, -7 86807 187206 1 E+0 1. 

4 6607252749O3E+OI. -7 992194878060E+OI, 4 83801 1470419E+01. 

3 427335034189E+01. 3 172440837729E+01. 

2 81 1 822072 577E+01, 2 86238023 8 584 E+0 1. I 964 1716328 18E+01. 

I 987402867 145E+01, l 331822394620E+01. 

3 1 8766405 1246E+00. 

u- 4.2743336 1 3270E-03. -2.202585497 121 E+0 1. -2 645259999024E+01, - 
6 063 3 76525407E+0 1 7 8048444608 1 6E+0 1 , 

4 62359873007 IE+01, 4 0071 98328577E+01. 4 83801 14704 19E+01. 

3 419581771 260E+OI. 3 1642I2737107E+0I. 

2 8 10181 328B56E+0I, 2 86098 5929624E+01. 1 963 703285971 E+0 1, 

1 987333933870E+01, 1 33 20499292 77E+0 1, 

3 18782204493 7E+00. 

u- 4 27433308008 IE-03. -2 2088933251 19E+01. -2 643193837210E+01. - 
6 0976571 19546E+01, -7 86807 1 87206 1 E+0 1 . 

4 66072 52 7 4903 E+0 1, -7 992I94878060E+01. 4 83801 14704I9E+01. 

3 427335034 1 89E+0 1 , 3 1 7244083 7729E +0 1 . 

2 81 1822072577E+01. 2 862 3 802 38 584 E +01. 1 964I71632818E+OI. 

1 987402867 145E+01, I 33 1 822394620E+01 . 

3 1 8766405 1246E+00. 

u- 4 2743336 13270E-03. -2 202585497121E+01. -2 645259999024E+OI. - 
6 063376525407E+01, -7 8048444608 16E +01. 

4 623 59873007 1 E+0 1. 4 007198328577E+01. 4 83801 1470419E+01, 

3 4 19581771 260E+01. 3 1642 12737 1 07E+O1, 

2 810181 328856E+01, 2 S60985929624E+01, 1 963703285971 E+0 1. 

1 987333933870E+0I, 1 3 3 20499292 77E +01. 

3 1 8782204493 7E+00, 

/3640 nmi 

u- 4 27433361 32 70E-03, -2.20258549712 IE+01, -2 645259999024E+01. - 
6.063376 52 5407E+01 , -7 8048444608 16E +01 . 

4 623598730071 E+0 1. 4 007I98328577E+01, 4 83801 1470419E+01. 

3 4I958177I260E+01. 3 1642 12737 107E+01, 

2.821252839576E+01, 2 870574049722E+01. 2 966858499644E+0 1 , 

2 987805852482E+01, 2 330297094262E+01. 

3 186562801 2 10E+00. 

u- 4 274329468 128E-03. -2 207572 164824E+01. -2 6480825 383 27E+01. - 
6 1 24029854938E+0 1 . -7 94926 1 93 5 1 83 E+0 1 , 

4 65408689 1647E+01. -7 857401 45351 5E+01, 4 8 183073 19575E+01, 

3 4697861 332 16E+01, 3 .2064 10790 1 20E+O1, 

2.818501 554 105E+01. 2.868246733 162E +01. I 96608404072 IE+01 . 

1 9876883 74684E+0 1, 1 330730628828E+01 , 

3 18687623 7782E+00, 

u- 4 2743 15061 53 3E-03. -2 083199291 760E+01, -2670226024393E+01. - 
6.040374 503005E+0), -7.667961 1 17423E+01, 

4 41 2923823 1 29E +01. -7 99750005 193 7E +01. 4 8666595944 50E+01, 

3 469786 1 3 32 16E +01, 3 20641 0790 120E+0 1 . 

2 818501554 105E+01, 2 868246733 162E +01. 1 96608404072 IE+01, 

1 98768I374684E+01, I 330730628828E+0 1 . 

3 18687623 7782E+00. 

u- 4 274306536628E-03. -2 01961 827$ 164E+01. -2 68157780I909E+01, - 
6 011 972982 534E+0 1 . -7 609800789898E+0 1 . 

4 386786285 560E+01, -7 997954289628 E+0 1. 4 86622954 521 5E +01. 

3 4685264223 58E+01, 3 205 3 94268403 E+0 1. 

2 81 8293957303E+O1 , 2 868071 OI8069E +01. 1 96602528762 5E+01. 

1 987679537281 E+0 1, 1 330757 147006E+01. 


3 1 16895964 76SE+00, 

Amt yet 

u- 4. 27430963 7309E-03. -2 00057970461 34E+01, -2 67667972 5 598E+0 1 . - 
6 0223 13281 598E+01. -7 44496246 7778E+01. 

4 421749631901E+01. 4 OO4018S14OI6E+OI, 4 869629766 1 4 7E+0 1 . 

3 4697861 3321 6E+0 1, 3 30641 0790 120E+0 1 . 

2 818501 554I05E+OI, 2 .868246733 1 62E+01. 1 76608404072 IE+01. 

1 9876883 74684E+01, 1 330730628828E+0 1 . 

3 1 8687623 7782E+00, 

Amt yet 

u- 4 2743098612I8E-03.-2000721394000E+01, -2 676680228 124E +01. - 
6 018145869216E+O1, -7 43225276001 2E+OI, 

4 3990U537991E+01, -7.9877233 13527E+01. 4 8604522721 71E+0 1, 

3 468987 1981 99E+01. 3 305758391 144E+01, 

2 818356164077E+01. 2.868 12945987 IE+01. 2 066066004 205 E +01, 

1 98768251 5094E+01. I.330747760206E+01, 

3 1 86888963 287E+00, 

Amt yet 

u- 4 2743 209971 54 E-03, -2 00769443 5983 E+0 1. -2 .67570651 1 402E+01, - 

5 85933101 7995E+01, -7216723574025E+01, 

4 3 12579688950E+OI, 4 1204568 594 72E+01. 4 7899388875 1 2E+01, 

3 4313393B9942E+01. 3.275264592107E+OI, 

2 81 17441391 13E+01, 2 86302 1528 797E+01. 2 06 1248469 139E+01. 

1 985550239908E+01, 1.33 153698723 3 E+0 1, 

3 187476745 181E+00, 

$ 

*•* •••*•••••• 

C 

c TRAJECTORY SIMULATION INPUTS 

C 

€•••• 

PJGENDAT 

TITLE - OH* LOX/LH2 BIMESE 28 5 deg ddivenn* 60 Klb *, 

EVENT - 10. / FIRST EVENT 
WGTSG - 1362398. 

FESN - 500, / FINAL EVENT NUMBER 
MAXT1M - 3000. / MAXIMUM TIME 
ALTMAX - 10000000 , / MAXIMUM ALTITUDE ALLOWED 
ALTMIN - -5000 . / MINIMUM ALTITUDE ALLOWED 
PRNC - 0, /setup plot file 
c PRNCA “ 0, /setup delimited plot file 
MONX - 4HDYNP, 4HFAZB, 4HASMG. 5 H ALPHA, 6HGAMMA1. i 
MONITOR Q A NORMAL FORCE 

PRNT(97) - 5HXMAXI. 5HXMAX2, 5HXMIN2, 6HALTTTO, 4HMACR 
4HDYNP.5HXMAX5, 

4HAMYB, 5HTTMYB, 4hfezb. / 

K3U1D(1) - 1./ INERTIAL EULER ANGLE OPTION 
IGUID(2) - 0. / SAME STEERING OPTION FOR ALL CHANNELS 
1GUUX4) - 1. / CUBIC POLY WITH CONSTANT TERM SET BY INPUT 
PITPQ 1) - 0 . / INITIAL PITCH RATE 
Y A WPC( I ) - 0 . / INITIAL YAW RATE 
ROLPC( I ) - 0 , / INITIAL ROLL RATE 

NPC( 1 ) - 2. / PRINT CONIC BLOCK AT EACH PHASE CHANGE 
PINC - 20. / PRINT INTERVAL EVERY 20 SECONDS 
NPC(2) - I, / FOURTH ORDER RUNGE KUTTA 
DT - I, / INTEGRATION STEP SIZE 

NPC(3) - 4. / PLANET RELATIVE INPUT ON VELOCITY VECTOR 
VELR - .1. / INITIAL VELOCITY AT LAUNCH PAD 
GAMMAR - 90 . / INITIAL PLANET RELATIVE FLIGHT PATH ANGLE 
AZVELR - 0 . / INITIAL AZIMUTH ANGLE OF VELR 
tzl -67. 

NPQ4) - 2, / PLANET RELATIVE INPUT ON ALTITUDE VECTOR 
ALTITO - 0 . / INITIAL ALTITUDE 
GDALT - 0 , / INITIAL ALTITUDE 
GDLAT - 28 5, / LATITUDE OF LAUNCH SITE 
LONG - 280 0, / LONGITUDE OF LAUNCH SITE 
LONG! - 280 0, / LONGITUDE EAST OF PRIME MERIDIAN 
NPC<5) - 5, / 1976 US STANDARD ATMOSPHERE 
NPQ6) - 0, / NO WINDS 
NPC(7) - I. / ACCELERATION LIMU SET 
ASMAX -30,/ LIMIT ACCELERATION THROUGHOUT TRAJECTORY 
LREF - 162 6, /BODY LENGTH FROM NOSE TO BASE 
SREF - 3500 7, /REFERENCE AREA FOR ORB ITER ONLY 
NPC(8) - 2. / CL, CD Jk CM TABLE INPUTS 

NPC(9) - 1 . / ROCKET ENGINE WITH THRUST TABLE AND ISP VAC 
NPQ 10) - 0, / STATIC TRIM IN PITCH 

GXP(1) - 160 2, /X -LOCATION OF ENGINE GIMBAL MEASURED FROM 
NOSE 

ITRIM - 1. /TRIM WITH ENGINE DEFLECTIONS ONLY 
IENGT - 1, / CALCULATE ENGINE INCIDENCE ANGLES FROM STATIC 
TRIM EQUATIONS 

NPC( 1 1 )- 0, / NO FUNCTION AL INEQUALITY CONSTRAINTS 
NPQ12) - 3, /CALCULATE CROSSRANGE AND DOWNRANGE BASED 
ON ORBIT 

ALTREF - 1000, / 100 NMI CIRCULAR MILE REFERENCE ORBIT 
AZREF - 0 . / AZIMUTH REFERENCE 
NPC< 1 3) - 0, / DO NOT JETTISON PROPELLANT 
NPC(14) - 0. / NO HOLDDOWN BEFORE LIFTOFF 
NPQ 1 5) - 0. / DO NOT CALCULATE AEROHEATTNG 
NPQ 16) - 0. / USE OBLATE EARTH GRAVITY MODEL 
c NPC( 1 7) - 2. / JETTISON BOOSTER 
NPC( 18) - 0. / DO NOT TERMINATE TRAJECTORY 
C NPC( 19) - 0. / DO NOT PRINT INPUT CONDITION SUMMARIES 
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NPQ20) - 0, / DO NOT USE ANY SPECIAL DT CALCULATION 
NPC(21) - 0. / DO NOT CALCULATE FUEL AND OXIDIZER WEIGHTS 
AND VOLUMES 

NPQ22) - 0, / DO NOT CALCULATE THROTTLING 
NPC<23) - 0. / DO NOT COMPUTE VELOCITY MARGINS 
NPQ24) - 0. / DO NOT COMPUTE ANY PARAMETER INTEGRALS 
NPC(2S) - 3, / COMPUTE VELOCITY LOSSES AND PRINT AT EACH 
TIME 

NPC<26) - 0, / NO SPECIAL AEROHEATTNG CALCULATIONS 
NPC(27) - 0, / DO NOT INTEGRATE ENGINE FLOWRATES 
IWPF( I ) - I, / INTEGRATE BOOSTER FLOW RATES 
NPC<28) - 0, / TRACKING STATION OPTION NOT USED 
NPC(29) -0,/DO NOT COMPUTE VACUUM IMPACT POINTS 
NPC<30) - 0, / USE N-STAGE VEHICLE WEIGHT MODEL 
IENGA(I>- I./ THROTTLE ORB ITER 
lENGMF(l) - I, / ORB ITER ROCKETS ARE ON 
IWDF(1 )- \J CALCULATE ORB ITER FLOW RATES AS A TABLE 
LOOKUP 

NENG - I, / THE NUMBER OF THRUSTING ENGINES 
NPQ3 1 ) - 0. / NO VERNAL EQUINOX, SUN SHADOW, SUN ANGLE 
CALCULATIONS 

NPQ32) - 0, / NO PARACHUTE DRAG 

NPQ33) - 0, / NO ARC LENGTH CALCULATIONS 

NPC(34) - 0, / NO KEPLER IAN STATE CALCULATIONS 

NPC(35) - 0, / NO SENSED VELOCITY INCREMENT CALCULATIONS 

NPC<36) - 0, / DO NOT ACTIVATE SUNLIGHT OPTIONS 

NPC(37) - 0, / DO NOT ACTIVATE DATE OPTION 

NPC(38) - 0, / NO ATMOSPHERIC TURBULENCE 

NPC<39) - 0, / NO ATMOSPHERIC GUSTS 

NPC(40) - I,/ 

S 

PSTBLMLT 

TVC I M - 3 335 1 , / Thrust multiplier for Orbiter 
WDIM * 3 3351, / Flow rite multiplier for Orbiter 
AE1M “ 3 3351, / Exit »re» multiplier for Orbiter 

S 

PSTAB 

TABLE - 5HTVC1T, 0,53 1 500 , / VACUUM THRUST OF Orbiter engine 
S 

PSTAB 

TABLE - 4HWD I T,0, 1 1 99 77 , 1 FLOW RATE Of Orbiter engine 
S 

PSTAB 

TABLE - 4HAE1T, 0,33 558, / EXIT AREA OF Orbiter engine 
$ 

PSTAB TABLE - 4HXCGT, 1 .6HWEIGHT.2, 1,1,1, 

2006450,123 51,134870,1 15 50, 

S 

PSTAB TABLE- 5HXREFT.0, 1 09 90, 

S 

PSTAB TABLE-5HZ REFT, 0,0, 

•include ‘/home/esdl l/5tooley/post/BIMESE/PT2PT/bimese aero' 

ENDPHS- 1. 

S 

ISgendat 

event- 1 5,critr-4htime,value-07 1 
iguid(4>=2, 

S 

IStblmlt 

$ 

IStabtablc-4hpiRj,6htime ,8,1 ,1,1, 

10, 90, /all inertial pitch angles are ind van 
30, 90, 

50, 90, 

100,90, 

150, 90, 

200, 90, 

280, 90. 

320, 90. 
endpht-1, 

S 

ISgendat 

title - Oh* Ramp up alpha*, 

event- 1 8, critr-5hwprop, value-- 1 166600, /Maximum allowable propellant 

iengmf( I )-0, /turn all engines off 

iguid-0.0.3, /use aero angles with dalpha set in next phase 

dalpha-40, /independent variable 

endphs - 1, 

S 

ISgendat 

title - Oh ’alpha reaches specified value*, 
event- 1 9,cntr-6hgammai ,val ue-0, 
npc( 1 7>-2. /weight jettison option 

wjett-4925, /vent residuals 

endphs - 1, 

S 

ISgendat 

tide - Oh ’Begin Descent* , 
evcnt-20,critr-6hgammai,vaJue-0, 
dt -4, /increase step size during long descent 
iguid-0,0,2, /use aero angle tables 
npc(l7)-2, /weight jettison option 

wjett-4295, / vent residuals 


xmax(5)-0, / met max gamma monitor 

S 

IStblmlt 

S 

I Sub taWe-6haJphaO,4hvdr, 10,1, -1,1, 
25000. 40, 

1 8000, 40, /Hypersonic trim at 40 degrees 
1 6000. 0, /control imparted to optimizer 
12000 . 0 . 

10000, 0, 

8000, 0. 

6000, 0, 

4000, 0, 

2000 , 0 , 

0 , 0 . 
endphs- 1. 

S 

ISgendat 

title - Oh* Dead nation acheived*. 
event-500, criti-Shgdalt, value-50000, 
ENDPHS - I. 

ENDJOB- 1, 

ENDPRB - I, 

S 
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* ••••••••••• 

C LOX/LH2 BIMESE W/RS2 100 Engine (revised 2/1/99) 

C WB-003-C Bimese 
C T/WH .5 / atrip on bunt time“10* 

C to 100 wni 

C 215 DEGREES INCLINATION 
C LATEST WB-OOi AERO ADDED 5/6/97 

* * 

c 

SRCKM - 4, / PROJECTED GRADIENT METHOD 
c IPRO - -I , / PRINT FINAL TRAJECTORY ONLY 
I OF LAG - 0. / ENGLISH IN - ENGLISH OUT 
CONEPSO) - 19 9. / CONVERGENCE TOLERANCE 
CONEPS<2) - 0000 1, / MIN ALLOWED PERCENTAGE CHANGE 
CONEPS<3)- 00001,/ 

CONEPS<4)« 00001./ 

CONEPS<5)« 00001. / 

PCTCC-001./ 

IDEB-0, /No Detail 
C 

* 

c OPTIMIZATION VARIABLE 

******** 

c 

OPTVAR - 6HDPRNG 1 . / MAXIMIZE FINAL BURNOUT WEIGHT 
OPT - I, / MAXIMIZATION FLAG 

OPTPH - 500. / PHASE AT WHICH WEIGHT IS MAXIMIZED 
WOPT - 000003. / OPTIMIZATION VARIABLE WEIGHTING 
MAXJTR - 10, / MAXIMUM ITERATIONS ALLOWED 

C 

"”**'*' 

c 

c CONSTRAINT VARIABLES 

C 



c 

NDEPV - 5, / NUMBER OF DEPENDENT VARIABLES USED AS 
CONTROLS 

DEPVR ■ 5HXMAX3.SHXMAX1 . 5HXMAX2, 5HXMIN2. 5HXMAX4, 
DEPVAL - 4, 1000 .. 379165., -379165 .40, / INSERTION ORBIT VALUES ♦ 
CONSTRAINTS 

DEPTL - I, I. 100, 100. 3, /DESIRED ACCURACIES 

DEPPH - 500, 500,500, 500, 500 y EVENTS 
IDEPVR - 1. I. 1,-1, 1, 

C 



c 

c 

c CONTROL VARIABLES 

C 

* 

c 

tabi( 1 ) MhpiMhp»Mhpitt4hpitt,41ip<tt,4hpitl,4hp«tt,4hpitl,4hp»tt 
tabl(10) 

-6halph*l,6haJphat,6halphat.&he)phat.6haJph*t.6hllphil,6haJphat.6half>hal, 

tably< 1 p* 1 ,2.3,4,4.5,6,7,8, 
tably< 10^3.4,5,6.7,8.9, 30. 

mndv - 20, An this deck banking angles are not needed to meet constraints 
indvr ■ 

5htabl 1 ,5htaW2,5btabl3,5ht»bl4.5hubl5,5htabl6,5htaW7.5htabll.5htabl9, 
indvrtIO) 

-6foabi )0,6htabl 1 1 ,6htabl 1 2,6htab) 1 3,6htab( 1 4,6htabl 1 5,6htabl 1 6,6htabi 1 7,6heta 

pc'. 

indvr{19) - 6hpitpc2,6hpitpc2, 
indph - 15.15.15.15,19, 19.19.t9.19, 
indph(10) - 200.200.200,200.200,200^00.200,14,15.16, 
u- -I 14270571 7169E+-00, -4 5523240337 1 IE-01, *2.7052 19506659E-0I . - 
7 014759690919E-OI, -2 3 15063 5 5 1074 E-01, 

-2 1706244 1410 IE-01,- 1 03 1046299574E -01, -1 722001909I34E-OI, - 
1 2 17415*3 5263E-01, -1 1 5001 5 16041 5 E-0 1 . 

-1 909311515345E-03, -l .655436517730E4)l, 1 6291 12269404 E +01, 
u- -1 41611 0012522E+00, -3736 1729 1 702 5E-0 1,-3 21*9921 5 1 B57E-0L- 
7.294691349305E-01, -2 70111 5735413E-01. 

-2 393 3 594041 TIE-01, -1 0917131 2997 5E-01, -1 17011 20UI76E 4)1 . - 
I 212 566472966E-0 1,-1 91 2 704401 3 01 E-0 1 . 

-1 91036533593 IE-03, -1 7 199453924 16E-01, 1 62741150 1636E+0 1, 
u- -1 379159497410E+00, -5 721219749266E-OI. -4 210 54603 3644E -01, - 
7 203025104 5O4E-0 1 . -2 297361915200E-01. 

-I 934959071 753E-OI, -9 119963122101E-02. -1 501943760030E-01, - 
1 122 15964462 1 E-01, -1 671 S65591504E-01, 

-1 909074 29 7634E-03, -1 .657015 165317E-01. 1 631041929051E+01, 
u- -1 56306026 1213E+O0, -4.204227379936E-01. 4 1*01 11729701 E-0 1. - 
6.322 1003 15I56E-0], -2 1796331 671 76E-01, 

-2 03160961 2770E-01, -1 0649 1 071 3472E-01, -1 7441 1051 6730E-01, - 
1 26149392220 1 E-01 . -2 01271 4659367E-0 1 , 

-1 9I2442257671E-03. -1 13919S474341E-01, 1 629175054744E+01 . 
u- -1 59*11 5017013E+00, -3 966502 5 163 73 E-01, -7 444247 193000E-01, - 
5 04*1251 13924E-01.-1 7096624 5 76 5 5 E-01, 


-M3 5901 1 54233E-0 1 , -1 064996* 1 0070E-0 1 , - 1 763 1 599794*1 E-0 1 . - 
1 24791941 1693E-01, -1 951570735172E-01, 

-1 91 1633*24 151E-03.-1 733406327192E-01. 1 630146 599166E+01, 
u- -l 379159497410E+00, -5 721219749266E-01. -4 2*054603 3644E-0 1 , 

-7 20302 5104 504E-01 . -2 2973619152Q0E-0I. 

-I 934959071 753E-01, -9 119963122101E-02, -I 501943760030E-01, - 
1 122 1 5964462 1 E-01, 

-I671S65591504E-01, 

-l 909074297634E-03.-I 657015165317E-0I, 4 072197I16392E+01. 

4 .04649924160 1E+00, 4 09*27064091 E+OO, 

4 0279901 76507E+01, 

3.0279901 76507E+OU.327990176507E+OU 3279901 76507E+OU 327990176 
507E+0I. 

u- 9 242S01759115E-01, -1.3795291 74790E+00.-5 552979221 5S4E-01, 

3 701639152951E-OI.-2.234290616519E-O1.-1 143221 353295E-01, 
-9.06I212975613E*02,-1.231I95620693E-OI, -9.2603 1793173 IE-02, 

1 923 200134 1 22E-0 1 1 904291597594E-03 ,- 1 .2516156551 56E-0 1 . 
4.043547255215E+0I, 5.105 144151214E-02. I 011 93971 05 56E+01, 
4.070377031 501E+01, 2.79901 7650700E-01. 2 3279901 76507E +01 . 
2.3279901 76507E+01, 2.3279901 76507E+01, 

u- 4 417 1 5*363 725E-0 1 .-9 2042451 1 3696E-0 1 ,4 07 1 037724961E-0 1 . 

3 1779943*3 543E-01 .-2.203605593222E-0 1 1 *1535*42 1 1 99E-01, 
-9.9M943222223E-02.-2 201I1610I115E-01.-9 9762524463 77E-02, 

2 4 707 3 189286 5 E-0 1 ,-l 9036923703 77E-03 ,- 1 2457069 1 461 1 E-0 1 , 
3.910694125902E+01, 

3.927560866427E+OI, 3 94* 1726042 14E+01, 3 715I10000107E+OI, 

3 0352 14171 299E+OI, 2 50969461 7097E +01, L01747212S524E+O1. 

5 176752190923E+00, 

u- 4 25 1 147856104E-01 ,-9 53337O436527E-0I.4 43*372 1 37141E-0 1 . 

-3 1 32763765653E-0 1 ,-2 53 14693 36151E-0 1 .-2 4 1 074 1 277909E-0 1 . 

-2 006457*7545 1 E-0 1,-2 230 1 04441959E-0 1 .-2 003 398 252216E-0 1 , 

2 456209939122E-01.-1 903790332027E-03.-1 246720147 1 56E-01 . 

7 375410999323E+01, 5 5729 1 673 7096E+01. 5 7623 1 7541 649E +01. 

3 7I5110000107E+01, 3 035214I71299E+O1, 2 50969461 709 7E +01. 

I 017472121524E+01, 5 I76752190923E+00, 

u- 4 2537l6364710E-0l,-9 5OI70*212094E-OI,4 4I55I2693240E-01. 

-3 92974364 5432E-0 1 ,-2 72977 1 963170E-0 1 ,-2 709495* 1 2173 E-0 1 , 

-2 *05* 1 093000SE-0! ,-2 92919152*2 50E-0 1 .-2 903 1 3*205 567E-0 1 . 

-2 956346101950E-0 1 .-2 9037*9321 1 39E-03.-1 24671 1 3591 53E-03. 

7 374426557 195E+01, 5 572916737096E+OI, 5 7623 17541649E+01, 

3 7151I0000107E+01, 3 0352141 71 299E+01. 2 50969461 7097E+01. 

1 01 7472121 5 24E+01 , 5 1 7675219092 3 E+OO, 
u= 4 2303 1 4401179E-0 1 ,-9.77613 3296223E-0 1 ,-!.622991$05297E-0 1 , 

-3 9697744 1 2M6E-0 1 ,-2 743121542664E-0 1,-2 71 *929 1 043 72E-0 1 , 

-2 *1 2004636312E-0 1 , -2 93 1 3 *74425 1 4E-0 1 . -2 905 1 1 709*23 7E-0 1 , 

-2 956974209 105E-01.-2 9037956 1 6494E-03.- 1 .2467 1 1 55767 1 E-03, 

7 3467*1 025037E+01, 5 572916737096E+01, 5 7623 17541 649E+01. 

3 715110000107E+01. 3 035214171299E+01, 2 5096946 1 709 7E +01. 

I 01 7472821524E+01, 5 1 767 5 2*90923 E+OO, 
u- -0 1 .-9 795743 774426E-02.4 .6369 1 1 0660 3 0E -02. 

-3 9724 7694 9904E-0 1,-2 744101 1 49 1 7 5 E-0 1,-2 7 1 9519950960E-0 1 . 

-2 * 12442497406E-01, -2 93 14926733 59E-0 1,-2 905259201340E-01, 

-2 95701 97631 32E-0 1,-2 903 79601094 3 E-03, -2 24671 1 572561E-03, 

3 5457611 1 2036E+0 1 . 4 072916737096E+0 1 , 3 0623 1 754 1 649E+0 1 . 

3 715110000107E+01, 3 03521417I299E+OI, 2 5096946 17097E +01, 

1 0 1 7472128524E+01 , 5 1 76752*90923 E+OO, 

u- -1 6*7*1 51365ME-01.-1 71386332M43E-01.-2 6233750263 1 2E-01, 

2 1 42177002561E-0 1 , 2 3332*000 1 M3E-02.-2 1 3 1 926223*61 E-0 1 . 

-5 22906922266 1 E-01 , *5 1042 1 1100240E-0 1 .-7 2013771 73091E-0 1 , 

-4 526091954050E-01.-2 923294513265E-03.-2 241*5 1770100E-03, 

4 0520641 3071 3E+02, 3 95*1627236* 3 E+Ol, 3 737 1266301 63 E+Ol, 

3 715110000437E+01, 3 0352I4171299E+O1, 2 5096946 1709 7E +01. 

1 0 1 747282* 524 E-K) 1 , 5 176752190923E+00, 

u» -12742*545*1 16E-03.-2 14412S175449E+01.-2 35907901 1735E+01,- 
2 35907901 1 735E+01, 

-5 S196497ll**9£+Ol,-9.24907767324*E+OI,-6 916956274567E+01. 

-9 251437443429E+01 ,-9 990424 731464E+0 1 , 3 910694 1 25902E+0 1 , 
3927560166427E+O1, 3.0411 726042 14 E+Ol. 3 0151I0000107E+01, 

2 0352141 71 299E+01, 2.30969461 7097E+O1. I.3I7472121524E+01. 

3 1 76752190923E+00, 

u- 4 274279039759E-03.-2 242011663 503E+01.-2 961 1 123 70361 E+Ol. - 
2 961 11237036IE+OI, 

-5 311151516036E+01.4 997944*09427E+01,-9 03925129 1 462E+0 1, 

-9 1660371761 llE+01,-9 99042473 M64E+01 . 4 0000133 11904E+01, 

4 000000031 147E+01. 3 091 193196437E+01. 3 0472*0667995 E+OL 

2 0474 1 » 722 764E +01 . 2 32014401 1004E+01 . 1.9 15381 126792E+01 . 
0317S17262700E+01, 

u- 4.274123 722437E-03.-2 MS746217263E+O1, -4 534742240M2E+01. 

-5 191633542759E+01.-5 91 46250 19722E+0 1.4 540636263795E+OI, 

-9 399519192613E+01.-9 753 162 536695 E+Ol. -9 19042473M64E+01, 

3 70235I9S1527E+01, 3.1135913026I4E+01, 3 0046W231522E+01, 

2 9*1 1 2 1 307979E+0 1 . 2 02340593 1241E+0 1 . 2 299 1 03 797495E+0 1 . 

1 3 1 1M2107970E+01. 3 1 7716059 1 164E+00, 

u- 4 274072572495E-03. -2 S16641501994E+01, -3 34701*1 1 1551E+01. - 
6.1373252624ME+01, -6.251 257715392E+01, 

-9 64091440 34 3 5E+0 1, -9243020614162E+01, 4 5214731 SI277E +01. - 
9 1111611673 1 6E+OI, 3 70 1 363143 799E +01. 

3 1*3 0395 18934E+0 1 . 3 004256333 114E+OL 2 9*79404*75 52E +01. 

2 023272349109E+01. 2 29*990539951 E+Ol, 

1 3I113O959560E+OI. 3 177159559009E+00, 
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u- 4 274071353 USE-03, -2 559443750022E+01, -3 2SI3530I6857E+O1, - 
6 59944 78463 59E+0 1 . -6 192 1992747 1 3E+01, 

-9.48287 1 652585E+0 1 , -9 3884794 1 408 2E+0 1, 4 7693 20) 8094 7E+0 1 . - 
9. It 39307 1 5438E+0 1 . 3 7023436219 1 1E+0 1 , 

3.883586650997E+01, 3 004676 545 896E +01, 2 9881 19751687E+OI, 
2.023404774361 E+0 1. 2 299 1028063 93 E+01, 

1 31W42703764E+O1, 3 1 77860 5 82265 E +00, 

u- 4 274071 2631 63E-03, -2 4786440643 $0E +01, -3 370478499887E+01, - 
6 74586 7697088E+0 1 , -6.22 1448221 52 7E +01, 

•9 49744749 1 014E+01, -9 17931249881 1E+01, 4 54326734 1824E+01, - 
9 8784233 36566E+01, 3 69933 1642590E+OI. 

3.881 88898096 7E+01, 3 003389812349E+O1, 2 9875596 12285E+01, 

2 022992988482E+01 , 2 298 744345 56 7E+0 1, 

1 3 188057 18188E+01, 3 1778 7071 740 7E+O0, 

u- 4 274077804089E-03, -2 55BI38848863E+01, -3.290303563780E+OI, - 
6.7506691401 2 1E+01, -6 230268249569E+01, 

-9 5845 1 7380663 E+0 1, -9 2987691 13329E+01, 4 613704570232E+01, - 
9.88267081 1828E+01, 3.701 6893 56 138E+01, 

3. 883 22 1589361 E+01, 3.00439524 1 854E+01, 2 98 7999772823 E+0 1. 

2 0233I6336693E+OI, 2 .2990281 13317E+01, 

13188 348993 23E+0 1, 3.177859896802E+O0, 

U- 4.274084 1 I5086E-03. -2 870265383497E+01, -3 584 10775745 1 E+0 1 . - 
6.7504605 19355E+01, -6 25I798358462E+01, 

*9.6853075 1 4406E+01, *9. 151669193368E+OI, -8 399830725023E+0] , - 
9.88267081 1828E+01, 3 71 7852965609E+01, 

3 8 597 5 868664 9E +01, 3 040374025 202E +01, 2 9464963 7623 3E +01, 

2.0501 18579078E+01, 2 278636373 149E+01, 

1.3 1 8834899323 E+0 1, 3 1 778 59896802 E+00, 

u- 4 2740488 17827E-03, -2.3223084 13291 E+0 1, -3 75531530I5I9E+0I, - 
6 8 246260 3 005 8 E+0 1 , -6.28923 13 34789E+01 , 

-9.7586721 25897E+01, 4 61 95766453 3 3E+01, 4 00546294 98 59E +01. - 
9 88267081 1828E+01, 3 7 14636055494 E+0 1, 

3 745881691 163E+01, 3 0785 1 4598 76 7E+0 1, 2 8662 156964 98E +01. 

2 089786522557E+O1, 2.237834366289E+OI, 

1 .3 1 8052340387E+01, 3 .1777904383 15E+00, 

u- 4 2 7404763 3424E-03, -2 318557591 186E+01, -3 764258952893E+01, - 

6 9735562 17845E+01, -6 32629391 8536E+01, 

-9.851 192368307E+O1, 4 53 1 469692692E+0 1 , -7 913532432629E+01, - 
9 88267081 1828E+OI, 3 70597050682 IE +01, 

3.7428I9197097E+0I, 3 0760766 5 7770E +01, 2 8652207 1244 7E+0 1 . 
2.089062 1 638 S3 E+0 1, 2 23 72495343 OOE+Ol, 

1.31 79885S8446E+0 1 , 3 1 777848 5981 4E +00. 

u- 4.274046680880E-03, *2 3 1627802 1563 E+0 1. -3 7710151 14126E+01, - 
7.073 52967 1 566E+0 1 . -6 3 52084040 1 78E+0 1 . 

-9 93283961 2 1 1 7E+01 , 4 4829 1 1 2573 1 5E+01 . -7 852625085563E+0 1 , - 
9.88267081 1828E+01, 3 704 22 507742 6E +01, 

3 7428 191 9709 7E +01. 3 076076657770E+01, 2 8652207 1244 7E +01. 
2.089062 1 63853E+01, 2 237249534300E+01, 

1.31 7988558446E+0 1 , 3 1777848598 1 4E +00, 

u- 4.274025 107497E-03, -2 226835 16204 1 E+0 1, *3 36059892726 1 E+0 1 . - 

7 301232341 105E+O1, -6 4 188248 3 8908 E+0 1, 

-1.013325675975E+02, 4 4499600982 14E+0 1 . -7 6S8616776559E+01, - 
9 876270635844E+01, 3 7008 72 5008 76E +01. 

3 74087453595 7E+01, 3 07444 7078 595 E+0 1, 2 864 564604 7 48E +01. 

2 088563378577E+OI , 2. 23 68 5 3046228 E +01. 

1 .3 1 794 52809 12E +01 , 3 1777809 1880 5 E +00,0 5.-1 0834837048 17E+00, - 
6 299342591 1 20E-01, 

PERT - 13*1 OE-4. / PERTURB 

8 

C* 

C 

C TRAJECTORY SIMULATION INPUTS 

C 



PSGENDAT 

TITLE - OH* LOX/LH2 BIMESE 28 5 deg delivering 60 KJb *. 

EVENT - 10, / FIRST EVENT 
FESN - 500, / FINAL EVENT NUMBER 
MAXT1M - 5000 , l MAXIMUM TIME 
ALTMAX - 10000000 , / MAXIMUM ALTITUDE ALLOWED 
ALTMIN - -5000, / MINIMUM ALTITUDE ALLOWED 
PRNC - 0, /setup plot file 
PRNC - 0, /setup plot Tile 
PRNC A - 0, /letup delimited plot file 

MONX - 4KDYNP, 4KFAZB. 4HASMG. 5H ALPHA. 6HGAMMA1, / 
MONITOR Q A NORMAL FORCE 

PRNT(97) - 5HXMAXI, 5HXMAX2, 5HXMIN2. 6HALTTTO, 4HMACH, 
4HDYNP.3HWD 1 ,3HWD2, 5 h wpcp 1 , 5h wprp2, 

4HAMYB, 5HTTMYB, 5HPSTOP, / 

K3U1D(1) - I, / INERTIAL EULER ANGLE OPTION 
IGU1D(2) - 0. / SAME STEERING OPTION FOR ALL CHANNELS 
IGUID(4) - 1, / CUBIC POLY WITH CONSTANT TERM SET BY INPUT 
PITPC( 1 ) - 0 , / INITIAL PITCH RATE 
YAWPC( I ) - 0 , / INITIAL YAW RATE 
ROLPQ I ) - 0, / INITIAL ROLL RATE 

NPC(I) - 2, / PRINT CONIC BLOCK AT EACH PHASE CHANGE 
PINC - 20, / PRINT INTERVAL EVERY 20 SECONDS 
NPC(2) - l , / FOURTH ORDER RUNGE KUTTA 
DT - I, / INTEGRATION STEP SIZE 

NPC(3) - 4, / PLANET RELATIVE INPUT ON VELOCITY VECTOR 
VELR - 1 , / INITIAL VELOCITY AT LAUNCH PAD 
GAMMAR - 90 , / INITIAL PLANET RELATIVE FLIGHT PATH ANGLE 
AZVELR - 0 . / INITIAL AZIMUTH ANGLE OF VELR 


azl - 90, 

NPC<4) - 2, / PLANET RELATIVE INPUT ON ALTITUDE VECTOR 
ALTITO - 0 , / INITIAL ALTITUDE 
GDALT - 0 , / INITIAL ALTITUDE 
GDLAT - 28 5, / LATITUDE OF LAUNCH SITE 
LONG - 280 0, / LONGITUDE OF LAUNCH SITE 
LONG1 - 280 0, / LONGITUDE EAST OF PRIME MERIDIAN 
NPC(5) - 5, / 1976 US STANDARD ATMOSPHERE 
NPC(6) - 0. /NO WINDS 
NPC(7> - 1 , / ACCELERATION LIMIT SET 

ASMAX -30./ LIMIT ACCELERATION THROUGHOUT TRAJECTORY 
LREF - 162 6. /BODY LENGTH FROM NOSE TO BASE 
SREF - 3500 7, /REFERENCE AREA FOR ORB ITER ONLY 
NPC(8) - 2. /CL. CD* CM TABLE INPUTS 

NPC{9) - 1, / ROCKET ENGINE WITH THRUST TABLE AND ISPVAC 
NPC( 10) - 0, / STATIC TRIM IN PITCH 

GXP( 1)- 160 2, /X -LOCATION OF ENGINE GIMBAL MEASURED FROM 
NOSE 

ITRIM - I, / TRIM WITH ENGINE DEFLECTIONS ONLY 
IENGT - 1, / CALCULATE ENGINE INCIDENCE ANGLES FROM STATIC 
TRIM EQUATIONS 

NPC( 1 1 ) - 0, / NO FUNCTIONAL INEQUALITY CONSTRAINTS 
NPQ12) - 3, / CALCULATE CROSSRANGE AND DOWNRANGE BASED 
ON ORBIT 

ALTREF - 100 0./ 100 NM1 CIRCULAR MILE REFERENCE ORBIT 
AZREF - 0 , / AZIMUTH REFERENCE 
NPC( 1 3) - 0, / DO NOT JETTISON PROPELLANT 
NPC( 14) - 0, / NO HOLDDOWN BEFORE LIFTOFF 
NPC( 1 5) - 0. / DO NOT CALCULATE AEROHEATING 
NPC( 16) - 0, / USE OBLATE EARTH GRAVITY MODEL 
NPC( 1 7) - 2, / JETTISON BOOSTER 
NPC( 1 8) - 0, / DO NOT TERMINATE TRAJECTORY 
C NPC( 1 9) - 0, / DO NOT PRINT INPUT CONDITION SUMMARIES 
NPC(20) - 0, / DO NOT USE ANY SPECIAL DT CALCULATION 
NPC(2 1 ) - 0, / DO NOT CALCULATE FUEL AND OXIDIZER WEIGHTS 
AND VOLUMES 

NPC(22) - 0. / DO NOT CALCULATE THROTTLING 
NPC(23) - 0, / DO NOT COMPUTE VELOCITY MARGINS 
NPC(24) - 0, / DO NOT COMPUTE ANY PARAMETER INTEGRALS 
NPQ25) - 3, / COMPUTE VELOCITY LOSSES AND PRINT AT EACH 
TIME 

NPC(26) - 0, / NO SPECIAL AEROHEATING CALCULATIONS 
NPC(27) - 0, / DO NOT INTEGRATE ENGINE FLOWRATES 
IWPF( 1 ) - 1,1./ INTEGRATE BOOSTER AND STRAP FLOW RATES 
NPC(28) - 0, / TRACKING STATION OPTION NOT USED 
NPC(29) - 0, / DO NOT COMPUTE VACUUM IMPACT POINTS 
NPC(30) - 3, / USE N-STAGE VEHICLE WEIGHT MODEL 
IENGA( 1 ) - 0, 1 . / THROTTLE ORB ITER, BUT NOT STRAP ONS 
IENGMF(l)- 1,1, /ORB ITER AND STRAP ROCKETS ARE ON 


iwdf - 

1,1, /uae lip vac 

ncngl - 

l. 

nengh - 

2, /highest number of engines 

mentnk - 

1,2, /map each engine to a specific tank 

wprp(l) - 

37648. /weight of SRB propellent 

witpd(l) - 

376480, /dry weight of SRB 

wstpd(2) - 

1181050. 

men ftp - 

L2. 


NPC(3 1 ) - 0, / NO VERNAL EQUINOX SUN SHADOW. SUN ANGLE 
CALCULATIONS 

NPC(32) - 0, / NO PARACHUTE DRAG 

NPC(33) - 0, / NO ARC LENGTH CALCULATIONS 

NPC(34) - 0, / NO KEPLER1AN STATE CALCULATIONS 

NPC(35) - 0. / NO SENSED VELOCITY INCREMENT CALCULATIONS 

NPC<36) - 0. / DO NOT ACTIVATE SUNLIGHT OPTIONS 

NPC(37) - 0, / DO NOT ACTIVATE DATE OPTION 

NPC(38) - 0, / NO ATMOSPHERIC TURBULENCE 

NPC(39) - 0, / NO ATMOSPHERIC GUSTS 

NPC(40) - 1,/ 

S 

PSTBLMLT 

TVC2M - 3 3351, / Thru* multiplier for Orbiter 
WD2M - 3.3351, / Flowrate multiplier for Orbiter 
AE2M - 3 335],/ Exit area multiplier for Orbiter 
CDM -12, /Table multiplier for drag becu&ie of tirapom 
$ 

PSTAB 

TABLE - 5HTVC2T, 0.53 1500 ,/ VACUUM THRUST OF Orbiter engine 
S 

PSTAB 

TABLE - SHTVC IT, 0,1 120300, /VACUUM THRUST OF mp-ons engine 
S 

PSTAB 

TABLE - 4HWD2T.0.1 199 77, / FLOW RATE OF Orbiter engine 
S 

PSTAB 

TABLE - 4 HWD IT, 0,40 70, / FLOW RATE OF atra|«ni mgine 
S 

PSTAB 

TABLE - 4HAE2T, 0,33 551, / EXIT AREA OF Orbiter engine 
S 

PSTAB 

TABLE - 4HAE1T, 0,56 869, / EXIT AREA OF ttrmp-ofw engine 
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PSTAB TABLE - 4HXCGT.I.6HWEIGHT.2,M.l, 

2006450.123 51,134870,115 50, 

S 

PSTAB TABLE-5HXREFT.0.I09 90. 

S 

PSTAB TABLE- 5HZREFT.0,0. 

•indude Thome/wdl l^iooiey/poet^lMESE^T2PT/btmae aero’ 

ENDPHS- I. 

S 

ISgendal 

event-14,l,critr - Shame , 

value -20, /chwtgt throttle to acceleration limiting 
xmax(5) - 0, /neat throttle to can make sureetape(l)iifi1 set to above one 
etapc(l)- 5, 
endphs - 1. 

S 

ISgendal 

event- 1 $.cntr-»hd me, valueH)7, endphs- 1 , 
ipiid(4)-0, 

S 

ISgendal 

event- 1 6,critr-4htim«,velue-20.endph*- 1 , 

S 

Uganda! 

event- 1 8, 1 ,critr-4hasmg. 

value -2 999, /change throttle to acceleration limiting 
npo(7) - ), / limit acceleration 
umax -30./ limit acceleration to 3 g*» 
npc<22)-2, ^ 

ienga( 1) - 0,1, / ttoottlc whiter using acc limit, but not atrip one 
etapdl)- 1, 
endphs - I, 

S 

ISgendal 

event- 19, 1 .critr-Shwprp 1 .valueH), 

value - 0, /turn on new engine when 1 * stage fuel is empty 
netpi -2, 

iengmfll) *0, /strap ons off 

npc(7)- I./ limit acceleration 
umax - 3 0. / limit acceleration to 3 g't 
npc(22)-2, 

ienga(l )- 0,1. /throttle obiter using acc limit 
etapdD- I. 

S 

pUblmlt 
edm * 10. 

S 

l$ub table-4hpitt,l,6htime ,1, 1.1,1. 

10, 90, /all inertial pitch angles are ind van 
30, 90. 

50. 90. 

115, 90, 


150. 90. 

200. 90. 

2*0. 90. 

320. 90. 

endphs-i. 

S 

ISgendal 

title — 0h*Ramp up alpha*. 

event- 100, critr-5hwprp2, value--909679. /Maximum allowable propdlani 

iengmfl I )-0, 0. Aura all engines off 

iguid-0,0,3, Aiae aero angles with dalpha aet in next phase 

dalpha - 40, /independent variable 

endphs - 1, 

S 

ISgendal 

title - 0h*alpha reaches specified value*, 
event- 1 50,cfitr-6hgammai,value-0.mdl - 9, 
npc(17)-2, / weight jettison option 

wjett-4925, / vent residuals 

andphs - I. 

S 

ISgendal 

title - Oh* Begin Descent*. 
event-200, critr-dhgammai, value-0, 
dt -4, /increase step sire during long descent 
igmd-0,0.2. /use aero angle tables 
npc(17)-2, /weight jettison option 

wjett-4295, / vent residuals 

xmax(5>-0. / reset max gamma monitor 

$ 

IStbimlt 

S 

I Stab tsWe-6halphat,l.4hvelr. 10,1, -1,1, 

25000. 40, 

18000, 40, /Hypersonic trim at 40 degrees 
1 4000, 0, /control imparted to optimizer 
12000 , 0 , 

10000, 0. 

8000, 0. 

6000. 0, 

4000, 0. 

2000 , 0 , 

0. 0, 
endphs- 1. 

S 

ISgendat 

title - Oh* Destination acheived*. 
event- 500, critr-5hgdali,value-50000, 

ENDPHS - I. 

ENDJOB- l. 

ENDPRB - 1, 

S 
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Point-to-point Bimese with SRM 
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Fast Package Advantage 
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Estimated Market Elasticity 
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Recurring Cost Prediction 



Accounting only for recurring costs , it is clear tnai me 
configurations with strap-on SRB thrust-augmentation 
cannot meet the $1 M/flight price target 



Non-Recurring Cost Summary 
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